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FOREWORD 

This  is  the  Final  Report  covering  work  performed  by  Thiokol  Chemical 
Corporation,  Reaction  Motors  Division,  Denville,  New  Jersey  under  Air 
Force  Contract  AF04(611)-ll6.s0.  The  internal  report  number  is  Reoort 
RMD-5809-F.  The  work  was  administered  under  the  direction  of  Capt. 
William  H.  Summers,  Rocket  Propulsion  Laboratory,  Edwards  Air  Force 
Base,  under  Project  No.  3148. 

The  research  effort  reported  herein  was  conducted  during  the  period 
1  June  1966  to  30  November  1966  under  RMD  Project  b809.  This  work  is 
an  extension  of  previous  work  performed  bv  Thiokol-RMD  under  Contract 
No.  AF04( 611) -9946  and  reported  in  AFRPL-TR-65-257,  dated  February, 
1966. 


The  principal  investigator  of  the  effort  described  herein  is  Mr. 

T,  F.  Seamans.  Other  contribu*orr  to  the  program  are  Dr.  M.  Vanpee 
and  Messrs.  G.  A.  Brown,  B.  E.  Dawson.  F.  G.  Hoffman,  A.  F.  Lum 
and  R.  L.  Storms.  Dr.  Vito  D.  Agosta,  Professor  of  Aerospace 
Engineering,  Polytechnic  Institute  of  Brooklyn,  served  as  a  eonsultant 
principally  for  the  heat  transfer  aspects  of  the  vaporisation  model. 

This  techiucal  report  has  been  reviewed  and  is  approved. 

W.  H.  Ebelke,  Colonel,  USAF 
Chief,  Propellant  Division 
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ABSTRACT 

By  experimental  and  theoretical  means,  a  workable,  fundamental  model 
of  hypergolic  ignition  in  space-ambient  engines  has  been  de.  w loped  by 
refining  an  ea iMac  rncdei  in  three  Important  areas-  The  model  a?pl<»*  to 
propellant  combinations  whose  dominant  ignition  reactions  are  gas  phase 
reactions,  specifically  /hydrazine-type  fuels.  The  model  is  based 

on  both  physical  kinetics  of  propellant  droplet  evaporation  and  overall 
kinetics  of  ignition  reactions.  The  refinements  made  to  the  earlier  model 
are:  (a)  incorporation  of  propellant  flashing  effects  within  the  injector 
manifold  volume  following  valve  opening,  (b)  accounting  for  p re -ignition 
heat  transfer  between  the  thrust  chamber  walls  and  the  vaporizing  vapor/ 
drop  system,  and  (c)  accounting  for  formation  of  a  reaction  intermediate  by 
N2O4/MMH  during  the  ignition  delay  period.  For  two  different  engine 
configurations,  the  resulting  computerized  mode  predicted  pre-ignition 
chamber  pressure  histories  and  ignition  delay  times  that  agree  well  with 
experimental  values.  The  N2O4/MMH  reaction  intermediate,  which  was 
found  to  be  a  mixture  of  at  least  four  compounds  the  principal  ones  of 
which  are  MMH*  H2O  and  MMH*  HNO3,  evolves  heat  upon  its  formation 
and  consequently  causes  markedly  shorter  ignition  delays  than  would 
otherwise  occur.  From  thermochemical  determinations,  it  appears  that 
the  reaction  intermediate  plays  an  important  role  in  severe  ignition 
pressure  spiking  and  possibly  oxidizer  manifold  explosions  as  well. 
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SECTION  I 
INTRODUCTION 


A  major  problem  of  attitude  control  engines  operating  on  nitrogen 
tetroxide  and  hydrazine* type  fuels  is  the  repeatability  and  reliability 
of  the  ignition  of  these  hypergols.  In  high -response  engines  capable  of 
very  short  minimum  impulse  bits,  ignition  pressure  spikes  occur  at 
space  conditions.  The  amplitude  of  these  spikes  can  attain  many 
thousands  of  pounds  per  square  inch  and  cause  structural  failure  of  the 
thrust  chamber.  Thic  is  especially  true  for  those  chambers  fabricated 
of  a  brittle  material  such  as  molybdenum. 

Ignition  pressure  spikes  can  also  trigger  combustion  instability 
which,  if  the  pulse  is  cf  sufficient  duration,  can  likewise  cause  structural 
failure  due  to  excessive  local  heating.  Even  with  short  pulses,  the  impulse 
developed  is  different  with  unstable  combustion.  Consequently,  impulse 
reproducibility  is  compromised,  leading  to  potential  control  difficulties. 

Another  phenomenon  occurring  during  the  ignition  interval  concerns 
explosions  within  the  oxidizer  injector  manifold  which  can  also  lead  to 
failure  of  the  hardware,  especially  the  propellant  valve.  The  genesis  of 
the  oxidizer  injector  explosions  is  believed  to  originate  during  tail-off  of 
the  previous  pulse  and  involves  propellant  vaporization  and  condensation 
and  possibly  the  formation  and  subsequent  rapid  reaction  of  a  reaction 
intermediate  within  the  manifold. 

Hardware  "fixes"  and  propellant  additives  to  alleviate  these  ignition- 
related  problems  have  been  only  partially  effective  (1-3).  A  need  exists, 
therefore,  for  a  mere  fundamental  study  of  ignition  of  these  propellants  at 
reduced  pressures  in  order  to  provide  information  from  which  to  develop 
rational  corrective  and/or  preventive  measures. 

This  report  covers  a  fundamental  study  performed  by  Thiokol 
Chemical  Corporation,  Reaction  Motors  Division,  during  a  six-months 
program  sponsored  by  Edwards-RPL  under  Contract  AF0-l(6ll)-l  1630. 

The  overall  objective  of  the  program  is  to  refine  and  verify  a  fundamental, 
mathematical  model  of  hypergolic  ignition  delay  in  space-ambient  engines. 

The  initial  development  of  the  model,  for  N^O^hydrazine-type 
fuels,  was  performed  at  Thiokol-KMD  under  Contract  AF04(6i  l)-9946. 

Part  II.  Reference  1  is  the  final  report  of  the  six  months  program.  It 
was  found  that  at  low  ambient  pressures  the  dominant  ignition  reactions 
of  N^gO^hydracinc-type  fuels  are  gas  phase  reactions  and  strongly 
pressure  dependent.  It  was  necessary,  therefore,  for  the  model  to  in¬ 
clude  pre -ignition  chamber  pressurisation  due  to  propellant  vaporization, 
as  well  as  chemical  kinetics  of  ignition  reactions. 
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The  resulting  model  consisted  of  two  coupled  sets  of  equations.  In 
one  set,  the  pre-ignition  pressure -time  history  in  a  thrust  chamber  is 
determined  taking  into  account: 

a)  propellant  physical  properties, 

b)  rate  of  propellant  injection, 

c)  physical  kinetics  of  propellant  droplet  evaporation  and  vapor 
condensation, 

d)  heat  transfer  between  propellant  vapor  and  drops, 

e)  drop  residence  time, 

f)  propellant  vapor  efflux  through  the  nozzle,  and 

g)  thrust  chamber  geometry. 

In  the  other  set,  the  chemical  delay  time  as  a  function  of  pressure, 
temperature  and  vapor-phase  composition  is  determined.  The  needed 
overall  kinetic  factors  for  the  vapor-phase  ignition  reactions  of  these 
hypergols  were  experimentally  determined  during  the  earlier  program. 

Verification  of  the  ignition  delay  model  was  precluded  by  two 
unexpected  findings:  appreciable  pre-ignition  heat  transfer  between 
injected  propellants  and  thrust  chamber  walls,  and  formation  of  a  reaction 
intermediate  at  sub-ignition  pressures.  Both  of  these  affect  the  pre-ignition 
pressure  rise  in  a  thrust  chamber  and  thereby  the  ignition  delay  time  as 
well. 


The  objectives  of  »he  present  program  are  to  refine  and  then  verify  a 
computerized  mathematical  procedure  for  calculation  of  hypergolic 
propellant  ignition  delay  in  space  imbient  engines,  including  the  effects  of 
pre-ignition  heat  transfer  and  formation  of  a  sub-ignition  reaction 
intermediate.  The  program  consists  of  three  closely  interrelated  phases 
whose  specific  objectives  are: 

Phase  I  -  Experimental  evaluation  of  pre-ignition  heat  transfer 
between  injected  propellants  and  thrurt  chamber  walls 
for  N^O^/MMH  and  two  thrust  chamber  configuration. 

Phase  II-  Investigation  of  the  effect  of  pre-ignition  formation  of 
the  N^O^/MMH  reaction  intermediate  on  ignition  delay 
and  ignition  pressure  spikes  including  rate  of  formation, 
heat  of  formation,  shock  sensitivity  and  contribution 
to  pre -ignition  pressure  transients. 

Phase  III-  Refinement  of  the  mathematical  ignition  delay  model 

developed  by  previous  AF04{6l  l)-9946  contract  efforts 
to  include  the  effects  of  pre -ignition  beat  transfer  and 
formation  of  a  sub-ignition  reaction  intermediate. 
Verification  of  the  refined  model  by  calculating 
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ignition  delays  for  NgO^/MMH  in  various  different 
thrust  chamber  configurations  and  comparing  the 
results  with  test  values. 

It  has  been  necessary  to  expand  the  Phase  I  effort  to  include  the  effect 
on  chamber  pressurization  of  propellant  "fLashing"  within  the  injector 
volume  at  valve  opening.  The  short  but  finite  times  required  for  the 
propellant  flows  to  reach  nominal  steady  state  levels  slow  the  pre -ignition 
pressure  rise  in  a  thrust  chamber  and  thereby  affect  ignition  delay  times. 

The  results  of  the  theoretical  and  experimental  program  are  discussed 
in  the  following  sections. 
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|  SECTION  II 

SUMMARY 

The  overall  objective  of  the  program  reported  herein  is  to  refine  and 
verify  a  fundamental,  mathematical  model  of  hypergolic  ignition  delay  in 
space-ambient  engine*.  The  model  applies  to  propellant  combinations 
whose  dominant  ignition  reactions  are  gas  phase  reactions,  specifically 
/hydrazine- type  fuels.  The  model,  which  has  been  refined  in  three 
important  areas  in  the  present  program,  is  based  on  both  physical  kinetics 
of  propellant  droplet  evaporation  and  overall  chemical  kinetics  of  gas-phase 
ignition  reactions. 

This  reports  the  final  report  of  a  six-month  program  which  was  per¬ 
formed  in  three  interrelated  phases.  The  purpose  of  Phase  1  wan  to  experi¬ 
mentally  evaluate  pre-ignition  heat  transfer  between  thrust  chamber  walls 
and  the  vaporizing  vaoor/drop  system  in  two  thrust  chamber  configurations. 
The  propellant  combi  ation  of  concern  in  this  phase  and  throughout  the 
program  is  ^O^/MMH.  Phase  I  was  expanded  to  include  evaluation  of 
propellant  "flashing"  effects  within  the  injector  volume  following  valve 
opening.  Both  the  heat  transfer  and  the  initial  propellant  flashing  within 
the  injector  volume  affect  pre-ignition  pressure  rise  in  a  thrust  chamber 
and  thereby  affect  the  ignition  delay  time  as  well. 

Phase  II  was  an  investigation  of  a  reaction  intermediate  formed  by 
N^O^/MMH  at  sub-ignition  conditions.  The  existence  oi  the  reaction  inter¬ 
mediate  had  been  established  previously  bn.  it  remained  for  the  present 
program  to  uitermine  its  compojition,  properties,  stability,  heat  of 
formation  and  rate  of  formation  in  order  to  incorporate  formation  of  the 
material  during  an  ignition  delay  period  into  the  fur  .lament  al  hypergolic 
ignition  model. 

The  purpose  of  Phase  III  was  to  refine  the  mathemaii  al  ignition  delay 
model  developed  by  previous  AF04(6l  l)-994&  contract  efforts  to  include 
pre-ignition  heat  transfer  and  formation  of  the  sub- ignition  reaction 
intermediate.  Propellant  flashing  effects  within  the  injector  volume  were 
included  as  well.  The  refined  tnodel  was  then  to  be  verified  by  comparing 
calculated  and  experimentally  measured  engine  ignition  delays  for  various 
thrust  chamber  configurations. 

Tests  were  conducted  with  five  50 -lb  thrust  ergi  it  during  the  program 
at  a  simulated  altitude  of  18$£00  ft.  The  ranges  of  chamber  design  para¬ 
meters  covered  are:  design  chamber  pressure.  20  to  200  psia:  L*.  10  to  50 
inches;  and  contraction  ratio,  3.  2  to  9.  0.  It  was  plannee  to  use  two  injectors 
during  the  program:  a  one-on-one  mpinging  stream  doublet  injector  for  hot 
Brings  and  a  self-impinging  doublet  injector  for  cold  flow  tests  with  each 
propellant  individually.  The  purpose  of  the  cold  flow  tests,  in  which  chamber 
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pressurization  due  to  propellant  vaporization  was  measured,  was  to  provide 
a  basis  against  which  tc  compare  computed  vaporization  pressurization 
curves  for  each  propellant.  The  self -impinging  injector  was  used  in  order 
to  obtain  the  effects  of  stream  interaction  on  spray  break-up  for  a  singly 
flowed  propellant.  This  injector  maintains  the  same  manifold  volume  and 
liquid  injection  velocity  as  one  propellant  side  of  the  one-on-one  doublet 
injector  used  for  hot  firings.  The  two  injectors  also  have  the  same  im¬ 
pingement  angle  and  impingement  length.  Of  necessity,  however,  the 
internal  surface  area  of  the  self -impinging  injectoi  is  larger  and  the 
propellant  path  is  not  as  direct  as  would  be  desirable. 

Despite  the  similarities  between  the  injectors,  chamber  pressuriza¬ 
tion  due  to  vaporization  of  each  propellant  singly  flowed  proceeded  more 
slowly  with  the  self -impinging  injector  as  a  consequence  of  its  higher 
internal  surface-to-volume  ratio  and  the  more  tortuous  propellant  path. 

The  one-on-one  doublet  injector  was  therefore  used  for  subsequent  cold 
flow  tests  and  for  all  hot  firings. 

Experimental  evidence  of  important  pre-ignition  reactions  was  obtained 
in  two  way 8  from  the  engine  tests  with  the  one-on-one  doublet  injector.  In 
addition  to  chamber  pressure  histories  which  were  measured  in  all  engine 
tests,  chamber  wall  temperature  histories  were  measured  in  several  of  the 
cold  flow  and  hot  firing  tests.  With  regard  to  the  pressure  measurements, 
it  was  found  that  the  pre-ignition  chamber  pressure  of  a  hot  firing  increases 
markedly  more  rapidly  and  to  higher  levels  than  the  simple  sum  of  the 
pressure  histories  of  each  reactant,  singly  flowed.  Tae  chamber  wall 
temperature  histories,  which  were  obtained  with  a  flush  mounted,  high 
’-esponse,  surface  junction  thermocouple,  show  that  the  walls  cool  in  the 
case  of  singly  flowed  propellants  whereas  they  heat  up  during  the  ignition 
delay  times  of  hot  firings. 

The  effects  of  propellant  flashing  within  the  injector  volume  on  pre¬ 
ignition  chamber  pressurization  were  incorporated  into  the  ignition  delay 
model  by  empirically  determining  time -dependent  coefficients  of  injected 
mass  flowrate  for  each  propellant.  For  the  feed  systems /propellant  valves  / 
injector  assembly  used  during  the  program,  a  "best  fit"  was  obtained  with 
times  from  initial  entry  of  propellant  into  the  thrust  chamber  to  attainment 
of  full  steady  state  flowrates  of  1.  6  and  0.  8  msec  for  N2O4  and  MMH,  respec¬ 
tively.  Using  linearly  increasing  mass  flowrates  during  these  time  intervals, 
good  agreement  between  calculated  and  experimental  chamber  pressure 
histories  during  the  flow  build-up  period  was  obtained  for  each  propellant. 

The  predominant  mechanism  for  pre-ignition  heat  transfer  between 
thrust  chamber  walls  and  the  vaporizing  vapor/drep  system  was  determined 
by  analysis  tc  be  conductive  heating  of  propellant  droplets  upon  impingement 
with  chamber  walls.  Other  mechanisms  considered  were  radiant  heating 
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and  convective  heating  of  propellant  vapor.  The  effects  of  conductive 
droplet  heating  were  incorporated  into  the  ignition  delay  model  by 
empirically  determining  time -dependent  coefficients  for  the  time- 
varying  temperatures  of  the  drops  of  each  propellant.  The  time 
dependent  characteristic  of  the  coefficients  is  required  to  account  for 
the  experimentally  determined  fact  that  chamber  wall  surfaces  quickly 
cool  during  vaporization  of  each  propellant.  The  wal .  temperatures 
level-off,  however,  as  the  walls  become  coated  with  condensed  phase 
propellant  and,  thus,  the  coefficients  used  to  account  for  wall  heating 
are  held  independent  of  time  after  being  initially  time -dependent. 

In  practice,  three  coefficients  had  to  be  empirically  evaluated  to 
.account  for  wall  heating.  They  are  related  to  the  maximum  heating  rate 
during  a  run,  the  minimum  (steady  tetate)  rate,  and  the  time  daring  which 
the  rate  decreased  from  the  maximum  to  the  minimum.  Conductive 
heating  of  each  drop  is  permitted  in  the  computer  program  only  after  the 
drop  has  had  time  to  travel  from  the  injector  to  the  chamber  wall.  Values 
of  the  coefficients  giving  good  agreement  between  calculated  and 
experimental  chamber  pressure  histories  were  determined  for  and 

MMH  in  each  of  two  engine  configurations.  It  remained  therefore  to 
incorporate  into  the  ignition  delay  model  the  formation  of  the  N^O^/MMH 
reaction  intermediate  during  the  ignition  delay  period. 

Reaction  intermediates,  or  "residues",  formed  by  N2O4 /hydrazine- 
type  fue’  have  been  observed  under  various  conditions  in  several 
laboratories  and  during  engine  testing,  particularly  at  simulated  altitudes. 
Little  was  known  about  the  material,  and  therefore,  a  number  of  analyses 
and  determinations  were  made  on  the  liquid  reaction  intermediate  which 
forms  from  vapor  N2O4/MMH  at  low  pressure  and  room  temperature. 

The  initial  tests  on  the  material  were  performed  to  determine  its 
stability  and  sensitivity.  The  material  was  found  to  be  less  stable 
thermally  and  to  impact  than  monomethyl  hydrazine  nitrate,  HNO3, 

synthesized  in  the  laboratory  from  MMH  and  HNO3.  The  explosive  strength 
of  the  N2O4/MMH  reaction  intermediate  was  found  to  be  close  to  that  of 
lead  azide. 

The  density  of  the  liquid  material  was  determined  at  three  temperatures. 
Its  density  at  25°C  is  1.  18  gr/cc.  Vapor  pressure  measurements  were 
attempted  but  these  were  unsuccessful  due  to  the  slow  but  significant  de¬ 
composition  of  the  material.  Nevertheless,  it  could  be  ascertained  from 
the  reaction  intermediate  production  runs,  in  a  low  pressure  apparatus, 
that  the  vapor  pressure  of  the  material  is  less  than  0.  5  mmHg  at  room 
temperature. 
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Various  chemical  analyses  were  performed  on  several  samples  of 
the  material.  From  the  results  of  elemental  analyses,  an  empirical 
formula  for  the  reaction  intermediate  was  computed  to  be  C4H24N90£. 
However  the  fresh  material  appears  to  be  a  mixture  of  at  least  four 
compounds,  the  principal  ones  being  MMH*  H2O  and  MMH-  HNO3.  The 
effect  of  the  gradual  decomposition,  referred  to  above,  is  apparently  to 
reduce  the  percentage  of  MMH-  HNO3  in  the  sample  and  to  convert  the 
monohydrate  of  MMH  to  the  dihydrate,  MMH-  2H2O.  Also  CO2  is  produced 
and  it  is  believed  to  exist  in  the  degraded  samples  as  CH3N2H2COON2H4CH3. 

The  heat  of  formation  of  fresh  reaction  intermediate  was  determined 
from  measurements  of  the  heat  of  combustion  of  the  material  with  oxygen. 
The  heat  of  formation  for  the  liquid  material  was  found  to  be  -0.  509  kcal/gr, 
using  the  empirical  formula  given  above.  Then,  for  an  overall  stoichio¬ 
metry  given  by: 

3N02(g)  +  4CH3N2H3{g)  — *►  C4H24N9C>6(1)  +  Nz(g) 
the  heat  of  this  reaction  for  formation  of  the  reaction  intermediate  during 
an  ignition  delay  period  is  calculated  to  be  AHr=  -66.  2  kcal/mole  of  MMH. 
This  value  is  only  21%  of  the  total  heat  of  complete  combustion  on  a  fuel 
basis. 

To  determine  the  rate  of  formation  of  the  reaction  intermediate,  a 
low  pressure  flow  tube  apparatus  was  used  in  which  the  reactant  vapors 
were  greatly  diluted  with  nitrogen.  The  temperature  profile  along  the  tube 
axis  was  measured  by  fine  wire  Pt-10%  Rh/Pt  thermocouples  positioned  at 
known  distances  from  the  mixing  point  of  the  reactants.  The  overall  order 
of  reaction  for  formation  of  vapor  phase  reaction  intermediate  was  found 
to  be  approximately  two  from  testa  at  elevated  (75°C)  initial  temperature 
and  various  total  pressures  (30  to  60  ramHg).  In  the  room  temperature 
tests,  condensation  of  the  reaction  intermediate  occurred  and  thus 
activation  energy  for  formation  of  reaction  intermediate  could  not  be 
extracted  from  the  effect  of  initial  temperature  on  temperature  rise  along 
the  flow  tube.  The  activation  energy  must  be  low,  however,  because  the 
reaction  proceeds,  and  proceeds  rapidly,  at  room  temperature.  Further, 
the  activation  energy  must  be  less  than  that  of  the  ignition  reactions  since 
the  material  forms  at  conditions  under  which  ignition  does  not  occur  (in 
the  flow  tube,  for  instance). 

Formation  of  the  reaction  intermediate  during  an  ignition  delay 
period  was  incorporated  into  the  ignition  delay  model  by  adding  a  heat 
source  and  a  vapor  mass  sink.  The  heat  source  accounts  for  the  heat 
evolved  as  the  material  forms  and  the  vapor  mass  sink  accounts  for 
consumption  of  vapor  phase  reactants  as  the  condensed  phase  reaction 
intermediate  is  produced. 
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In  the  calculations  for  ignition  delay  in  an  engine,  an  activation 
energy  of  4.  5  kcal/mole  was  used  for  the  formation  reaction  of  the 
reaction  intermediate.  Since  the  activation  energy  for  ignition  of  these 
propellants  is  5.  2  kcal/mole  (determined  in  a  previous  program),  the 
selected  value  of  4.  5  kcal/mole  for  the  activation  energy  for  formation 
of  the  reaction  intermediate  is  reasonable. 


Ignition  was  determined  analytically  by  summing  over  successive 
50  microsecond  time  intervals  of  the  solution  the  fractional  degree  of 
ignition  computed  for  each  time  interval.  The  fractional  degree  of 
ignition  for  each  time  interval  is  the  ratio  of  the  length  of  the  time  interval 
(50  microseconds)  to  the  ignition  delay  time  calculated  for  the  conditions  of 
pressure,  temperature  and  composition  of  the  reactant  vapors  prevailing  at 
that  time.  The  analytical  expression  used  to  calculate  the  ignition  delay 
time  for  each  time  interval  (derived  in  a  previous  program)  is: 
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where  Pe  =  pf  +  pox  and  a  =  pox .  pox,  pf  and  T  are  obtained  for  each 
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time  interval  from  the  pre-ignition  chamber  pressurization  analysis  which 
also  takes  into  account  formation  of  the  reaction  intermediate.  The 
fractional  degrees  of  ignition  are  summed  over  successive  time  intervals 
until  the  fraction  becomes  unity,  indicating  ignition.  The  ignition  delay 
time  in  an  engine  is  that  time  in  the  computer  solution  at  which  the  fraction 
becomes  unity. 


A  series  of  calculations  was  made  in  which  Ajnt,  the  frequency 
factor  for  formation  of  reaction  intermediate,  was  varied  in  order  to 
obtain  for  one  engine  configuration  (design  chamber  pressure  of  200  psia) 
the  best  fit  between  (1)  calculated  and  experimental  pre-ignition  chamber 
pressure-time  histories  and  (2)  calculated  and  experimental  ignition  delay 
times.  A  very  good  fit  was  obtained  for  Ajnt  =  1.  7  x  10  11  cal.  cc/mole2-sec. 
This  value  was  then  used  in  a  computer  run  for  a  second  engine  configura¬ 
tion  (design  chamber  pressure  of  75  psia)  with  the  result  that,  again,  a 
good  fit  was  obtained.  For  each  engine  configuration,  the  calculated  and 
experimental  ignition  delay  times  agreed  to  within  1%. 

The  calculations  indicate  that  the  reaction  intermediate  plays  an 
important  role  in  ignition  of  these  hypergols.  Although  formation  of  the 
reaction  intermediate  removes  mass  from  the  vapor  phase  and  thereby 
tends  to  retard  pre-ignition  chamber  pressurization,  the  heat  which  i« 
evolved  more  than  compensates  for  the  vapor  mass  removal  and  finally 
results  in  markedly  shorter  ignition  delay  times.  A  further  consequence 
of  the  evolved  heat,  however,  is  that  the  reactant  temperatures  are 
substantially  raised  during  the  ignition  delay  period  to  levels  at  which, 
first,  condensation  of  the  reaction  intermediate  will  no  longer  occur, 
second,  the  composition  of  subsequently  formed  reaction  intermediate  will 
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likely  differ  from  the  composition  of  material  formed  earlier  in  the  ignition 
delay  period,  and  third,  the  material  will  undergo  decomposition  or 
possibly  reaction  wiJi  additional  oxidizer.  Thus,  although  considerable 
insight  has  been  obtained  into  the  many  events  occurring  during  an  ignition 
delay  time,  further  study  of  the  latter  part  of  the  ignition  delay  period  is 
needed. 

The  computed  results  predict  that  many  of  the  smaller  N2O4  drops 
cool  evaporatively  to  a  point  at  which  MMH  vapor  will  tend  to  condense 
upon  them.  This  previously  unanticipated  phenomenon  raises  a  question 
as  to7 what  chemical  reaction  will  occur  at.the  oxidizer  liquid/fuel  vapor 
interface.  The  reverse  situation  was  found  not  to  apply,  i.  e.  fuel  drops 
do  not  cool  sufficiently  to  cause  condensation  of  oxidizer  vapor  upon  them. 

From  the  thermochemical  determinations  performed  on  the  reaction 
intermediate,  it  is  found  that  ixie  reaction  between  the  condensed  phase 
material  and  additional  oxidizer  is  very  energetic  and  is  a  likely 
contributor  to  severe  ignition  spikes.  The  calculations  of  ignition  delay 
in  an  engine  indicate  that  sufficient  oxidizer  exists  in  the  vapor  state  at 
ignition  to  react  with  all  of  the  computed  amount  of  reaction  intermediate 
formed  during  the  delay  period.  Therefore,  for  an  increased  understand¬ 
ing  of  severe  ignition  pressure  spikes,  it  is  necessary  to  determine  the 
energy  release  rate  of  the  reaction  between  oxidizer  and  the  reaction 
intermediate. 

The  present  hypergolic  ignition  delay  model,  though  not  complete 
in  every  detail,  is  nevertheless  suitable  for  developing  design  information 
pertaining  to  ignition  phenomena.  For  example,  the  effects  of  propellant 
leads/lags  and  various  transient  flow  increase  schemes  for  each  propellant 
from  time  zero  can  be  analytically  evaluated  as  to  their  effect  on  ignition 
delay  time  and  the  amount  of  reaction  intermediate  formed.  This 
information  can  assist  engine  development  programs  by  economically 
defining  specific  development  objectives. 
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SECTION  III 

PHASE  I  -  PHYSICAL  KINETICS  ASPEC1S  OF  PRE¬ 
IGNITION  CHAMBER  PRESSURIZATION 


The  hypergolic  ignition  delay  model  developed  previously  under 
Contract  AF04(6l  1 ) - 994 6  is  based  on  both  physical  kinetics  of  droplet 
evaporation  and  chemical  kinetics  of  igniiion  reactions.  The  resulting 
model  consists  of  two  cou  tied  sets  of  equations  (1).  In  one  set,  the 
pre-ignition  pressure-time  history  In  the  chamber  is  calculated  taking 
into  account  the  physical  properties  of  the  propellants,  the  rate  of 
propellant  injection,  propellant  evaporation  and  condensation,  heat 
transfer  between  vapor  and  droplets,  droplet  residence  time,  rocket 
engine  chamber  geometry,  and  propellant  vapor  efflux  through  the 
nozzle.  In  the  other  set,  the  chemical  delay  time  as  a  function  of 
pressure,  temperature  and  vaj  ur-phase  composition  is  determined. 

In  this  section  of  the  report,  we  are  concerned  with  the  first-mentioned 
set  which  deals  with  chamber  pressurization. 

It  was  found  from  the  initial  model  development  effort  that  the 
vaporizing  vapor/drop  system  is  non-adiabatic  (i.  e.  ,  significant  heat 
transfer  with  thrust  chamber  walls  occurs),  and  that  a  short  but  finite 
time  is  required,  following  flow  initiation,  for  the  propellants  to  attain 
their  nominal  steady-state  flowrates.  Both  phenomena  affect  pre¬ 
ignition  chamber  pressurization  and  thereby  the  ignition  delay  time  as 
well. 

In  this  section  of  the  report,  the  experimental  and  analytical  work 
performed  to  incorporate  these  phenomena  into  the  mahematical  model 
is  described.  In  the  following  sub-sections,  the  equations  describing 
pre-ignition  chamber  pressurization  due  to  propellant  vaporization  are 
reviewed  (from  Ref.  1);  then  the  experimental  program  which  includes  both 
cold  flow  and  hot  firing  tests  is  discussed.  Finally  in  this  section,  the 
analytical  results  to  include  injector  flashing  and  wall  heating  effects  in 
the  model  are  given. 

A.  THEORY  -  CHAMBER  PRESSURIZATION  DUE  TO  PROPELLANT 

VAPORIZATION 

The  injected  propellant  spray  is  represented  mathematically  by  a 
three-drop-size  spray  model.  The  three  drop  sizes  are  obtained  by 
applying  a  logarithmiconormal  distribution  (with  a  specific  geometric 
standard  deviation  depending  on  the  type  of  injector)  to  a  given  mass 
median  drop  size  and  then  selecting  three  radii  such  that  30%  of  the 
injected  weight  flow  is  in  drops  of  the  smallest  radius,  40%  in  the  inter¬ 
mediate  size  drops,  and  30%  in  the  largest  drops  (4).  Thus  there  are 
three  classes  of  drops,  1  — i  ~3.  The  number  of  drops  in  each  class  is 
obtained  from 
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M  3  W1 

Nl  "  4*  r\  /*i  where  Wt  =  0.  3  Wj  n\ 

W,  =  0.  4  W: 

W,  =  0.  3  Wj 

and  Wj  =  wAtjyj,  the  mass  of  propellant  enjected  during  one  time  interval. 

Pressurization  of  a  thrust  chamber  is  treated  mathematically  as  a 
sequence  of  steady  state  processes  in  very  short  time  intervals.  At  the 
start  of  each  time  interval,  a  new  three-class  ensemble  of  drops  enters  the 
thrust  chamber.  These  undergo  evaporation  during  the  time  interval  as  do 
the  drop  ensembles  which  entered  previously.  At  any  time,  each  class  of 
drops  has  a  unique  radius,  temperature  and  fraction  of  an  individual  drop 
that  is  frozen. 


Each  drop  is  identified  by  i,  j  and  N  numbers  which  indicate 
respectively  the  class  to  which  the  drop  belongs,  the  time  interval  in  which 
the  drop  entered  the  thrust  chamber,  and  the  time  interval  in  question. 

Thus  the  possible  values  of  i,  j  and  N  are:  l^i&3,  ls&js&N,  and 
N  =  1,  2.  .  . 


Based  on  the  kinetic  theory  of  gases,  the  mass  evaporated  from  a 
single  drop  in  one  time  interval  is  given  by  (5); 

TSMgT 


Gij=f  pv(Tdij)N-r  pg(  yN-iMri/N.i  ^ 


^(^dij^N- 1 


(2) 


where  Py^dij^N-l  is  the  vapor  pressure  corresponding  to  the  temperature 
of  the  i**1  size  drop  of  the  j  "  ensemble  at  the  end  of  the  previous  time 
interval. 


The  assumptions  involved  are: 

1.  The  number  of  molecules  leaving  the  liquid  surface  per  unit  time 

{when  the  ambient  pressure  is  below  the  liquid  vapor  pressure) 
is  the  same  as  the  number  which  impinge  on  the  same  surface 
when  the  liquid  is  in  equilibrium  with  its  vapor. 

2.  The  gas  evaporates  at  the  same  temperature  as  the  liquid  drop 
from  which  it  evolves. 

3.  The  molecular  weight  of  the  vapor  and  liquid  is  the  same. 

4.  Ideal  gas  law. 

The  total  mass  evaporated  during  one  time  interval  depends  in  part  on 
the  total  surface  area  of  all  the  drops  in  the  system.  The  total  surface  area 
depends  in  turn  on  the  number  and  radius  of  the  drops  in  each  class. 
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Multiplying  Eq.  2  by  N. ,  the  number  of  drops  in  each  class  as  given  by 
Eq.  1,  and  then  summing  over  the  i's  and  j's,  one  obtains  for  the  total 
mass  evaporated  during  one  time  interval: 


N  3  N  3 

G  =  £  x  G.  ;Ni  =  2  2 
N  j*l  i=l  1J  1  j=l  i=l 


W'Wn-I 


-P  (T  )  1 

g  g  N-iJ 


3CiaW^j  /  Mjc 

/Vrij^N-  1 


(3) 


where  =  0.  3,  C7  =  0.  4  and  C3  =  0.  3  The  density  of  the  liquid,  ,  is 
assumed  constant. 


Condensation  onto  chamber  walls  can  occur  if  the  gas  pressure  exceeds 
the  propellant  vapor  pressure  evaluated  for  the  wall  temperature.  The 
mass  of  vapors  condensing  on  chamber  walls  in  one  time  interval  is  given 
by: 


G 

w 


=  IPg(V 


N-l 


-  P  (T  )]  a  A  At 

V  w  c 


N 


f  Mgc 

>I2"r(Vn-i 


(4) 


but  with  the  constraint  (h6G  ^00. 

w 


The  mass  efflux  through  the  nozzle  during  one  time  interval  for 
constant  k  for  the  gas  is  given  by 


mno.  '  (P*»N-1  A*A,N 


/hMgc 

V  R(TgL 


R(T«Vl-lVk+l 


hii 


(5) 


The  mass  of  vapor  in  the  chamber  at  the  end  of  the  time  interval 
is  obtained  from  the  perfect  gas  law. 


(mc) 


N 


:  VCM  fPg 


(6) 


N 


The  temperature  of  the  gas  in  the  chamber  at  the  end  of  the  N- —  time 
interval  is  obtained  by  taking  a  mass  weighted  average  of  the  temperature, 
which  includes  the  temperature  of  the  gas  generated  due  to  evaporation 
during  the  time  interval  plus  the  temperature  of  the  gas  left  from  the 
previous  time  interval  minus  the  temperature  of  the  gas  which  condenses 
on  the  chamber  walls  and  that  which  passes  through  the  nozzle.  For  a 
single  specie  system  and  a  constant  vapor  heat  capacity, 

N  3  Qn 

2  2  Gjj( Tdij)N- 1  4  (Tg)  j  S N- 1 “(^w)n  “  ( ^nos^wl  '  t 

r  _ _ _ SL  (7, 

^  Gjj  f  ( mc^. j  -  (Gw)n  -  (mno,)N 
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where  Qpj  ie  the  total  heat  transferred  from  the  vapors  to  all  the  drops  in 
the  system  during  one  time  interval.  is  derived  subsequently. 

A  mass  balance  of  the  system  yields  the  gas  pressure  in  the  chamber 
at  the  end  of  the  time  interval. 

(mc^”(rnc^_|  (8) 

The  terms  on  the  left  hand  side  are  given  by  Eqs.  3,  4,  and  5  while  the 
right  hand  side  is  Eq.  6  evaluated  at  the  end  of  the  N  and  N-l  time 
intervals. 


The  gas  pressure  in  the  chamber  at  the  end  of  the  N^time  interval 
is  therefore: 


*gN'TgN 


&V  '  <Gw'n 


-(m  ) 

noz 


"V 


(9) 


where  T^jq  is  evaluated  from  Eq.  7. 


New  drop  radii  and  drop  temperatures  must  be  calculated  at  the  end 
of  every  time  interval  to  account  for  the  effect  of  evaporation  during  the 
time  interval  on  the  size  and  temperature  of  each  drop.  A  new  radius 
for  each  drop  is  obtained  from  a  mass  balance  on  the  drop  together  with 
Eq.  2,  giving: 


(ruV(r>J>N-lD' 


r-  3[pv(Tdij)jj_j  -  Pg^Tg^N-l^  °  AtN 

r^-c — 

*-  A^VpN-l  \ 

Pwk(Tdij)N  J 

'/) 


(10) 


The  new  drop  temperature  is  obtained  from  an  energy  balance  on  the 
drop,  including  the  effect  of  heat  transfer  from  the  vapors  to  each  drop. 
The  energy  reaching  the  surface  of  a  drop  from  the  gas  during  one  time 
interval  is,  according  to  (4): 


qv  *>N  •‘.Wrlj)f).|ZA.Nt(T|)N  i  -(Tdij)N  il 


<U> 


where 


h  = 


7  -  *  0  PS 

e*"1  *  *  *  ^vhAt^rj^ 

- &. - (2+0.  6  Pr1/* Re*/*] 

^rU*N-l 


-1 


Pr  » 


(12) 


(13) 
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Re  = 


2(rij)N-l  vd  M<Pg)N-l 


HgR<Tg)N-l 

G|j  in  Eq.  12  is  given  by  Eq.  2.  Equation  14  omits  drag  effects. 

An  energy  balance  on  the  drop  is  then 

r'»JVc”I  cpijjTd ij^N-1  *  ^Tdij'Nj  4  9Vjj  ^fN 

from  which  the  new  drop  temperature  at  the  end  of  the  time  interval  is 


(Tdij)j^  =  (^dij^N-i 


3[PV(  Tfjij)  jyj.j  -  Pg^g^N-  J  Q  'r‘j^N- J  Xt,AtN  J  ^7 


P\  cpi  <r^k 


dV)'N-l 


wheie 


3hl<Tg»N-l  -<Tdl)lN.ll  'Vn-I  Za‘N 

A  cpi 


The  temperature  within  the  drop  is  assumed  to  be  uniform.  The  heat 
of  vaporization,  *e,  the  heat  capacity  of  the  liquid,  Cpj,  and  the  liquid 
density,  /^j,  are  assumed  constant.  Values  of  these  quantities 
corresponding  to  the  mid-temperature  of  the  range  encountered  are  used. 

A  drop  being  cooled  by  evaporation  can  arrive  at  its  freezing  point. 

Once  Tjjj  equals  Tpp,  then  Tjjjj  remains  constant  as  the  drop  freezes. 

The  fraction  of  the  drop  that  freezes  in  one  time  interval  due  to  evaporation 
and  including  vapor-droplet  heat  transfer  effects  is  given  by 

V«  =  X1JN  I  w/l  XF  +  ^Vjj  AtN  (18) 


The  solid  fraction  of  the  drop,  or  its  "quality"  Xjj,  increases  in  successive 
intervals  until  the  drop  is  frozen  solid.  Xjj  is  given  by 


,  f  3C‘i  .  CP>  I 


with  the  constraint  that  OdXjj<l.  Gjj  and  •  are  given  by  Eqs.  2  and  17 
respectively.  The  summation  over  successive  time  intervals  for  the  drop 
in  question  begins  when  the  temperature  of  the  drop,  T^..,  reaches  Tpp. 
The  summation  ends  when  the  drop  is  completely  frozen.-*  i.  e.  Xjj  =  \. 

Once  the  drop  is  frozen  solid,  further  evaporation  (sublimation)  causes 
the  temperature  of  the  drop  to  decrease.  Equation  16,  modified  to  reflect 
the  solid  state  of  the  droplet,  becomes  applicable  again. 
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The  +ota.l  energy  reaching  the  surface  of  all  drops  in  the  system 
in  one  tim  interval,  i.  e.  C-g  of  Eq.  7  is: 


QN  i! 


N  3 


M  N  3 
iAtNNis  AS 


^irMrj:)*  _ZNiALj{T. 


where  Nj,  the  total  number  of  drops  in  each  class,  is  given  by  Eq.  1 
and  h  and  Z  are  given  by  Eqs.  13  and  12,  respectii,  iy. 

The  above  equations  together  with  the  vapor  pressure  equation 
for  the  propellant  under  consideration  have  been  programmed  for 
an  IBM  360  Model  50  computer.  The  length  of  the  time  intervals 
Atpj,  used  in  calculations  ranges  from  20  to  50  microsec  .nds 
depending  on  the  particular  propellant  and  n  oto  freon  etry  mder 
consideration. 


Typical  gas  pressure  and  temperature  histories  calculated  by 
the  above  equations  are  given  in  Fig.  1  for  N?C  j  and  MMH  in  a 
specific  engine  configuration  (Configurat  i  #1).  The  calculation 
assumes  a  step  function  increase  in  flov  te  at  time  rero  and  an 
adiabatic  vapor/drop  system. 

In  the  computer  program,  tne  acti*.  lifetime  of  each  drop  is 
limited.  The  drop  residence  time  used  i  defined  simply  as  the 
.cme  required  for  a  drop,  moving  at  :  constant  velocity  equal  to 
the  liquid  injection  velocity,  to  trave'  .  rom  the  injector  orifice  to 
a  point  on  the  chamber  wall  and  then  to  the  nozzle  throat.  The  point 
on  Lie  chamber  wall  is  obtained  from  the  direction  of  the  injected 
propellan.  -  .ream  and  the  radius  of  the  chamber.  Computer  runs 
in  which  cop  residence  times  were  varied  by  a  factor  of  two 
showed  mat  this  factor  is  unimportant  tor  an  adiabatic  vapor /drop 
system.  The  reason  for  this  is  that  once  a  drop  has  been  in  the 
system  for  one  residence  time  period,  its  temperature  is  such  that 
its  vapor  pressure  is  close  to  the  gas  pressure  orevailing  at  that 
t  me.  accordance  with  Eq  the  amount  of  evaporation  (or 
condensation)  from  thaf  time  onward  is  small. 

A*so,  in  the  compute:  program,  a  drop  is  made  inactive  before 
one  residence  time  period  lapses  if  it  freezes  solid,  i.  e.  *ij  =  i-  " 
is  mi  tied  that  due  »c  the  low  conduct*.'.*!/  of  the  solid,  an  excessive 
temperature  gradient  exi  s  in  a  solid  drop  which  markedly  reduces 
the  evaporation  rate  at  the  surface  of  the  drop. 

B.  EXPERIMENTAL  PROGRAM  -  THRUST  CHAMBER  TESTS 

1.  Test  Apparatus  and  Instrumentation 
a.  Experimental  Facilities 
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All  thrust  chamber  tests  were  performed  in  a  large  stainless 
steel  vacuum  tank  to  minimize  the  increase  in  ambient  pressure 
due  to  propellant  vaporization  or  combustion.  The  vacuum  tank, 
which  is  shown  in  Fig.  2,  is  7  ft  in  diame*  r  and  25  ft  long,  pro¬ 
viding  a  volume  of  approximately  1000  ft^.  The  tank  has  four  8  in. 
diameter  ports  along  each  side  for  instrumentation  and  observation 
purposes. 

The  pumping  system  consists  of  a  Kinney  KD  780  pump  having 
a  625  cfm  pumping  capacity  at  a  pressure  of  1  mmHg  and  a  Roots 
6000  blower  which  has  a  pumping  capacity  of  3250  cfm  at  a  pressure 
of  10'1  mmHg.  All  thrust  chamber  tests  were  conducted  at  an 
ambient  pressure  of  0.  3  mmHg  and  were  less  than  120  msec  dura¬ 
tion,  so  the  ambient  pressure  was  essentially  constant  for  the 
entire  test  period.  Schlieren  movies  of  the  flow  issuing  from  an 
injector  were  taken  at  a  higher  pressure,  5  mmHg,  to  permit 
better  definition  of  the  initial  gas  propellant  flow  which  is  difficult 
to  detect  at  the  lower  ambient  pressures. 

The  oxidizer  and  fuel  propellant  systems  consisted  of  a  300  cc 
stainless  steel  tank,  Jamesbury  safety  valve,  Fox  propellant 
solenoid  valve,  and  associated  hand  valves  and  tubing  as  shown  in 
Fig,  3.  Each  set  was  mounted  on  a  separate  plate  which  contained 
provisions  for  all  necessary  external  connections  such  as  pressuriza¬ 
tion  and  vent  systems.  The  propellant  %'alves  and  chamber  were 
mounted  under  the  plates  with  the  injector  facing  downward  as  shown 
in  Fig.  4.  Also  visible  in  this  photograph  are  the  actuators  for  the 
Jamesbury  safety  valves  which  are  mounted  out  of  view  on  the  top 
of  the  plates. 

b.  Experimental  Hardware 

The  50  lb  thrust,  attitude  control  thrust  chamber  consisted  of 
an  injector-valve  assembly,  a  chamber  flange,  a  transparent 
chamber  section,  and  a  stainless  steel  nozzle.  An  exploded  view  of 
a  thrust  chamber  assembly  is  shown  in  Fig.  5,  although  the  particu¬ 
lar  chamber  length  shown  in  the  photograph  was  not  tested  in  the 
present  program.  The  Kistler  water-cooled  adapters  shown  in  this 
picture  were  replaced  with  the  uncooled  adapters  shown  in  Fig.  4. 

The  diaphragms  of  the  flush-mounted  transducers  were  protected 
with  RTV  silicone  rubber.  For  short  run  times  (several  hundred 
milliseconds),  the  coaling  eliminates  thermal  effects  on  the  trans¬ 
ducers  from  either  cold  propellants  or  hot  combustion  gases  which 
would  otherwise  lead  to  erroneous  pressure  signals. 
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Five  chamber  and  nozzle  configurations  were  chosen  from  the 
27  possible  combinations  available  from  the  previous  programs 
under  Contract  AF04(6ll)-9946.  The  design  parameters  for  the 
five  selected  configurations  are  detailed  in  Table  I.  The  trans¬ 
parent  chamber  sections  were  fabricated  from  commercially 
available  sizes  of  acrylic  tubing  cut  to  the  desired  length. 

TABLE  I 

THRUST  CHAMBER  CONFIGURATIONS 


Configuration  No. 

1 

2 

3 

4 

5 

Design  P  ,  psia 

200 

75 

20 

200 

200 

Chamber  I.  D.  ,  in. 

1. 25 

1.  875 

3.  50 

0.  75 

1.  25 

Contraction  Ratio 

9.  00 

7.  90 

7.  68 

3.  22 

9.00 

L*,  in. 

10.  3 

10.  8 

13.0 

10.0 

50.0 

Chamber  Length,  in.  ^ 

0.  25 

0.  25 

0.  25 

2.  21 

4.  69 

Nozzle  Throat  Area,  sq.  in. 

0.  137 

0.  351 

1.  253 

0.  137 

0.  137 

(a)  Length  of  transparent  section  only. 


Each  stainless  steel  nozzle  section  includes  a  3/4  inch  long 
cylindrical  portion  whose  internal  diameter  is  the  same  as  the  inside 
of  the  respective  transparent  chamber  section.  Kistler  pressure 
transducers  were  flush  mounted  in  the  cylindrical  portion  of  the  nozzle 
sections.  The  volume  of  this  section,  as  well  as  that  of  the  convergent 
portion  of  the  nozzle,  was  included  in  the  L*  calculations  for  sizing 
the  transparent  sections. 

The  entire  chamber  assembly  was  held  together  with  1/4  inch 
threaded  rods  which  were  insulated  from  the  metal  nozzle  to  keep 
the  Kistler  transducers  ungrounded.  Buna-N  O-rings  sealed  the 
flanges  and  the  transparent  chamber  sections. 

The  propellant  solenoid  valves  are  1/4  in.  venturi  valves  manu¬ 
factured  by  the  Fox  Valve  Development  Co.  ,  Inc.  These  stainless 
steel  valves  have  a  pintle  with  a  Teflon  poppet  which  provides  a  seal 
in  the  entrance  of  the  venturi  section.  The  propellant  volume  in  the 
valve  downstream  of  the  seat  is  quite  small.  The  valves  were  mounted 
directly  to  the  back  face  of  the  injector  to  minimize  dribble  volumes. 

A  high  voltage  pulse  circuit  was  used  to  reduce  the  valve  opening 
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time  to  about  1  msec. 

Two  injectors  were  used  during  the  program:  a  one-on-one 
doublet  injector  and  a  self-impinging  doublet  injector.  Originally, 
the  one-on-one  doublet  injector  was  to  be  used  for  hot  firings 
whereas  the  self-impinging  injector  was  to  be  used  for  obtain¬ 
ing  chamber  pressurization  data  for  vaporization  of  singly 
flowed  propellants.  This  data  would  be  the  basis  against  which 
computed  chamber  pressurization  histories  (without  chemical 
reaction)  would  be  compared.  The  reason  for  using  a  self- 
impinging  injector  for  singly  flowed  propellants  was  to  duplicate 
the  effects  of  stream  interaction  on  break-up  and  atomization 
of  the  injected  propellants  in  hot  firings. 

Section  views  of  the  two  injectors  are  shown  in  Fig.  6, 
with  the  propellant  valves  outlined.  The  throats  of  the  venturi 
valves  are  approximately  in  line  with  the  rear  face  of  the  in¬ 
jector.  The  single -element  doublet  injector  consists  of  two 
#53  (0.  0595  in.  dia. )  drill  holes  impinging  at  a  90°  angle  0.  14 
inches  from  the  injector  face.  The  self-impinging  injector  was 
designed  to  duplicate  one  propellant  side  of  the  single  -element 
doublet  injector,  but  provide  stream  interaction  effects  on  spray 
break-up.  Therefore  the  impingement  length,  impingement 
angle,  injector  dribble  volume  (0.  0021  cu.  in.  downstream  of 
venturi  hole),  and  liquid  propellant  injection  velocity  (69.  0  and 
71.  5  ft/sec  for  N2O4  and  MMH,  respectively)  are  identical  in 
both  injectors.  The  internal  surface  area  is  unavoidably  larger 
in  the  self-impinging  injector  and  the  propellant  path  is  not  as 
direct  as  would  be  desirable.  The  effect  of  these  differences  is 
discussed  in  Section  III.  B.  2. 

c.  Instrumentation 

In  addition  to  conventional  instrumentation  required  for 
propellant  pressurization  and  system  monitoring,  special 
instrumentation  was  employed  to  determine  ignition  delays  and 
monitor  pressure  transients  during  the  ignition  delay  period. 

This  instrumentation  included  high  speed  schlieren  movies  of 
the  propellant  impingement  zone,  photomultiplier  tube  to  detect 
ignition,  and  high  response  pressure  transducers  and  recording 
equipment. 

The  purpose  of  the  schlieren  system  was  to  investigate 
propellant  stream  characteristics  prior  to  ignition.  Although 
the  schlieren  instrumentation  could  not  be  used  to  detect 
propellant  entry  during  the  actual  thrust  chamber  tests  because 
of  the  poor  optical  properties  of  the  plastic  chambers,  it  was 
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used  to  study  the  injection  characteristics  of  both  injectors  without 
a  chamber  installed.  Times  could  be  determined  from  valve  signal 
to  vapor  entry*  liquid  entry,  and  steady-state  liquid  flow. 

A  two-mirror,  parallel-beam  schlieren  system  was  used 
with  a  Fastax  high-speed  camera  to  obtain  suitable  time  resolution. 

The  essential  components  of  the  system  included  a  Unertl  Model 
BH  6  Ncrmal  and  Color  Schlieren  Source,  a  pair  of  eight-inch 
front- surface  parabolic  mirrors  of  64-inch  fecal  length,  knife-edge, 
and  the  Fastax  camera  capable  of  approximately  16,  000  pictures 
per  second  using  split-frame  optics.  Time  resolution  when  the 
valve  opened  was  approximately  7  frames  (14  pictures)  per  msec 
as  determined  by  a  1000  Hz  timing  light.  The  Fastax  camera  and 
knife  edge  are  visible  in  Fig.  2. 

A  flame  detector  consisting  of  an  RCA  1P28  photomulti¬ 
plier  tube  was  us>.  i  to  determine  the  time  at  which  ignition  occurred. 
The  tube  is  sensitive  to  wavelengths  from  about  2200A  to  600oX. 
but  no  filter  was  required  since  the  tank  was  not  illuminated  for  the 
thrust  chamber  tests.  It  could  sense  ignition  through  the  plastic 
chambers.  A  short -duration,  single-flash  strobe  light  was  also 
detected  by  the  photomultiplier  tube,  as  well  as  by  the  camera,  for 
a  time  reference  to  correlate  the  various  instrumentation. 

Chamber  pressure  was  measured  by  Model  701A  and  603A 
Kistler  pressure  transducers  and  Model  566  Kistler  charge  amplifiers. 
The  two  transducers  were  mounted  in  the  same  plane  at  the  nozzle 
end  of  the  thrust  chamber  as  shown  in  Fig.  4.  The  7G1A  transducer 
with  a  rise  time  of  6  microseconds  and  high  sensitivity  was  used  to 
detect  chamber  pressure  rise  due  to  propellant  entry  and  vaporization 
prior  to  ignition.  The  603A  transducer,  with  less  sensitivity  but  with 
a  resonant  frequency  over  400,000  Hz  and  a  rise  time  of  1  micro¬ 
second,  was  used  for  the  pressure  transients  Lt  and  after  ignition. 

Two  additional  603A  transducers  in  the  propellant  lines  upstream  of 
the  Fox  propellant  valves  (Fig.  4)  measured  the  feed  line  pressure 
transients  when  the  valves  opened. 
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Propellant  temperatures  were  monitored  at  the  propellant 
valves  using  metal  sheathed  copper -constantan  thermocouples. 

A  bare  bead  thermocouple  on  the  back  face  of  the  injector  was 
used  to  measure  initial  hardware  temperature.  All  temperatures 
were  recorded  prior  to  a  run  on  a  single  Brown  strip  chart 
recorder  using  a  selector  switch. 

The  transient  data  were  recorded  on  a  Tektronix  Type  555 
Dual-Beam  Oscilloscope  with  type  CA  amplifiers  to  permit 
simultaneous  display  of  four  parameters.  The  four  parameters 
selected  differed  for  the  cold  flow  tests  and  the  ignition  tests. 

For  cold  flow  of  one  propellant,  the  parameters  were  the  701A 
chamber  pressure  transducer,  the  603A  valve  inlet  pressure 
transducer  for  the  respective  propellant,  valve  current,  and 
also  the  photomultiplier  tube  signal  when  the  schlieren  movies 
were  taken.  For  ignition  tests,  the  701A  chamber  pressure 
transducer  was  used  to  measure  pressure  rise  prior  to  ignition 
due  to  both  propellants  vaporizing  and  reacting,  the  603A  chamber 
pressure  transducer  was  used  for  ignition  spikes  and  transients, 
the  photomultiplier  tube  for  ignition  detection,  and  the  combined 
valve  currents  for  valve  operation.  Oscilloscope  sweep  rates 
varied  from  1  to  20  msec /cm  depending  on  the  test  information 
desired. 


2.  Experimental  Results 

Experimental  chamber  pressure-time  curves  for  vaporiza¬ 
tion  of  singly -flowed  propellants  in  thrust  chambers  of  various  con¬ 
figurations  were  measured  to  serve  as  the  basis  against  which  to 
compare  the  computed  vaporization  pressurization  curves.  Two 
injectors  were  used  in  the  tests:  a  self-impinging  doublet  injector 
and  a  single  element  one-on-one  doublet  injector. 

a.  Self -impinging  Injector  vs.  1 -on- 1  Doublet  Iniector. 

The  self-impinging  injector  was  used  in  order  to 
obtain  the  effects  of  stream  interaction  on  spray  break-up  for 
a  singly-flowed  propellant.  Although  the  self-impinging  injector 
maintains  the  same  manifold  volume,  impingement  angle,  liquid 
injection  velocity  and  impingement  length  as  the  single  element 
doublet  injector,  vaporization  chamber  pressure  rise  times 
were  markedly  different  for  the  two  injectors,  as  shown  in  Fig. 

7.  The  figure  gives  typical  oscillograms  of  tocts  v.'ith  cnch 
propellant  singly  flowed  through  each  injector  in  two  different 
engine  conf.gurations.  In  each  oscillogram,  the  middle  curve  at 
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Self-Impinging 
Injector _ 

Run  No.  3 


Propellant 


1  -on-  1  Injector 
Run  No.  53 


Figure  7.  Oscillograms  oi  inn 
Inj.  cturs  in  Each  oi 


Flowed  N7O4  and  N'.MH  Through  Two 
o  Engine  Configurations 
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the  left  side  (start  of  trace)  is  the  chamber  pressure  as  measured 
by  a  Kistler  701A  transducer.  The  upper  curve  gives  propellant 
valve  iulet  pressure  (Kistler  603A  transducers)  and  the  lower 
curve  is  propellant  valve  current.  Deflection  of  the  latter  is  down¬ 
ward  in  the  case  of  N2O4  and  upward  in  the  case  of  MMH.  Pertin¬ 
ent  data  and  scale  factors  for  the  eight  oscillograms  are  given  in 
Table  II. 

With  both  engines  and  with  both  propellants,  chamber  pressure 
rise  times  are  seen  to  be  significantly  slower  with  the  self-imping¬ 
ing  injector.  The  four  chamber  pressurization  curves  for  Engine 
Configuration  No*-*?  are  plotted  on  a  common  scale  in  Fig.  8.  The 
longer  rise  times  of  the  self-impinging  injector  are  pronounced 

The  transient  flow  characteristics  of  each  injector  with  each 
propellant  singly  flowed  were  determined  from  high  speed  schlieren 
movies  (approx.  14,000  pictures  per  second).  Because  of  the  poor 
optical  quality  of  the  acrylic  thrust  chamber  sections,  the  schlieren 
movies  were  taken  without  chambers  or  nozzles  attached  to  the 
injectors.  A  summary  of  the  results  is  given  in  Table  III  which 
gives  the  average  times  from  valve  signal  to  emergence  of  vapor, 
vapor/liquid  mixiure,  all  liquid,  and  steady  state  liquid  flow.  The 
transition  from  vapor  flow  to  vapor/liquid  mixture  flow  is  sudden 
and  is  "caught"  in  a  single  picture  of  the  movie.  By  contrast,  the 
transitions  to  all  liquid  flow  and  to  steady  state  liquid  flow  are 
gradual.  Consequently  the  times  from  valve  signal  to  these  latter 
transitions  are  subject  to  an  error  of  perhaps  one  millesecond. 
Nevertheless,  it  is  clear  that  steady  state  liquid  flow  is  achieved 
more  quiekly  with  the  one-on-one  doublet  injector.  Although  the 
manifold  volumes  are  the  same  for  each  propellant  in  the  two  in¬ 
jectors,  the  more  tortuous  path  in  the  self -impinging  injecto- 
(Fig.  6)  and  its  higher  internal  surface -to-volume  ratio  (which  tends 
to  inert  se  propellant  flashing)  retard  attainment  of  steady  state 
liquid  flow.  The  reduced  mass  throughput  do-ing  the  flow  build-up 
period  with  the  self -impinging  injector  causes  a  slower  pressure 
rise  in  a  thrust  chamber  and  would  lead  to  a  longer  ignition  delay 
time.  Thus,  the  self-impingin^  injector  flow  did  not  properly 
simulate  the  hot  firing  flow  of  one  propellant  through  the  one-  n-one 
injector.  The  latter  injector,  then,  was  used  subsequently  for  both 
cold  flow  tests  of  singly  flow  \  .opellants  and  for  hot  firings  as 
well. 
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TABLE IU 

SUMMARY  OF  SCHLIEREN  CHARACTERIZATION  OF  INJECTORS 


Time  from  Valve  signal 

to:* 

Propellant  Injector 

Vapor 

msec 

Vapor/ 

Liquid 

msec 

All 

Liquid 

msec 

St,  St. 
Liquid** 
msec 

Run  Nos. 

n2°4 

1 -on-i 

- 

4.  6 

6.  1 

7.  3 

44,  44a, 45 

self 

1.  8 

4.  2 

7.  0 

H.  5 

11-14,  16 

impinging 

MMH 

1 -on- 1 

3.  6 

3.  9 

- 

7.  6 

46,47 

?elf 

2.  0 

4.  5 

6.  7 

12.  0 

17-19 

impinging 

*Valve  cpenirg  ti meo* 0.  8  msec  for  each  valve  based  on  valve  current 
traces. 

**St.  St.  =  Steady  State. 

o.  Cold  Flo\  d  Kot  Firing  Pre  Ignition  Chamber  Pressurization 

Experimental  charnher  p  ssure -time  curves  were  measured  for 
each  propellant  singly  flowed  using  the  one-on-one  doublet  injector 
and  all  five  engine  configurations  (Table  I).  Similar  curves  for  hot 
firings  in  the  five  engines  were  also  obtained.  Figures  9  and  10  give 
typical  oscillograms  for  each  propellant  singly  flowed  and  for  a.  hot 
firing  in  each  of  the  five  engine  configurations.  Corresponding 
oscillograms  in  each  figure  differ  only  in  the  time  base,  Fig-  9  having 
the  faster  sweep  rates. 

For  all  cold  flow  tests  in  the  two  figures,  the  upper  curve  at  the 
left  side  {start  of  trace)  of  each  oscillogram  is  valve  inlet  pressure 
(Kistler  603A  transducer).  The  middle  curve  gives  the  chamber  pressure 
rise(Kistlcr  701A  transducer)  due  to  vaporization  of  the  propellant.  The 
lower  curv'-  is  valve  current. 

For  the  hot  firings,  the  two  upper  curves  are  chamber  pressure. 

Thu  first  of  these  two  to  rise  is  the  signal  from  the  701A  transducer. 

When  ignition  occurs,  this  trace  goes  off  scale  and  is  not  seen  there¬ 
after.  The  second  of  the  two  upper  curves  is  the  signal  from  a  603A 
trr  nsducer  also  flush  mounted  in  the  thrust  chamber.  The  gam  for 
this  transducer  is  reduced  so  that  chamber  pressure  during  combus¬ 
tion  is  recorded.  The  lower  two  traces  in  the  hot  firings  are  photo- 
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rfe°4  Cold  Flow 


MMH  Cold  Flow 


n2o4/mmh 

Hot  Firing 


Engine 

Config. 


#1 


Run  No.  53  Run  No.  55 


Run  No.  37 


73  .  75  77 


Time  Base  =  5  msec/cm  for  Runs  73  and  75 
=  2  mspc/cm  for  All  Other  Runs 


Figure  9.  Oscillograms  (Faster  Sweep)  of  Cold  Flow  and  Hot  Firings  in 

Five  Engine  Configurations 
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N?Oa  Cold  Flow 


Engine  Run  No.  52 
Conti  g. 


MMH  Cold  Flow 
Run  No.  54 


n2o4/mmh 

Hot  Firing 
Run  No.  58 
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multiplier  signal  and  a  combined  valve  current  signal.  The  latter  is 
the  one  that  remains  on- scale  throughout  the  runs.  When  the  photo¬ 
multiplier  receives  light  from  ignition  or  combustion,  the  signal 
deflects  off-scale  downward.  Pertinent  data  and  scale  factors  for 
the  oscillograms  of  the  figures  are  given  in  Table  IV. 

Ignition  delays  and  pressure  in  the  chambers  at  ignition  were 
obtained  for  the  five  engine  configurations  from  these  hot  firing 
oscillograms  and  those  of  repeated  runs.  The  results  are  summarized 
in  Table  V.  As  noted  in  the  table,  the  ignition  delays  are  times  from 
valve  signal  to  ignition.  The  valve  current  traces  in  the  cold  flow  runs 
indicate  a  valve  opening  time  of  less  than  one  millisecond.  Because 
the  pressure  in  the  chamber  just  prior  to  ignition  is  rising  rapidly 
generally,  the  pressure  at  ignition  is  difficult  to  specify  with  good 
precision  (for  example,  see  Run  No.  61  in  Fig.  9).  In  a  few  cases, 
the  pressures  at  ignition  as  determined  from  repeated  runs  varied  by 
as  much,  as  20%. 

The  two  cold  flow  chamber  pressurization  curves  and  the  hot 
firing  pre-ignition  chamber  pressurization  curve  for  Engine  Configura¬ 
tions  1  and  2  are  plotted  on  common  scales  in  Figs.  11  and  12  res¬ 
pectively.  In  each  figure,  the  dashed  curve  is  simply  the  sum  of  the 
cold  flow  chamber  pressurization  curves  for  each  reactant.  A  com¬ 
parison  of  these  sum-of-partial-pressures  curves  ana  their  respec¬ 
tive  hot  firing  pre -ignition  curves  shows  good  agreement  for  the  first 
millisecond  of  the  rising  pressures.  A  gradual  divergence  then  begins, 
the  hot  firing  pre-ignition  curves  rising  more  steeply  than  the  sum-of- 
partial-pressures  curves.  For  both  engine  configurations,  the  slopes 
of  the  two  hot-firing  pre-ignition  curves  increase  with  time  until  ignition 
occurs,  which  causes  a  discontinuity  in  the  curves  (Figs.  9  and  10). 

The  divergence,  which  reaches  100-1  30  mmHg  at  ignition,  indicates 
appreciable  pre-ignition  reactions.  The  divergence  and  its  cause  are 
the  subjects  of  Section  IV  of  this  report.  In  the  present  section,  we 
are  concerned  with  analytically  predicting  the  cold  flow  chamber 
pressurization  curves  of  each  propellant  singly  flowed. 
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TABLE  IV 

PERTINENT  DATA  AND  SCALE  FACTORS  FOR 
OSCILLOGRAMS  OF  FIGURES  9  AND  ’0 


Propellant  and  _ Scale  Factors 


Run  No. 

Propellant 

Hardware  Temp. 

— °FTT° — 

Pr-7Q1A 

mmHg/cm 

Pr-603A 
psi/  cm 

Sweep  Rate 
msec/ cm 

37 

NTO/MMH 

77 

103 

75.5 

2 

38 

NTO/MMH 

77 

103 

75.5 

5 

41 

NTO/MMH 

77 

103 

30.2 

5 

42 

NTO/MMH 

77 

103 

30.2 

2 

48 

NTO 

78 

51.5 

- 

10 

49 

MMH 

78 

51.5 

2 

50 

MMH 

78 

20.6 

- 

10 

51 

MMH 

78 

20.6 

- 

2 

52 

NTO 

79 

103 

- 

10 

53 

NTO 

79 

103 

- 

2 

54 

MMH 

80 

20.6 

_ 

10 

55 

MMH 

80 

20.6 

- 

2 

56 

NTO 

80 

103 

- 

10 

57 

NTO 

80 

103 

- 

2 

58 

MMH 

81 

20.6 

- 

10 

59 

MMH 

81 

20.6 

2 

61 

NTO/MMH 

81 

103 

75.5 

2 

64 

NTO/MMH 

81 

206 

151 

5 

65 

NTO 

82 

20.6 

• 

10 

66 

NTO 

82 

20.6 

- 

2 

67 

MMH 

83 

5.2 

_ 

10 

68 

MMH 

83 

5.2 

- 

2 

69 

NTO/MMH 

83 

51.5 

7.6 

5 

70 

NTO/MMH 

83 

51.  5 

7.6 

2 

72 

NTO 

72 

103 

- 

10 

73 

NTO 

73 

51.  5 

5 

74 

MMH 

74 

20.6 

• 

10 

75 

MMH 

74 

10.  3 

- 

5 

76 

NTO/MMH 

74 

103 

75.  5 

5 

77 

NTO/MMH 

74 

103 

75.5 

2 
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TABLE  IV  -  cont'd 

PERTINENT  DATA  AND  SCALE  FACTORS  FOR 
OSCILLOGRAMS  OF  FIGURES  9  AND  10 


Parameters  Common  to  All  Runs  of  Figures  9  and  10  : 


Initial  Ambient  Pressure 
Propellant  Tank  Pressure  -  N2O4 
Propellant  Tank  Pressure  -  MMH 
Steady  State  N2O4  Flowrate 
Steady  State  MMH  Flowrate 
Injector:  1-on-l  doublet 
Scale  Factors: 

Valve  Inlet  Pressure  -  Ox 

Valve  Inlet  Pressure  -  Fuel 


=  0.  3  mmHg 
=  740  psig 

=  710  psig 
=  0.120  lbm/sec 
=  0.0751bm/sec 


83  psi/cm 
7  1  psi/  cm 
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TABLE  V 

SUMMARY  OF  EXPERIMENTAL  IGNITION  DELAYS 
AND  CHAMBER  PRESSURES  AT  IGNITION 


APPROX.  AVERAGE 

ENGINE  AVERAGE  CHAMBER  PRESSURE 


CONFIGURATION  NO. 

IGNITION  DELAY* 

AT  IGNITIOr 

1 

msec 

5.7  -  0.2 

mmHg 

260 

2 

6.6  +  0.  2 

230 

3 

9.4  +  0.4 

70 

4 

6.0  t  0.  2 

330 

5 

6.7  1  0.  2 

130 

*  Valve  signal  to  ignition, 


Chamber  Pressure, 
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0  2  4  6  8  10  12  14 


Time,  msec 


Figure  12.  Experimental  Chamber  Pressurization  Curves 
for  Cold  Flow  and  Hot  Firing  -  Engine 
Configuration  No.  2. 
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c.  Cold  Flow  and  Hot  Firing  Chamber  Wall  Temperature  Histories 

Experimental  thrust  chamber  wall  temperature  histories  were 
measured  for  singly  flowed  N2O4  and  MMH  in  Engine  Configurations 
l  and  2.  Wall  temperature  histories  were  also  obtained  for  hot 
firings  in  Engine  Configuration  No.  1. 

The  measurements  were  made  using  a  high- response  (approxi¬ 
mately  10  microsecond  response  time)  chromel-alumel  Nanmac 
surface  junction  thermocouple  flush  mounted  in  the  cylindrical 
portion  of  the  nozzle  section  and  diametrically  opposite  the  Kistler 
701A  transducer.  The  cold  junction  of  the  thermocouple  was  at 
32°F.  Oscillograms  of  cold  flow  tests  with  each  propellant  in  the 
two  engines  are  given  in  Fig.  13. 

In  the  oscillograms,  the  upper  trace  is  valve  current,  the 
middle  trace  is  chamber  pressure  (Kistler  701A),  bottom  trace  is 
the  inside  surface  wall  temperature  as  given  by  the  flush  mounted 
Nanmac  tnermocouple.  The  polarity  of  the  recorder  is  such  that  an 
upward  deflection  of  the  thermocouple  trace  means  a  decrease  in 
temperature.  Time  bases  for  the  oscillograms  are  2  msec/cm  for 
the  left-hand  column  and  20  msec /err.  for  the  right  hand  column. 

Other  pertinent  data  and  scale  factors  are  given  in  Table  VI. 

Although  the  same  tank  pressures  were  used  as  in  the  earlier 
runs  given  in  Figs.  9  and  10,  the  flow-rate  of  the  oxidir-r  is  lower  in 
the  present  runs  as  indicated  by  the  chamber  pressurization  curves 
for  in  Fig  13  versus  the  c  or  responding  curves  of  Figs.  9  and 

’  0.  Apparently  partial  plugging  c'  the  oxidizer  passages  occurred 
during  the  interval  between  the  two  groups  of  tests.  The  reduced 
NiC>4  flowrate  however  does  not  negate  the  wall  temperature  trends 
measured  in  these  runs. 

In  the  case  of  a  singly  t'lowea  propellant  through  the  one-on-one 
doublet  injector,  the  propellant  initially  impinges  on  the  chamber  wall 
at  a  localized  area.  Chamber  wall  cooling  is  therefore  non-uniform 
as  may  be  seen  by  comparing  the  top  two  oscillograms  of  Fig.  1  3 
with  the  two  directly  below.  The  nozzle  section,  and  therefore  the 
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Engine 

Config.  Propellant  Run  No.  83 


Run  No.  85 


Nanmac 
T/C  Location 


Opposite 

Propellant 

Impingement 

Side 


89 


On  Propellant 

Impingement. 

Side 


2 


MMH 


92  91 


Scope  Sweep  Rate 

2  msec/cm  20  msec/cm 


On  Propellant 

Impingement 

Side 


Figure  13.  Experimental  Chamber  Wall  Temperature  Histories 
for  Singly  Flowed  N2O4  and  MMH  in  Engine  Configurations  1  and  2 


PERTINENT  DATA  AND  SCALE  FACTORS  FOR  OSCILLOGRAMS  OF  FIGURE  13 
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thermocouple,  could  be  rotated  relative  to  the  injector  in  90° 
steps.  In  Runs  83  and  85,  the  thermocouple  was  located  on  the 
side  of  the  chamber  opposite  the  localized  impingement  area  and 
it  registered  only  a  slow,  gradual  decrease  in  wall  temperature. 

In  Runs  87  and  88  as  well  as  in  the  remaining  runs  of  Fig.  13, 
the  thermocouple  was  located  on  the  same  side  of  the  chamber  as 
the  principal  localized  impingement  area.  In  these  cases,  wall 
temperature  begins  to  decrease  very  soon  after  the  pressure 
begins  to  rise  in  the  chamber.  The  extent  of  wail  cooling  near 
the  localized  impingement  area  is  markedly  greater  with 
as  propellant  than  with  MMH. 

The  oscillograms  with  the  slowe  .  sweep  rates  indicate  that 
the  temperature  of  the  wall  near  the  localized  impingement  area 
decreases  fairly  quickly  to  a  level  which  remains  approximately 
constant  until  shutdown.  With  a  given  propellant  the  nearly 
constant  wall  temperature  during  a  run  is  slightly  lower  for  a 
lower  design  chamber  pressure,  as  may  be  seen  by  comparing 
Runs  89  with  88  for  and  91  with  94  for  MMH. 

At  shutdown,  wall  temperature  undergoes  further  cooling,  as 
shown  most  noticeably  in  Run  88.  The  increased  wall  cooling 
during  tail -off  is  attributed  to  vaporization  of  propellant  which  had 
coated  the  wall  surface  during  the  run  time.  This  behavior  is 
observed  to  a  greater  or  lesser  extent  for  both  propellants  and  both 
engine  configurations  tested. 

The  wall  temperature  histories  of  Fig.  13  indicate  significant 
cooling  even  during  typical  ignition  delay  times.  The  start  of  wall 
cooling  lags  only  slightly  (approx.  1/2  msec)  the  start  of  chamber 
pressure  increase  due  to  propellant  vaporization.  However,  cooling 
of  the  walls  is  non-uniform  over  the  inside  surface  of  the  chambers. 
The  curves  suggest  also  that  the  walls  become  coated  with  propellant 
during  the  runs,  the  propellant  acting  as  insulation  to  slow,  together 
with  the  cooled  walls,  subsequent  heat  addition  to  the  vaporizing 
vapor/drop  system. 

Wall  temperature  histories  were  also  measured  during  hot 
firings  of  Engine  Configuration  No.  1.  Oscillograms  for  these 
runs  are  shown  in  Fig.  14  together  with  oscillograms  of  appropriate 
cold  flow  tests  from  Fig.  13.  In  the  oscillograms  for  the  hot  firings, 
the  upper  trace  is  a  combined  valve  current  trace,  the  middle  trace 
is  wall  temperature  (Nanmac  thermocouple)  and  the  bottom  trace  is 
chambe-  pressure  (Kistler  701A).  As  before,  a  decrease  in  tempera¬ 
ture  causes  and  upward  deflection  of  the  thermocouple  trace.  Also, 
the  N2O4  flowrate  is  at  a  reduced  level  as  in  the  cold  flow  runs  of  Fig. 
13.  Pertinent  data  and  scale  factors  for  the  oscillograms  of  Fig.  14 
are  given  in  Table  VII. 
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Nanmsc  Thermocouple  Location  Relative  to  Injector  Orifices: 


Ox  Inlet  Side 
Ox  Fuel 


90°  from  Inlets 
Ox  Fuel 


Fuel  Inlet  Side 
Ox  .  Fu« 


701A 


Propellant 


Run  No.  83 


Run  No.  87 


Time  Base  =  1  msec/cm  for  Runs  97  and  98 
=  2  msec/cm  for  All  Other  Runs 

Figure  14.  Experimental  Cbamoer  Wall  Temperature  Histories  for  Cold 
Flow  and  Hoi  Firings  in  Engine  Configuration  No.  1- 


PERTINENT  DATA  AND  SCALE  FACTORS  FOR  OSCILLOGRAMS  OF  FIGURE  14 
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As  discussed  above,  the  cold  flow  oscillograms  of  Fit?.  14 
indicate  significant  wall  cooling  near  the  localized  impingement 
point  for  each  propellant  singly  flowed.  The  wall  temperature 
traces  for  the  hot  firings  of  Fig.  14  exhibit  a  quite  different 
behavior.  Regardless  of  the  orientation  of  the  thermocouple 
relative  to  the  propellant  injection  orifices,  the  wall  temperature 
trace,  during  the  ignition  delay  period,  initially  deflects  upward 
very  slightly  (slight  wall  temperature  decrease)  and  then  deflects 
noticeably,  though  gradually,  downward  (wall  temperature  increase). 

At  ignition,  a  discontinuity  in  the  downward  deflection  occurs, 
indicating  greater  wall  heating.  The  point  of  interest  here  is  that, 
wall  temperatures  decrease  in  cold  flow  tests  but,  in  hot  firing, 
wall  temperatures  decrease  very  slightly  initially  and  then  in¬ 
crease  during  the  remainder  of  the  ignition  delay  period.  The  wall 
temperature  increase  (obtained  at  all  three  thermocouple  locations) 
during  the  ignition  delay  period  is  clear  evidence  of  heat  evolution 
from  pre-ignition  reactions. 

C.  INJECTOR  FLASHING  EFFECTS  ON  PRE -IGNITION  CHAMBER 

PRESSURIZATION 

High  speed  schlieren  movies  readily  record  that  when  propellant 
emerges  from  an  injector  following  valve  opening,  the  propellant 
emerges  as  a  vapor  initially,  then  as  a  vapor /liquid  mixture,  and 
finally  as  all  liquid  flow.  The  reduced  mass  flow  during  this  flow  build¬ 
up  period  slows  the  vaporization  pressure  rise  in  a  thrust  chamber 
from  the  pressurization  rates  that  would  result  from  a  step-function 
increase  in  propellant  flowrate.  Injector  flashing  effects  must  there¬ 
fore  be  incorporated  into  the  chamber  pressurization  analysis  given  in 
Sec.  III.  A.  That  analysis  as  presented  assumes  a  step  increase  in 
propellant  flowrate  at  time  zero. 

A  rigorous  theoretical  treatment  of  the  problem  of  transient  flow 
increase  upon  valve  opening  involves  flash  vaporization,  feed  system 
dynamics,  complex  heat  transfer  problems ,  two  phase  flow,  etc.  and 
is  beyond  the  scope  of  the  present  six-months  program.  Nevertheless 
the  transient  flow  period  is  treated  empirically  for  the  specific  feed 
system/propellant  valves /injector  (1-on-l)  assembly  used  during  the 
program. 

As  a  first  approximation,  the  mass  flowrate  is  assumed  to  increase 
linearly  with  time  until  the  nominal  steady  state  flow  is  reached.  Then, 
in  Eq.  1  of  Sec.  III.  A,  the  mass  of  propellant  Wj  injected  during  one 
time  interval  Atjq  becomes: 
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and 


Wj  =  WggA^j^L- 
wj  =  w8SAtN 


if  0  —  t  *=  t 


88 


tSB**  i^°° 


where  tse  is  the  time  to  reach  the  nominal  steady  state  flowrate, 


(21a) 

(21b) 

^68- 


Since  the  mass  of  an  ensemble  (j  values)  of  drops  injected  in  one 
time  interval  must  remain  constant  (i.  e.  mass  of  drops  plus  mass  of 
vapors  generated  from  these  drops  is  constant)  as  the  ensemble  is 
"operated  on"  in  succeeding  time  intervals  (N  values),  Eq.  21a  is 
modified  for  use  in  the  computer  program  giving: 


W:  =  w86  (AtN)Z  ( if  0  <  N  < (JiltL  +  0)  (22) 

\  ^88  I  A^N 

and  (N  -  -^)<  )'  <  N 
AtN' 

where  9  -  drop  residence  time  expressed  in  number  of  time  intervals 
and  j'  =  N  -j  +  1.  (It  is  more  convenient  for  the  machine  calculations 
to  use  j'-values  rather  than  the  j-values  as  defined  in  Sec.  III.  A). 

The  effect  of  using  Eqs.  22  and  21b  in  place  of  W:  =  w  Atj^  of  Eq.  1 
is  shown  in  Fig.  15  for  and  MMH  in  Engine  Configuration  No.  1. 

The  two  solid  curves  in  the  figure  are  obtained  from  the  equations 
given  in  Sec.  III.  A  applied  to  Engine  Configuration  No.  1  (Table  I) 
with  the  experimental  conditions  given  in  Table  IV.  The  equations 
assume  a  step  function  increase  in  flowrate. 

The  two  dashed  curves  in  Fig.  15  are  obtained  using  Eqs.  22  and 
21b  to  describe  the  flow  increase.  Values  of  tsS  were  obtained  by 
curve  fitting,  1.  6  msec  and  0.  8  msec  (from  propellant  entry)  for 
N2O4  and  MMH.  respectively,  giving  adequate  agreement  between  the 
initial  portions  of  the  experimental  and  computed  chamber  pressuriza¬ 
tion  curves.  The  discrepancies  between  the  latter  portions  of  the 
experimental  and  computed  curves  are  due  to  wall  heating  effects  which 
are  described  in  the  following  section. 

D.  WALL  HEATING  EFFECTS  ON  PRE -IGNITION  CHAMBER 

PRESSURIZATION 

From  comparisons  of  experimental  and  calculated  chamber  pressuriza¬ 
tion  histories  for  singly  flowed  propellants,  it  was  ror«cluued  in  the 
previous  program  that  heat  is  transferred  from  the  chamber  walls  to 
the  vaporizing  vapor/drop  system  during  an  ignition  delay  period  (1). 

The  added  energy  leads  to  significantly  higher  vaporization  chamber 
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pressures  than  would  result  from  an  adiabatic  vapor/drop  system 
(Fig.  15). 

Jn  the  present  program,  three  possible  mechanisms  for  pre- 
ignition  heat  transfer  between  injected  propellants  and  thrust 
chamber  walls  have  been  considered:  radiant  heating,  convective 
heating  of  propellant  vapor,  and  conductive  heating  of  droplets 
upon  impingement  with  chamber  walls.  Because  of  the  low  tempera¬ 
tures  involved  (0°  to  85°F),  radiant  heating  is  insignificant.  Also 
convective  heating  of  propellant  vapors  is  relatively  minor,  amounting 
to  only  a  5  to  10°F  gaa  temperature  increase,  depending  on  engine 
configuration.  The  temperature  increase  was  calculated  by  con¬ 
ventional  methods  for  convective  heat  transfer  in  a  tube  (6). 

The  problem  of  conductive  heating  of  propellant  droplets  imping¬ 
ing  upon  chamber  walls  is  complex.  A  rigorous  theoretical  treatment 
would  take  into  account  drop  shattering  and  the  effects  thereon  of 
various  impingement  angles  and  velocities,  gradual  wetting  of  the 
chamber  walls  and  the  resulting  "insulation”  effect,  contact  area 
between  drops  and  the  walls,  varying  total  drop  surface  area  over  which 
vaporization  occurs,  varying  wail  surface  temperature,  etc.  Such  a 
ttwuWwu.,  ox  course,  is  beyond  the  scope  of  the  present  program.  The 
effects  of  conductive  droplet  heating  have  nevertheless  been  accounted 
for  by  empirically  determining  time -dependent  coefficients  for  the 
time-varying  drop  temperatures  of  the  three  drop-size  spray  model. 

The  time  dependent  characteristic  of  the  coefficients  is  required  to 
account  for  the  experimentally  determined  fact  that  chamber  wall 
surfaces  cool  during  the  propellant  vaporization  period  (Fig.  13). 

Drop  heating  is  caused  to  begin  in  the  computer  program  once  the 
initially  entering  drops  have  been  in  the  system  long  enough  to  travel 
to  the  chamber  walls.  Once  a  drop  reaches  the  walls,  its  temperature, 
which  has  been  and  continues  to  be  calculated  bv  Eq  16  of  Sec.  III.  A, 
is  increased  by  a  fraction,  f,  of  the  '•mperature  difference  between 
the  chamber  wall  and  the  theoretical  new  drop  temperrture,  i.  e. 

^Tdij^N  ~  <Tdij>N  +  f  (Tw  “  ^dij^J  (Z3) 

The  fraction  f  is  time -dependent  initially  to  reflect  cooling  of  the  chamber 
walls  (T^.  in  Eq.  23  is  held  constant  at  the  initial  hardware  temperature). 

The  fraction  f  therefore  starts  at  a  high  value,  f _ and  decreases  with 

time  to  a  lower  value.  f  j  ,  which  thereafter  remains  constant  i.  e 
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The  factors  required  by  Eqs.  23  and  24  to  empirically  account 
for  wall  heating  are  therefore  f  x,  tmjn,  t^  (the  time  interval 
during  which  f  goes  from  fmax  to  fmjn).  and  Djw  (the  pathlength  a 
drop  traverses  from  its  injector  orifice  to  impingement  on  the 
chamber  wall).  The  former  three  are  obtained  empirically  by 
curve  fitting;  the  latter,  Diw,  is  calculated  from  the  chamber 
geometry  and  the  propellant  injection  angle. 


The  effects  of  using  Eqs.  23  and  24  to  account  for  wall  heating 
are  shown  in  Figs,  lb  and  1  7  for  and  MMH  in  Engine  Configu¬ 

rations  1  and  2.  The  computed  curves  take  into  account  the  appro- 
Driate  engine  geometry  (Table  I)  and  the  conditions  of  the  appro¬ 
priate  experimental  runs  (Table  IV).  The  figures  show  reasonable 
agreement  between  the  experimental  chamber  pressurization  curves 
fr  r  each  propellant  singly  flowed  in  eacn  engine  and  the  corresponding 
calculated  pressurization  curves  which  include  wall  heating  as  well 
as  the  iinearized  flow  increase  discussed  above.  The  empirically 
determined  values  for  the  coefficients  and  factors  for  i  'hiding  wall 
heating  anu  linearized  flow  increase  in  the  calculated  curves  of  Figs. 

1  6  and  !  7  are: 

Engine  Configuration  No.  1  Engine  Configuration  No.  2 


N:°i 

MMH 

n2o4 

MMH 

t  s  g .  msec 

! .  6 

0.  8 

1.  6 

0  8 

^max 

0.  15 

0.  1  5 

0.  25 

0.  25 

*min 

0.  02 

O 

O 

O 

c 

0  03 

ldh *  m8Ct 

2.  5 

2.  5 

2.  5 

2  5 

The  calculations  show  that  many  N^O^  drops  evaporatively  cool  to 
the  freezing  point,  11.  8°F.  and  then  progressively  solidify  as  evapora¬ 
tion  continues.  In  contrast,  none  of  the  MMH  drops  reaches  the  freezing 
temperature.  -6 2°F.  In  the  drop  heating  calculations,  those  N2C>4  drops 
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Figure  17,  Experimental  and  Analytical  Chamber  Pressure  Histones 
for  Singly  Flowed  (top)  and  MMH  (bottom)  in  Engine 

Configuration  No.  2. 
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having  a  frozen  fraction  X-  (Eq.  19)  of  less  than  0.  6  when  they 
first  reach  the  chamber  wail  are  first  melted,  upon  impingement 
with  the  wall,  and  then  their  temperatures  are  augmented  (by  Eqs. 

23  and  24)  in  subsequent  time  intervals.  Those  NpO^  drops  having 
a  frozen  fraction  of  0.  6  or  greater  at  impingement  with  the  chamber 
walls  are  assumed  to  ricochet  off  the  walls  without  being  warmed. 

The  calculated  curves  which  include  drop  heating  and  linearized 
flow  increase  satisfactorily  reproduce  the  chamber  pressurization 
histories  of  the  singly  flowed  propellants.  Thus,  the  calculations 
satisfactorily  predict  the  dashed  curves  of  Figs.  11  and  12  for 
Engine  Configurations  1  and  2,  respectively.  In  addition,  the 
calculations  give  the  gas  temperature  history  in  each  engine  the 
vapor  mass  histories  of  each  propellant  in  the  cnambers,  time- 
dependent  propellant  vaporization  rates,  vapor  flowrates  through 
*he  nozzles,  drop  temperature  histories,  etc.  The  next  step  in 
the  development  of  the  mathematical  model  for  hypergolic  ignition 
delay  in  space-ambient  engines  is  to  account  for  the  disparity  be¬ 
tween  the  dashed  curves  of  Figs.  11  and  12  and  their  corresponding 
experimental  pre -ignition  curves  for  hot  firings  given  in  the  same 
Wo  figures.  This  step  requires  an  understanding  of  the  pre-ignition 
reactions  which  occur  during  an  ignition  delay  period  with  these 
propellants.  The  next  section  of  the  report  is  concerned  with  these 
pre-ignition  reactions,  the  reaction  intermediate  formed,  and  the 
method  by  which  to  incorporate  formation  of  the  reaction  intermediate 
into  the  ignition  delay  model. 
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SECTION  IV 

PHASE  H  -  CHEMICAL  KINETICS  ASPECTS  OF  PRE*  IGNITION 
CHAMBER  PRESSURIZATION 

In  the  case  of  nitrogen  tetroxide /hydrazine-type  fuels  as  propellants, 
a  condensed  phase  reaction  product  has  been  observed  during  engine 
testing  at  simulated  altitudes.  This  section  of  the  report  is  concerned 
with  the  role  of  this  reaction  product  in  hypergolic  ignition  delay.  The 
general  experimental  conditions  under  which  the  reaction  product  has 
been  observed  are  briefly  discussed  first.  Then,  the  explosive 
characteristics  of  the  material  <?  determined  from  standard  hazard 
evaluation  tests  are  descril  ^a.  Text,  some  physical  properties  of 
th^  material  are  given,  j  allowed  b--  discussions  of  the  composition  of 
the  liquid  reaction  product,  its  heat  J  formation,  and  its  formation 
kinetics.  Finally  in  thi^  section  of  die  report,  the  manner  in  which 
formation  of  the  reaction  product  affects  ignition  delay  time  is  analyzed. 

A.  OCCURRENCE  AND  APPEARANCE  OF  SUB -IGNITION  REACTION 
INTFRMEriAT  S 

i.  Laboratory  Occurrences 

The  overall  kinetic  factors  of  ignition  reactions  of  ^O^/hydra- 
zine-type  fuels  were  experimentally  determined  at  Thiokol-RMD  during 
a  previous  program  on  hypergolic  ignition  at  reduced  pressures  (Contract 
No.  AF04(£l  l}-9946.  Part  U).  It  was  found  that  a  condensed  phase 
reaction  product  formed  from  vapor-phase  reactants  under  conditions  of 
room  temperature  and  low  (sub-ignition)  pressures  (1).  The  material 
appear  ed  as  a  clear,  yellow,  viscous  liquid  of  low  vapor  pressure.  Sub¬ 
ignition  reaction  intermediates  were  found  for  N2O4  in  combination  with 
MMH,  50-50,  UDMH  and  N2H4.  The  quantity  of  material  formed  depended 
on  the  fuel:  MMH  and  50-50  giving  approximately  ten  times  as  much  as 
UDMH,  and  N2H4  yielding  only  very  small  amounts  under  the  conditions 
of  the  experiments. 

The  occurrence  of  a  reaction  intermediate  in  the  N7O4/N2H4  system 
has  been  observed  elsewhere  under  very  different  experimental  conditions 
(7,  8).  Using  infrared  absorption  spectroscopy,  Weiss  monitored  low 
temperature  reactions  of  solid-phase  N2O4/N2H4  at  atmospheric  pressure. 
Upon  gradual  warming  of  the  reactants  from  an  initial  temperature  of 
liquid  nitrogen,  Weiss  found  one  reaction  with  little  heat  evolution 
occurring  at  temperatures  as  low  as  -20l°F.  A  further  reaction,  in¬ 
volving  measurable  heat  evolution,  occurred  at  about  -75°F.  The  spectrum 
of  the  higher  temperature  reaction  was  found  to  agree  well  with  that  of 
hydrazinium  nitrate. 
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In  experiments  conducted  by  the  Bureau  of  Mines  (8),  mixtures 
of  granulated  N2O4/N2H4  at  liquid  nitrogen  temperature  were  per¬ 
mitted  to  warm  at  atmospheric  pressure.  Eventually,  while  the 
mixture  still  appeared  to  be  frozen,  a  violent  reaction  occurred  and 
often  a  yellow,  viscous,  liquid  residue  remained.  From  IR  spectra 
of  the  residue  and  other  spot  tests  it  was  concluded  that  the  residue 
contains  an  appreciable  percentage  of  hydrazinium  nitrate  which  is 
a  moderately  unstable  explosive. 

While  the  work  of  Weiss  and  of  Christos  et  al  corroborate  the 
existence  of  important  reaction  intermediates  in  N2O4/N2H4  systems, 
the  initial  conditions  of  the  reactants  in  their  experiments  differ 
greatly  from  those  in  an  engine.  In  Section  III  of  this  report,  it  is 
found  that  the  fuel  (MMH)  does  not  freeze  during  the  ignition  delay 
period  and  that,  although  the  smaller  droplets  Jo  freeze,  the 

minimum  oxidizer  drop  temperatures  at  near  0°F.  Also,  for  a 
propellant  feed  temperature  of  75°F,  the  minimum  average  gas 
temperature  in  a  typical  engine  is  about  25°F.  Finally,  the  gas 
pressure  increases  from  near  zero  of  space  to  100-300  mmHg, 
depending  on  engine  configuration,  during  the  ignition  delay  period. 

2.  Deduction  from  Experimental  Engine  Tests 

Indications  of  appreciable  pre -ignition  chemical  reactions  involving 
substantial  heat  release  are  found  in  Figs.  11  and  12.  In  these  figures, 
experimental  chamber  pressure -time  histor'es  of  singly  flowed 
propellants  in  Engine  Configurations  1  and  2  ar-’  plotted  together  with 
hot  firing,  pre-ignition  Pc-time  curves.  The  dotted  curve  in  each 
figure  is  simply  the  sum  of  .lit  pressure  curves  of  each  propellant 
singly  flowed.  Good  agreement  exists  between  the  sum -of -partial- 
pressures  curves  and  the  pre-ignition  hot  firing  curves  for  the  first 
millisecond  of  the  rising  pressures.  A  gradual  divergence  then  befins, 
the  hot  firing  pre -ignition  curves  ris-ng  rror  3  steeply  than  the  sum-of- 
partial  pressure  curves.  For  both  engine  configurations,  the  slopes 
of  the  two  hot  firing  pre-ignition  curves  increase  with  time  until 
ignition  occurs.  The  divergences  which  .each  100-130  mmHg  at 
ignition,  are  evidence  of  appreciable  pre-igni*jon  reactions 

A  further  indication  of  pre-ignition  reactions  comes  from  the 
chamber  wall  ter.  oerature  data  given  in  Section  111,  B.  2.  The 
oscillograms  of  Figs.  13  and  14  show  that  wall  temperatures  decrease 
in  the  case  of  singly  flowed  propellants  but  that,  in  a  hot  firing,  wall 
temperatures  begin  to  increase  early  in  the  ignition  delay  period. 


-53- 


(SNOraW)  H10N3T3AVM 


■uliilifiii 

■  ■■■■■■buhi 
pHiiiaamgi 

SKSIS! 

mmKMBsam 

HnvssnpHn 

RBssssmsisisi 

wpnsHri 

■  fifes  till 

■  M|$f§ 

jlHHiSSii'aSiiiBS 

B55b88b:s«!S{ 

EiSBBHSgSSiiiBi 


sssssssRissiiiSii 

tlSSKMKeSSHilliPHli 


SBRH8I 

■NKWI 


■B88I 

!b:s: 

■BSgSHSSSHIS! 


2  §32 


EH 

Mil 


jBBBSSSiiSSSSiiS! 


[■■SgSrsiss 
i!SBSSS:8SSSSSS8SiS! 
SjBMBMB88iia 
EgfilsSiSiESili'B! 

■■flKSI 


iSlaSSSSISS 

M 

IMHK1 

IBSSSSBSSSiB 


■Rj|Bg 

Biiiiii 


ilij 

IBS 


o  «  Q 


33NVWOSAV 


Figure  18.  IR  Spectra  of  Liquid-Phase  Reaction  Intermediates  of  N2O4/MMH 
Produced  in  an  Engine  (top)  and  in  a  Laboratory  Apparatus  (bottom). 
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3.  Engine-Produced  vs.  Laboratory-Produced  Reaction  Inter¬ 
mediates. 

IR  absorption  spectra  of  engine-produced  and  laboratory -produced 
reaction  intermediates  of  N2O4/MMH  are  given  in  Fig.  18.  The 
spectrum  of  the  laboratory -produced  material  is  taken  from  Reference 
1.  The  liquid-phase  sample  was  produced  from  an  undiluted,  stoichio¬ 
metric  mixture  of  the  reactants  in  a  pyrex  flow  system  at  low  pressure 
(6  mmHg)  and  room  temperature  (Sec.  IV.  A.  4,  below). 

The  engine-produced  material  and  its  spectrum  were  obtained  under 
another  program  at  Thiokol-RMD  (9).  The  100-lb  thrust  engine  was 
operated  in  a  pulse  mode  at  a  simulated  altitude  of  approximately 
125,  000  ft.  The  pulsing  rate  was  2  hertz  with  an  electrical  pulse  width 
of  10  msec.  hardware  and  propellant  feed  temperatures  were  60°F. 

Under  these  conditions  of  very  short  pulse  width  and  relatively  long 
off-time,  ignition  did  not  occur  during  the  160  pulses  of  the  test  sequence. 
Consequently,  sufficient  material  was  available  at  the  conclusion  of  the 
pulse  train  to  permit  collection  and  limited  analysis. 

As  shown  in  Fig.  18,  the  spectrum  of  this  sample  is  very  similar 
to  that  of  the  laboratory-produced  material.  The  only  significant 
difference  concerns  the  absorptions  between  1250  and  1500  cm"*. 

In  this  region,  the  laboratory  material  exhibits  one  broad  absorption 
band,  characteristic  of  a  nitrate.  In  contrast,  the  engine-produced 
material  has  a  weak  absorption  band  at  about  1325  cm"*  plus  another 
band  at  about  1460  cm"*. 

Other  features  of  the  two  spectra  are  essentially  identical,  including 
the  absorption  band  at  1200  cm”*.  This  band  is  generally  attributed  to 
nitrite;  however,  *io  nitrite  has  been  detected  by  wet  chemical  analyses 
of  the  laboratory-produced  material  (see  below,  Section  IV.  D.  1). 

It  should  be  pointed  out  that  the  laboratory-produced  sample  was 
protected  from  environmental  contamination  at  all  times  by  collecting 
and  maintaining  the  sample  under  a  dry  nitrogen  atmosphere.  In 
contrast,  ambient  air  was  used  to  bring  the  pressure  to  one  atmosphere 
in  the  22-ft  diameter  vacuum  sphere  before  the  engine  sample  could  be 
withdrawn.  Consequently,  some  contamination  of  the  engine  sample  could 
have  occurred.  Evidence  in  support  of  this  is  given  below  in  Sections  IV.  B.  1 
and  3.  Nevertheless,  despite  some  sample  contamination  or  degradation, 
the  IR  spectrum  of  the  engine-produced  sample  is  very  similar  to  the 
spectrum  of  the  laboratory-produced  material.  The  latter  is  therefore 
assumed  to  be  the  same  as  the  fresh  material  produced  in  an  engine  during 
an  ignition  delay  period. 
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4.  Production  of  Material  for  Analyses  and  Determinations 

The  N2O4/MMH  reaction  intermediate  used  for  chemical 
analyses  and  property  determinations  in  the  present  program  was 
produced  from  gas-phase,  undiluted  reactants  at  low  (sub -ignition) 
pressure  and  room  temperature.  The  experimental  apparatus  was  used 
initially  in  the  previous  hypergolic  ignition  program  and  is  fully  de¬ 
scribed  in  the  final  report  of  that  investigation  (1).  The  apparatus  is 
shown  in  Figs.  19  and  20  of  the  present  report,  the  latter  figure  being 
a  schematic  diagram.  Briefly,  the  apparatus  consists  of  vaporizing 
and  flowmetering  sections,  a  flow  reactor  and  a  low  pressure  hou°’ng 
with  observation  windows. 

Gaseous  N2O4  was  taken  directly  from  a  No.  3  shipping  cylinder 
(99.  5  mole  %  minimum  purity,  Matheson  Co.  ,  l^c.  )  immersed  in  a 
constant  temperature  bath  at  30°C.  The  heat  source  was  sufficient 
to  maintain  steady  flow  conditions  of  up  to  0.  2  gr/sec.  The  gaseous 
stream  was  metered  in  a  Fischer-Porter  2F-1/4-20-5  flowmeter 
(glass  float)  at  30°C  and  600  mmHg  pressure.  Under  these  conditions, 
the  oxidizer  exists  as  a  mixture  of  NO2  and  N2O4.  Therefore,  the 
flowmeter  was  calibrated  in  terms  of  mass  flow  by  collecting  the 
vapors  in  a  condenser  at  liquid  nitrogen  temperature  and  weighing 
the  mass  collected  during  measured  time  intervals. 

The  fuel,  MMH,  was  flowmetered  in  the  liquid  state  (08F-1/16- 
20-4/36,  sapphire  float)  and  then  vaporized  drop-wise  in  a  glass 
coil  heated  to  90°C.  The  vapors  expanded  into  a  one -liter  flask, 
which  was  also  heated  to  90°C,  and  then  passed  through  an  all-glass, 
heated  line  to  the  flow  reactor. 

The  flow  reactor  shown  schematically  in  Fig.  20  is  a  47-mm 
diameter  pyrex  tube  opened  at  one  end  to  the  ambient  pressure  of  the 
low  pressure  vessel.  The  flow  tube  is  housed  inside  a  metal  tube 
which  extends  downward  from  the  l  ft  diameter  x  3  ft  high  vacuum 
vessel  (Fig.  19).  The  pressure  in  the  latter  is  regulated  both  by  a 
throttle  valve  in  the  connecting  line  to  a  high  capacity  vacuum  pump 
(Kinney  KD  7  80)  and  by  an  appropriate  bleed  of  shroud  gas  (nitrogen) 
which  has  the  additional  function  of  diluting  the  corrosive  gases. 

The  procedure  for  making  the  reaction  intermediate  is  as  follows. 
The  apparatus  is  first  evacuated.  Then  a  bleed  of  nitrogen  shroud  gas 
of  about  1  CFM  is  admitted  into  the  vessel  and,  by  adjusting  the 
throttling  valve  in  the  vacuum  line,  a  pressure  of  about  0.  5  mmHg  is 
obtained  in  the  vessel.  The  desired  flowrates  of  fuel  and  oxidizer 
are  then  set.  The  ambient  pressure  in  the  chamber  is  then  gradually 
raised  by  bleeding  room  air  into  the  vacuum  line  through  an  auxiliary, 
downstream  valve.  If  the  pressure  is  raised  sufficiently  (5  to  IS  mmHg 
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depending  upon  reactant  mixture  ratio,  initial  temperature  and  flow 
tube  diameter)  ignition  occurs  at  a  well-defined,  reproducible, 
minimum  pressure. 

If  the  pressure  is  held  below  the  minimum  ignition  pressure 
smoke  is  seen  issuing  from  the  flow  tube.  Some  of  the  condensed 
phase  material  deposits  onto  the  walls  of  the  flow  tube  and  eventually 
collects  as  a  viscous  yellow  liquid  at  the  base  of  the  vertically  - 
oriented  flow  tube.  Six  to  eight  cc  of  material  are  collected  at  a 
time. 


After  shutting  off  the  propellant  flows,  the  vacuum  vessel  and 
flow  tube  are  filled  to  one  atmosphere  pressure  with  dry  nitrogen. 

The  sample  is  then  removed,  maintaining  a  nitrogen  atmosphere 
over  the  liquid  at  all  times. 

Two  N2O4/MMH  mixture  ratios  were  used  in  making  samples 
for  the  various  analyses  and  determinations  performed  during  the 
program.  The  stoichiometric  mixture  (O/F  =  2.  5  by  wt,  )  and  an 
oxidizer-rich  mixture  (O/F  -  5.0)  were  used,  the  latter  is  within 
the  composition  range  of  the  vapor -phase  reactants  in  an  engine 
during  an  ignition  delay  period  (O/F  of  simultaneously  injected 
liquid  propellants  of  1.  6).  The  lean  vapor-phase  mixture  results 
from  the  higher  volatility  of  the  oxidizer. 

B.  SENSITIVITY  CHARACTERIZATION  OF  REACTION  INTERMEDIATE 

Because  of  the  unknown  composition  and  characteristics  of  the 
N2O4/MMH  reaction  intermediate  at  the  outset  of  the  program,  the 
initial  tests  performed  on  the  material  were  of  the  hazard  evaluation 
type.  The  purpose  of  the  tests  was  to  signal  any  need  for  extra¬ 
ordinary  handling  procedures  and,  as  a  by-product,  to  provide  at 
least  qualitative  information  on  the  role  of  the  reaction  intermediate 
in  ignition  pressure  spiking. 

The  material  used  in  these  tests  was  made  from  the  stoichiometric 
mixture  of  N2O4/MMH. 

1.  Impact  Sensitivity 

The  impact  sensitivity  of  a  test  sample  is  obtained  in  drop 
weight  tests  by  dropping  a  known  weight  from  various  pre-set  heights 
onto  a  carefully  designed  cup  containing  the  material.  A  loud  report, 
puncture  of  a  retaining  steel  diaphragm,  flame  or  evidence  of  combus¬ 
tion  are  all  taken  as  signs  of  a  positive  test.  The  drop  height  which 
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yields  a  50%  probability  of  initiation  is  determined.  The  50%  height 
for  the  N2O4/MMH  reaction  intermediate  tested  is  10  inches  using 
a  2  kg  weight.  Table  VIII  gives  results  for  several  materials  for 
comparative  purposes. 

A  portion  of  the  engine-produced  sample  whose  IR  spectrum  is 
given  in  Fig.  18  was  subjected  to  drop  weight  tests  under  a  separate 
program  (9).  The  material  was  not  sensitive  in  the  drop  weight 
tester  at  a  drop  height  of  36  inches,  the  maximum  height  of  the 
apparatus.  The  stability  of  the  material  is  attributed  to  sample 
contamination  and/or  degradation,  discussed  above. 

For  comparative  purposes,  monomethyl  hydrazine  nitrate 
(CH3N2H3- HNO3)  was  prepared  in  the  laboratory  (9).  The  synthesis 
route  involved  the  direct  reaction  of  MMH  in  aqueous  solution  v/ith  a 
stoichiometric  quantity  of  HNO3  followed  by  dehydration  and  re¬ 
crystallization.  The  IR  spectrum  of  the  material  is  given  in  Fig.  21. 

TABLE  VIII 

DROP  WEIGHT  AND  MODIFIED  TRAUZL  BLOCK 
TEST  VALUES  FOR  VARIOUS  MATERIALS 


Material 

Drop  Weight 
50%  Height 
inch 

Trauzl  Block 
Expansion* 
cc  /  gr 

Nitroglycerir 

1 

27 

Hydrazine  nitrate,  N2H4-  HNO3  4 

20 

RDX  (explosive) 

7 

25 

98%  H202 

- 

15 

Lead  Azide,  Pb(N3)2 

3** 

10 

N^O^/MMH  Reaction 

Intermediate  10 

9 

ch3n2h3-  HNO3 

1  5 

- 

nh4cio4 

18 

0 

*  Corrected  for  volume  increase  due  to  blasting  cap  alone. 
♦♦Picatinny  Arsenal  Apparatus,  2  kg  weight  (10). 
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The  50%  drop  height  for  the  synthesized  compound  was  found  to  be 
15  inches,  indicating  that  the  laboratory -produced  N£0^/MMH 
reaction  intermediate  is  more  sensitive  to  impact  than  MMH*  HNO3. 

2.  Explosive  Strength 

Modified  Trauzl  Block  Tests  conducted  routinely  at 
Thiokol-RMD  use  a  sample  size  smaller  than  in  standard,  full- 
size  Trauzl  Block  tests.  Typically,  one  to  three  gram  samples 
are  used  in  a  lead  cylindrical  cup  of  given  dimensions.  A  standard 
No.  8  blasting  can  is  used  for  initiation.  Detonation  of  the  blasting 
cap  in  contact  with  the  Bample  causes  an  increase  in  internal 
cylinder  volume.  Corrected  for  the  blasting  cap  effect,  the  volume 
change  is  related  semi -quantitatively  to  the  explosive  strength  of 
the  test  sample.  The  modified  Trauzl  Block  tests  on  the  reaction 
intermediate  gave  a  corrected  expansion  of  9  cc/gr  for  a  2.  5  gram 
sample.  This  value  is  close  to  that  of  lead  azide  tested  in  the  same 
apparatus  (Table  VIII). 

3.  Thermal  Stability 

The  standard  ICRPG  thermal  stability  test  was  performed 
on  a  sample  of  reaction  intermediate  produced  in  the  laboratory. 

The  heating  record  is  shown  in  Fig.  22  together  with  similar  curves 
of  engine-produced  material  and  CHjN^Hy-  HNO3  synthesized  from 
MMH  and  HNO3  (9).  In  each  case,  the  heating  bath  was  programmed 
to  increase  at  a  rate  of  8°  C  per  minute.  The  curves  are  a  record 
of  the  differential  temperature  (Tsarrpje  -  Tbath)  plotted  against  the 
temperature  of  the  sample.  Because  of  the  polarity  of  the  recorder, 
positive  peaks  represent  exothermic  reaction.  The  initial  dip  in  each 
curve  is  due  to  the  lag  in  heating  of  the  sample  when  the  bath  heater 
is  activated. 

The  thermogram  of  the  reaction  intermediate  produced  in  the 
laboratory  indicates  that  two  exothermic  reactions  occur:  the  first 
starts  at  about  60°C.  the  second  at  about  110°C.  The  rate  and 
intensity  of  the  second  reaction  are  greater  than  that  of  the  first. 
i.l  a  sample  temperature  of  1  35°C  following  the  strong  exothermic 
reaction,  the  pressure  was  sufficient  to  burst  a  rupture  diaphragm 
rated  at  4700  psi. 

The  heating  curve  of  the  engine-produced  material  also  shows 
two  exotherms;  however,  they  are  displaced  to  slightly  higher 
temperatures  indicating  somewhat  greater  thermal  stabilstv  than 
the  laboratory -produced  material.  The  greater  stability  is  attri¬ 
buted  to  sample  contamination  during  retrieval  (Sec.  IV.  A.  3).  In 
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addition,  :he  curves  of  the  two  samples  show  that  the  intensities 
of  the  two  exothermic  reactions  are  reversed,  indicating  that  the 
reactions  are  not  sequential  reactions  of  the  same  material  but 
rather  are  due  to  two  materials  present  in  different  proportions 
in  the  two  samples- 

Th~  heating  curve  of  monomethyl  hydrazine  nitrate  prepared  in 
the  laboratory  is  significantly  different  from  those  of  the  reaction 
intermediates  (Fig.  22).  MMH'  HNO3  begins  to  decompose  exo- 
thermally  at  190°C.  The  greater  thermal  stability  of  this  compound 
over  that  ot  the  reaction  intermediate  does  not  preclude  the  com¬ 
pound  from  being  a  constituent  of  the  intermediates.  However,  the 
curves  of  Fig.  22  indicate  that  MMH1  HNOj  is  not  the  cause  of  the 
two  exotherms  at  lower  temperatures. 

C.  PHYSICAL  PROPERTIES  OF  REACTION  INTERMEDIATE 


Measurements  of  a  few  physical  properties  of  the  reaction 
intermediate  formed  by  stoichiometric  N^O^/MMH  in  the  laboratory 
apparatus  were  attempted:  liquid  density,  vapor  pressure,  freezing 
point  and  molecular  weight. 

1.  Liquid  Density 


The  density  of  the  material  was  determined  at  three 
temperatures.  The  results  obtained  are: 

T emp.  ,  °C  Density,  gr/cc 


1. 205 


17.  8 


1.  189 


25.  0  1.  180 

2.  Vapor  Pressure 

Vapor  pressure  measurements  of  the  material  at  temperatures 
from  0°  to  40°C  were  attempted  using  standard  techniques.  However, 
values  could  not  be  obtained  due  to  gradual  but  significant  decom¬ 
position  of  the  test  sample  throughout  the  temperature  range  investi¬ 
gated.  Although  the  decomposition  was  slow,  enough  gas  evolved  to 
obscure  the  actual  (low)  vapor  pressure  of  the  sample. 

Frjm  the  laboratory  production  runs  of  the  reaction  intermediate, 
ii  could  be  determined  that  the  vapor  pressure  of  the  material  is  list 
than  0.  5  mmHg  at  room  temperature.  This  value  is,  of  course,  well 
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below  the  vapor  pressures  of  the  reactants  (760  and  41  mmHg  at 
70°F  for  N2O4  and  MMH  respectively). 

Since  it  is  believed  that  the  reaction  intermediate  also  forms 
during  tail-off  of  a  pulse,  it  is  desirable  to  know  the  vapor  pressure 
of  the  material  at  various  temperatures  for  considerations  of  its 
removal  between  pulses  by  the  pumping  action  of  a  space  environ¬ 
ment,  In  this  regard,  although  the  decomposition  of  the  sample 
was  sufficient  to  preclude  accurate  vapor  pressure  determinations 
by  conventional  techniques,  on  a  gross  scale  the  decomposition  is 
very  slow.  For  example,  at  the  end  of  one  month,  the  remainder 
of  a  batch  of  material  which  had  been  kept  in  a  dry  box  qualita¬ 
tively  underwent  little  change  in  volume,  although  the  original 
clear  yellow  color  of  the  liquid  did  change  to  nearly  water-white. 

3.  Freezing  Point 

The  shape  of  a  cooling  curve  of  an  unknown  substance 
can  provide  information  on  whether  the  material  is  a  pure  compound 
or  a  mixture  of  compounds.  For  this  reason,  freezing  point 
measurements  of  the  laboratory  produced  reaction  intermediate 
were  attempted;  they  were  unsuccessful,  however,  due  to  the  very 
high  viscosity  of  the  material.  In  general,  very  viscous  materials 
tend  to  greatly  supercool  and,  when  this  occurs,  accurate  and 
reliable  freezing  point  data  cannot  easily  be  obtained.  In  the 
present  case,  freezing  points  ranging  from  -20  to  -50°C  were 
observed  but  the  data  are  not  deemed  reliable  due  to  the  viscosity 
effect.  For  example,  at  -10°  to  -20°C,  the  material  "poured"  a 
distance  of  one  inch  in  about  ten  minutes. 

4.  Molecular  Weight 

An  attempt  to  determine  the  molecular  weight,  or  average 
molecular  weight,  of  the  reaction  intermediate  by  the  freezing  point 
suppression  technique  in  benzene  as  solvent  was  unsuccessful.  It 
was  found  that  the  reaction  intermediate  is  not  soluble  in  benzene. 
Other  solvents  were  not  tested. 
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D.  COMPOSITION  OF  REACTION  INTERMEDIATE 

To  determine  the  composition  of  the  N2O4/MMH  reaction 
intermediate,  a  number  of  chemical  and  instrumental  analyses 
were  performed.  The  liquid  samples  were  produced  from  vapor 
N2O4  and  vapor  MMH  at  sub -ignition  pressures  (5-9  mmHg)  and 
room  temperature  in  the  laboratory  apparatus.  Two  mixture 
ratios  of  the  reactants  were  used  in  making  the  samples: 

O/F  =  2.  5  and  5.  0  by  wt.  The  former  is  the  stoichiometric, 
mixture  ratio  and  the  latter  is  within  the  composition  range  of 
the  vapor  phase  reactants  in  an  engine  during  an  ignition  delay 
period  (O/F  of  simultaneously  injected  liquid  propellants  of  1.6). 
The  lean  mixture  results  from  the  higher  volatility  of  the  oxidizer. 

The  initial  analyses  of  the  reaction  intermediates  were  made 
over  a  period  of  a  few  weeks  following  production  of  the  samples. 
Because  of  gradual  sample  decomposition,  which  was  learned 
subsequently  from  the  attempted  vapor  pressure  determinations 
(Sec.  IV.  C.  2),  the  results  of  the  initial  analyses  reflect  changes 
in  composition  due  to  sample  degradation.  These  results,  which 
suggest  that  the  composition  of  the  reaction  intermediate  depends 
upon  the  mixture  ratio  of  the  initial  reactants,  are  discussed 
subsequently. 

The  results  of  analyses  given  in  the  next  paragraphs  were 
obtained  within  two  days  of  production  of  the  sample  and  are  taken 
as  representative  of  freshly  produced  material  in  an  engine. 

1.  Results  of  Analyses  of  Fresh  Sample 

A  number  of  chemical  and  instrumental  analyses  were 
made  of  a  fresh  sample  of  the  N2O4/MMH  reaction  intermediate 
and  on  the  residue  auid  condensate  of  fractionation  experiments. 
The  sample  on  which  the  analyses  were  performed  was  made  at 
low  pressure  from  N2O4/MMH  at  a  mixture  ratio  of  5.  0  by  wt. 

The  infrared  spectrum  of  the  sample  is  given  in  Fig.  23.  It 
is  similar  to  the  spectrum  obtained  under  the  previous  program 
(1),  and  reproduced  in  Fig.  18,  except  that  free  MMH  is  indicated 
in  Fig.  23  by  the  shapes  of  the  absorptions  between  1150  and  800 
cm-1. 
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A  summary  of  the  results  of  the  chemical  analyses,  which 
were  performed  within  two  days  of  sample  production,  is  given 
in  Table  IX.  For  the  elemental  analyses,  the  percentages  of 
carbon,  hydrogen  and  nitrogen  were  determined;  the  values 
given  are  averages  of  two  determinations.  The  balance  of  the 
sample  is  assumed  to  be  oxygen. 

As  given  in  the  table,  the  specie  analyses  account  for  64.  1% 
of  the  sample.  However,  it  is  reasonable  to  assume  that  the 
measured  N-N  exists  as  MMH,  NH3  as  CH3  NH2  and  NO3"  as 
HNO3.  Adjusted  to  reflect  these  compounds,  the  specie  analyses 
then  account  for  85.  6%  by  wt.  ,  as  given  in  Table  X.  Approxi¬ 
mately  15%  of  the  sample  is  therefore  unaccounted.  No  nitrite 
(N02~)  and  no  CO2  were  detected  in  the  freshly  prepared 
sample. 

The  water  content  of  the  sample  (Tables  IX  and  X)  was 
determined  by  gas  chromatographic  analysis.  It  was  noted  during 
this  analysis  that  the  sample  also  contained  free  MMH.  The  amount 
of  the  latter,  however,  was  not  determined  quantitatively. 

Two  long  term  pumping  experiments  were  performed  on 
samples  of  the  fresh  material.  The  samples  were  maintained 
at  room  temperature  and  any  volatiles  removed  during  pumping 
were  trapped  at  liquid  nitrogen  temperature. 

In  one  experiment,  the  sample  was  pumped  on  for  40  working 
hours  at  a  pressure  of  2  mmHg.  The  loss  in  weight  of  the  sample 
was  0.  238  gr,  or  2.  51%  of  the  initial  weight.  IR  analysis  of  the 
collected  liquid  condensate  indicated  pure  MMH. 

Analysis  of  the  residue  after  pumping  gave:  C,  16.  53%  by 
wt.  ;  H,  9.  13%;  and  N,  40.  64%.  IR  indicated  that  the  residue  was 
similar  to  the  original  material  but  with  less  free  MMH.  Gas 
chromatographic  analysis  of  the  residue  indicated  no  significant 
change  in  H2O  content  from  the  original  material. 

A  second  sample  of  the  reaction  intermediate  was  pumped  on 
for  24  working  hours  at  a  pressure  of  0.  1  mmHg.  The  condensate 
and  residue  were  examined  by  IR  (Fig.  24).  The  condensate  was 
found  to  be  predominantly  MMH  but  with  some  nitrate  present  as 
indicated  by  N03*  absorption  at  1350  cm-1.  The  spectrum  of  the 
residue  is  similar  to  that  obtained  from  the  previous  pumping 
experiment. 
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TABLE  IX 

SUMMARY  OF  CHEMICAL  ANALYSES  OF  REACTION 
INTERMEDIATE  OF  N204/MMH  (C/F  =  5.  0  by  wt)* 


Elemental  Analyses  (duplicates): 
C 
H 
N 

O  (by  difference) 


Total 


1 6.  30  -  0.  09%  by  wt. 

8.60  i  0.10% 

42.79  *  0.95% 

32.  3 1% _ 

1 00.  0%  by  wt. 


Specie  Analyses 

Total  N-N  bond  functionality  (duplicate) 

Total  Atnmoniacal  functionality  (NH-j) 
(duplicate) 

NOj"  (duplicate) 

no2" 

CO-> 

ft# 

H20  (by  gas  chromatography) 

Total 


32.  75  -  0.  04%  by  wt. 

0. 24  i  0.  02% 

15.81  i  0.18% 
none  detected 
none  detected 

15.3% _ 

64. 1 0%  by  wt. 


♦All  analyses  performed  within  two  days  of  production  of  saxnpla. 
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TABLE  X 


ADJUSTED  SPECIE  ANALYSES  OF  N204/MMH 
REACTION  INTERMEDIATE 


Total  N-N  bond  functionality  as  MMH 
Total  ammoniacal  functionality  as  CH3  NH2 
N03‘  as  HNO3 

h2o 


co2 

no2 


53.  8 1%  by  wt. 

0.44% 

1  6.  08% 

15.3  % 
none  detected 
none  detected 


Total: 


85.  63%  by  wt. 
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Figure  Z4.  IR  Spectra  of  Condensate  (top)  and  Residue  (bottom)  of  Low  Pressure 
Fractionation  of  N2O4/MMH  Reaction  Intermediate. 
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Attempts  were  made  to  remove  water  from  a  sample  of  the 
original  material  by  drying  with  molecular  sieves  (Type  4A). 
Apparently,  the  water  is  tied  up  strongly,  for  it  was  not  removed 
by  this  method. 

2.  Comparison  of  Measured  and  Computed  Elemental 
Analyses 

To  gain  insight  into  the  15%  of  the  reaction  intermediate 
unaccounted  for  by  the  specie  analyses,  an  elemental  analysis 
was  computed  from  the  specie  analyses  given  in  Table  X.  The 
computed  elemental  analysis  is  compared  to  the  measured 
analysis  in  Table  Xi.  The  table  indicates  that  the  various  specie 
analyses  do  not  account  for  significant  amounts  (>19%)  of  carbon, 
nitrogen  and  oxygen.  Other  species  to  account  for  the  discrepan¬ 
cies  have  not  been  identified  to  date. 

3.  Possible  Composition  Derived  from  Specie  Analyses 

A  possible  composition  of  the  reaction  intermediate  can  be 
deduced  from  the  specie  analyses  as  follows.  First,  it  is  assumed 
that  all  of  the  CH3NH2  ( Table  X)  exists  in  the  reaction  intermediate 
as  CH3NH3NO3.  Next,  the  balance  of  HNO3  is  assumed  to  exist  as 
MMH‘  HNO3.  Finally,  all  of  the  H2O  exists  as  MMH-  H2O.  The 
derived  composition  of  the  intermediate  is  given  in  Table  XII.  The 
derived  composition  accounts  for  only  as  much  of  the  reaction 
intermediate  as  do  the  specie  analyses,  which,  as  stated,  are  the 
basis  of  the  derived  composition.  However,  despite  the  still 
unaccounted  15%,  the  composition  given  in  Table  XII  indicates  a 
small  but  significant  amount  (3.  2%)  of  free  MMH-  In  corroboration, 
the  IR  spectrum  of  the  original  material  indicates  free  MMH  aa  a 
constituent  and  free  MMH  was  also  indicated  by  the  gas  chromato¬ 
graphic  analysis.  Furthermore,  in  one  of  the  long  term  pumping 
experiments  (Sec.  IV.  D.  1),  the  condensed  volatiles  were  found  to 
be  pure  MMH  and  the  loss  in  weight  of  the  original  material  amounted 
to  2.  5%,  a  figure  which  compares  favorably  with  the  amount  of  free 
MMH  in  the  derived  composition  (Table  XU). 


COMPARISON  OF  MEASURED  ELEMENTAL  ANALYSIS  WITH 
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TABLE  XU 

REACTION  INTERMEDIATE  COMPOSITION  DERIVED 
FROM  ADJUSTED  SPECIE  ANALYSES^ 


Constituent 

ch3nh3no: 

MMH-  HN03 
MMH-  H20 
Free  MMH 

Total: 


Derived  Composition 
of  Intermediate 
by  wt. 

1.33% 

26.  28% 

54.  8  % 

3.  22% 

85.  63  % 


*  Table  X 

Despite  the  reasonableness  of  the  derived  composition  with 
regard  to  the  analytical  results,  an  anomaly  persists  concerning 
vapor  pressure.  The  vapor  pressures  of  MMH  and  MMH-  H,0 
are  very  similar,  being  about  41  mmHg  at  70°F.  Although  the 
vapor  pressures  of  MMH  HN03  and  CH3NH3NO^  are  not  known, 
they  are  expected  to  be  quite  low.  Assuming  that  the  latter  are 
negligible  and  that  no  complicating  molecular  associations  occur 
among  the  four  compounds,  the  vapor  pressure  of  the  four- 
component  mixtrre  can  be  loosely  estimated  by  Raoult’s  Law  and 
is  found  to  be  32  mmHg  at  70°F.  In  the  computation  which 
utilizes  the  derived  composition  given  in  Table  XII,  MMH  and 
MMH' H^O  are  treated  as  a  single  solvent  (after  computing  their 
mole  fractions)  since  their  vapor  pressures  are  the  same. 

The  very  great  discrepancy  between  the  estimated  vapor 
pressure  of  the  reaction  intermediate  and  the  experimentally 
determined  approximate  value  of  0.  5  mmHg  (Sec.  IV.  C-  2)  suggests 
for  one  thing  that  there  is  considerable  molecular  association  in 
the  mix'-ire  Cf  course.  Raoult’s  La  •>  assumes  equal  sizef  for  all 
molecules  and  this  condition  is  clearly  not  met.  However,  it  is 
considered  unlikely  that  the  size  effect  alone  accounts  for  the  nearly 
two  orders  of  magnitude  encountered  here. 


-73- 


AFRPL-TR-67-129 


Despite  the  insights  gained,  the  composition  of  the  freshly 
prepared  N2O4/MMH  reaction  intermediate  is  not  yet  completely 
resolved.  The  specie  analyses  suggest  that  the  material  is  a 
mixture  of  four  compounds:  MMH-  H^O,  MMH*  HNOj,  free  MMH 
and  CH3NH3NO3.  In  corroboration  of  this,  the  existence  of  free 
MMH  has  been  independently  substantiated  by  instrumental 
analyses  (IR  and  GC)  and  by  fractionation  experiments. 

The  four  compounds,  however,  account  for  only  85%  of  the 
reaction  intermediate.  The  balance  of  the  material  is  as  yet 
unidentified,  although  the  elemental  and  specie  analyses  indicate 
that  it  contains  about  14*.  by  wt.  C,  43%  N  and  43%  O.  This  gives 
an  overall  empirical  formula  for  the  unidentified  15%  of  approxi¬ 
mately  C4N|()Og. 

It  is  apt  ireni  at  this  point  that  the  composition  of  the  freshly- 
produced  reaction  intermediate  is  complex.  The  situation  is 
further  complicated,  however,  by  the  dependence  of  composition 
of  the  reaction  intermediate  on  mixture  ratio  of  the  initial  re¬ 
actants,  and  MMH.  This  conclusion  is  reached,  a*  will  be 

discussed,  based  on  the  results  of  analyses  which  follow  in  the 
next  sub-section. 

4.  Gradual  Decomposition  Effects  and  Composition  Dependency 
Upon  Initial  Reactant  Mixture  Ratio 

The  chemical  analyses  performed  Ir.i^aUy  in  the  program 
spanned  a  several  week  period  following  production  of  the  sample* 
Because  01  gradual  sample  decomposition,  which  was  subsequently 
learned,  the  results  of  the  initial  analyses  do  not  describe  the 
composition  of  freshly  prepared  material. 

The  chemical  analyses  discussed  in  the  preceding  sections 
were  made  within  two  days  of  sample  production;  the  ones  to  be  dis¬ 
cussed  in  this  section  aie  the  initial  analyses  made  on  "old"  samples. 
These  analyses  are  of  interest  for  two  reasons.  First,  a  compari¬ 
son  of  the  results  that  fol.ow  with  those  given  above  gives  some 
insight  into  the  effects  of  the  gradual  decomposition  of  the  material. 
Secondly,  since  the  initial  analyses  were  performed  on  two  samples 
made  from  two  different  mixture  ratios  of  N2O4/MMH,  the  degree 
of  constancy  of  composition  of  the  reaction  intermediate  with 
initial  reactant  mixture  ratio  can  be  ascertained. 
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The  results  of  the  analyses  of  both  samples  are  summarized  in 
Table  XIII.  In  the  specie  analyses,  the  total  N-N  bond  functionality, 
ammoniacal  functionality  and  NO3’  are  given  in  terms  of  MMH,  CH3NH2 
and  HNO3  respectively,  as  in  Table  X  for  the  fresh  sample  analyses. 

The  most  striking  difference  between  TabJes'XIH  and  X  concerns 
CO^  none  was  detected  in  the  fresh  sample  whereas  significant  amounts 
we '3  found  in  the  degraded  samples. 

It  is  most  probable  that  the  CO2  in  the  degraded  samples  exists  as 
methylhydrazinium-N-methylcarbazate,  CH^N^^COO^H^CH^.  Other 
constituents  and  their  concentrations  in  the  two  reaction  intermediates  can 
then  be  deduced  as  follows:  The  CH^NH^  (Table  XIII)  is  assumed  to  exist 
in  the  degraded  samples  (as  well  as  in  the  fresh  sample,  Sec.  IV.  D.  3 
above)  as  CH3NH3NO3.  Then,  the  balance  of  the  HNG3  is  assumed  to 
exist  as  HNO3.  Finally,  the  remainder  of  the  MMH,  after  allowing 

for  that  in  the  CO2  containing  compound  above,  exists  as  CHj^Hj’  Zl-^O. 

The  test  as  to  the  adequacy  of  this  model  for  the  compositions  of 
the  degraded  reaction  intermediates  lies  in  whether  the  "remaining  MMH" 
in  the  di-hydrated  state  satisfactorily  completes  the  material  balances 
between  the  derived  compositions  and  the  specie  analyses.  As  shown  in 
Table  XIV,  the  agreement  is  very  good.  The  derived  compositions, 
which  consist  of  the  same  four  compounds  for  each  degraded  reaction 
intermediate,  are  well  within  experimental  error  of  the  specie  analyses. 

It  is  to  be  noted  that  the  amount  of  H2O  obtained  by  the  specie  analyses 
is  just  sufficient  for  the  "remaining  MMH"  to  exist  as  MMH1  2H2O.  That 
is,  the  ratio  of  H2O  ■  to  -  "remaining  MMH"  is  the  same  for  both  samples 
despite  significant  differences  in  constituent  concentrations  (Table  XIV). 

Although  good  accountability  of  the  results  of  the  specie  analyses 
is  achieved  by  the  four-component  model  of  intermediate  composition, 
there  remains  an  unidentified  7-8%  portion  of  the  degraded  samples.  By 
comparison,  the  specie  analyses  of  the  fresh  sample  leave  an  unidentified 
15%  (Table  X).  Nevertheless,  because  of  the  consistency  of  the  derived-though- 
incomplete  analyses  of  the  degraded  samples  and  the  independent,  experi¬ 
mental  corroboration  of  one  constituent  in  the  fresh  sample,  it  seems 
reasonable  that  the  unidentified  portion  of  each  sample  is  a  fifth  (and  sixth?) 
constituent  and  will  not  nullify  the  constituent  compounds  deduced  above. 
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TABLE  XIII 

RESULTS  OF  CHEMICAL  ANALYSES  OF  DEGRADED  SAMPLES 
OF  N204/MMH  REACTION  INTERMEDIATES 


Elemental  Analyses: 

C 

H 

N 

O  (by  difference) 


Sample  A 

(N204/MMH  =  2.5) 

1 7.  4%  by  wt. 

9.4% 

32.  6% 

40.  6% 


Total  1 00.  0% 


Sample  B 

(N204/MiviH  =  5.0) 

15.5%  by  wt. 

8.5% 

30.  7% 

45.  3% _ 

1  00.  0% 


Specie  Analyses: 

Total  N-N  bond 
functionality  (computed 

in  terms  of  MMH)  50.  1%  by  wt.  49.  7%  by  wt. 


NO3"  (computed  in 

terms  of  HNO^)  12.5%  12,4% 


NO- 


noue  detected 


Total  ammoniacal 

functionality  (computed 

in  terms  of  CH3NH2)  0.  88% 

CC^  4. 4% 

h2°  25%  (t  3%) 


Total  92. 9% 


none  detected 


1.  00% 
8.4% 

20%  (t  3%) 
91.  5% 
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Although  the  differences  between  the  results  of  the  elemental 
analyses  of  the  two  degraded  samples  are  small  (Table  XIII),  they  are 
larger  than  experimental  error.  If  degradation  effects  are  similar  in 
the  two  samples,  then  the  original  compositions  of  the  samples  were 
different.  The  implication,  then,  is  that  the  reaction  intermediate 
composition  depends  upon  initial  reactant  mixture  ratio.  Furthermore, 
since  the  product  is  apparently  a  mixture  rather  than  a  single  compound, 
it  is  likely  that  product  composition  depends  also  upon  pressure.  Con¬ 
sequently,  the  reaction  intermediate  formed  early  in  au  ignition  delay 
period  may  be  different  from  the  material  formed  late  in  the  delay 
period  due  to  both  changing  reactant  mixture  ratio  and  pressure.  Then, 
since  shock  sensitivity,  thermal  stability,  detonability,  etc.  depend  on 
composition,  ignition  pressure  spikes  are  also  affected.  Further  work 
is  therefore  required  to  determine  more  completely  the  dependency  of 
product  composition  on  initial  reactant  mixture  ratio  and  on  pressure 
and  to  evaluate  resulting  changes  in  sensitivity  characteristics. 

E.  Heat  of  Formation  of  Reaction  Intermediate. 

One  requirement  for  incorporating  the  formation  of  the  reaction 
intermediate  into  the  fundamental  hypergolic  ignition  delay  model  is  the 
heat  evolved  when  the  material  forms.  The  heat  of  reaction  for  formation 
of  the  reaction  intermediate  can  be  calculated  if  the  heat  of  formation  of 
the  material  is  known.  The  heat  of  formation  was  therefore  determined 
experimentally. 

The  material  used  for  the  determination  was  made  at  low  pressure 
from  N2O4/MMH  at  a  mixture  ratio  of  5.  0  by  wt.  The  material  was  taken 
from  the  same  production  batch  as  the  samples  used  for  the  chemical 
analyses  reported  in  Se~  IV.  D.  1.  The  heat  of  formation  determination 
was  also  performed  within  a  few  days  of  production  of  the  material. 

1.  Method  of  Determination 

The  determination  was  made  by  measuring  the  heat  of  combustion 
of  the  sample  in  the  presence  of  pure  oxygen.  A  Parr  oxygen-bomb  adiabatic 
calorimeter  was  used.  It  was  calibrated  according  to  a  standard  NBS 
procedure  using  NBS  benzoic  acid  with  a  known  heat  of  combustion  of  6318 
cal/gr. 


Six  runs  were  made  in  determining  the  heat  of  combustion  cf  the 
reaction  intermediate.  Corrections  to  the  measured  heats  of  combustion 
were  applied  for  the  heat  liberated  upon  ignition  of  a  standard  10  cm  length 
of  fuse  wire  and  for  nitric  acid  formed  during  the  combustion.  The  nitric 
acid  correction  was  determined  by  titrating  the  washings  with  a  prepared 
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solution  of  Na£C03  where  1  ml  of  titrant  equalled  one  calorie.  Carbon 
dioxide  liberated  during  the  combustion  was  determined  by  a  COj>  traih 
of  anhydrone  and  preweighed  ascarite  absorption  tubes.  The  results 
are  reported  in  terms  of  %  C  recovered. 

2.  Experimental  Results 

The  calorimeter  constant  was  found  to  be  2443  cai/C°.;  ' 

The  experimental  data  of  the  six  runs  for  the  heat  of  combustion 
of  the  reaction  intermediate  are  given  in  Table  XV.  As  shown  in  the 
table(  the  average  heat  of  combustion  is  3.591  kcal/gr. 

From  the  elemental  analysis  of  the  sample  (Table  IX),  an 
empirical  formula  for  the  reaction  intermediate  is  found  to  be 
C  H  N  O  for  an  arbitrary  molecular  weight  of  1  gr. 

0*0136  0*0860  0*0306  0*0202 

Then,  an  equation  for  the  combustion  reaction  of  the  material  in  the 
presence  of  oxygen  is: 

C  M  N  O  (1)+J  O  (g)  =  K  HO  (1)  +L  CO  (g)+MN  (g) 

for  which  AHC,  the  measured  heat  of  combustion,  is  -3.591  kcal/gr. 

The  coefficients  K  and  L  in  the  equation  are  found  to  be,  in  units  of  mass, 
0.744  gr  and  0.598  gr,  respectively,  for  one  gram  of  reaction  intermediate 

The  heat  of  formation  of  the  reaction  intermediate  is  then  obtained 
from  AHfpr0(jUC£B  “AHf  reactants  wbere  AHC  =  -3.591  kcal  /gr 

and  the  heat  of  formation  of  H2O  (1)  and  CO2  (g  )  are  -3. 793  kcal  /  gr  and 
-2. 137  kcal  /gr,  respectively.  The  heat  of  formation  of  the  liquid  reaction 
intermediate  is  therefore  -0.509  kcal  /gr.  A  "per  unit  mass"  basis  has 
been  used,  rather  than  the  conventional  "per  mole"  basis,  because  the 
molecular  weight  of  the  material  is  unknown.  The  latter  can  be  postulated, 
however,  as  shown  below. 
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TABLE  XV 

EXPERIMENTAL  DATA  ON  HEAT  OF  COMBUSTION  OF 
N204/MMH  (O/F  =  5.0)  REACTION  INTERMEDIATE 


Run 

Sample 

Weight 

%  C  Recovered 
as  CO2 

A  H  Combustion 

1 

gram 

0.4269 

— 

cal/gr 

3686 

2 

0.7659 

— 

3672 

3 

1.0223 

15.35 

3679 

4 

1.2139 

14.  66 

3465 

5 

1.0140 

14.87 

3606 

6 

1.0630 

14.79 

3440 

Ave rage 


14.  92 


3591 
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3.  Heat  of  Reaction  for  Formation  jf  Intermediate 

For  calculating  the  heat  released  due  to  formation  of  the  reaction 
intermediate  during  an  ignition  delay  period,  the  overall  stoichiometry  of 
the  formation  reactions  is  required.  In  the  abser  ce  of  additional  informa¬ 
tion,  it  is  necessary  to  estimate  the  stoichiometry  as  follows. 

From  the  elemental  analysis  of  the  material  (Table  IX),  an 
empirical  formula  is  found  to  be  CH^  32N2.  25^1  486'  or  ve,y  nearly 
^4^25^9®6"  Since  the  hydrogen  and  ~arbon  must  come  from  MMH,  the 
H-to-C  ratio  must  be  6:1  assuming  no  gaseous  product  involving  either 
of  these  elements  is  forried.  Thus,  an  empirical  formula  of  C^H^NgO^ 
satisfies  the  H-to>C  rrdo  and  is  very  close  .o  the  empirical  formula  given 
by  the  elemental  analysis.  The  overall  stoichiometry  for  formation  of  this 
specie  is  given  by* 

3  N02(g)  +  4  CH3N2H3  (g)-»C4H24N906(l)  +  N2(g). 

The  heat  of  this  reaction  is  calculated  to  be  A  Hr=  -66.  2  kcal 
per  mole  of  MMH  using  the  value  of  the  heat  of  formation  of  the  inter¬ 
mediate  determined  above  and  7.  91  kcal/mole  and  22.  8  kcal/mole  for 
the  heats  of  formation  of  N02  (g)  and  MMH  (g)  respectively. 

For  complete  combustion 

10  NC2(g)  +  4  CH3N2H3(g)  4  C02(g)  +  12  H20(g)  4  9N2(g) 

for  which  the  heat  of  reaction  is  A  Hr=  -310.  0  kcal  per  mole  of  MMH. 

For  the  stoichiometry  given  above,  ti.e  heat  released  upon  formation  of  the 
intermediate  is  21%  of  the  total  heat  of  complete  combustion,  on  a  fuel 
basis. 

F.  KINETICS  OF  FORMATION  OF  REACTION  INTERMEDIATE 

Overall  formation  kinetics  are  also  required  for  incorporating  the 
reaction  intermediate  into  the  fundamental  hypergoiic  ignition  delay  model. 
An  acbabatic,  low  pressure,  flow  tube  system  was  used  in  which  the 
reactant  vapors  were  greatly  diluted  with  nitrogen.  The  dilution  technique 
was  used  in  order  to  provide  known  gas  velocities  and  stream  properties, 
particularly  heat  capacity,  along  the  flow  tube.  The  temperature  rise  due 
to  reaction  along  the  flow  tube  axis  was  measured. 

1.  Experimental  Apparatus  and  Instrumentation 

The  apparatus  consists  of  propellant  vaporizing  and  flowmetering 
sections  end  an  instrumented  flow  tube  which  exhausts  into  a  low  pressure 
housing.  The  apparatus  is  shown  photographically  in  Fig.  25  and 
schematically  in  Fig.  26. 
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Figure  25.  Apparatus  tor  Study  of  Formation  Kinetics  of 
Reaction  Intermediate.  Fuel  Vaporizing  and  Fiow.sr.ete ring 
Sections  (top)  and  Instrumented  Flow  Tube  (bottom). 
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Schematic  Diagram  of  Apparatus  for  Formation  Kinetics  Study 
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The  oxidizer  vaporizing  and  flowmetaring  section  of  the 
apparatus  were  used  for  both  the  kinetics  study  and  the  sample 
production  runs.  The  oxidizer  feed  system  is  described  in  Sec. 

IV.  A.  4.  The  oxidizer  transport  line  to  the  flow  tube  was  wire 
wound  to  permit  preheating  of  the  vapor  for  runs  at  elevated 
initial  temperatures.  The  gaseous  oxidizer  stream  was  pre-mixed 
with  the  diluent  nitrogen  stream  prior  to  entering  the  flow  tube. 

In  a  new  fuel  feed  system  built  for  the  kinetics  study.  MMH 
vapor  was  taken  directly  from  a  300  cc  stainless  steel  cylinder 
immersed  in  a  6S°C  temperature  bath  (Fig.  26).  For  the  highest 
flows,  two  cylinders  connected  in  parallel  were  used  in  the  single 
temperature  bath.  Steady  vapor  flows  of  up  to  0.  1  gr/sec  could 
be  maintained.  The  gaseous  stream  was  metered  by  a  calibrated 
nozzle  sized  to  produce  choked  flow  at  the  desired  flowrates.  The 
pyrex  nozzle  was  jacketed  permitting  hot  water  (65°C)  to  be  circu¬ 
lated  continuously.  Except  where  jacketed,  all  fuel  feed  lines  were 
wrapped  with  nichrome  heating  wire  (Fig.  25)  in  order  to  prevent 
condensation  of  the  fuel  vapors. 

Downstream  from  the  metering  nozzle,  a  two-way  stopcock 
directed  the  fuel  stream  either  to  the  flow  tube  or  into  a  line  which 
bypassed  the  flow  tube.  The  latter  was  used  while  setting  the  desired 
fuel  flowrates  (below).  When  an  established  fuel  flow  was  shifted  to 
the  flow  tube,  the  vapor  passed  first  into  a  manifold  (near-end  of  flow 
tube.  Fig.  25)  which  fed  eighteen  pyrex  tubes  of  2  mm  bore.  The  tubes 
distributed  the  fuel  vapors  across  the  inlet  of  the  flow  tube  (Fig.  27) 
to  promote  rapid,  complete  mixing  of  the  fuel  with  the  premixed  nitrogen/ 
oxidizer  stream. 

Diluent  nitrogen  gas  (Extra  Dry  Grade,  Matheson  Co.  Inc.  ) 
was  flowmetered  at  9.  8  psig  and  25°C  in  a  Fischer-Porter  FP-1/2-27- 
GiO/80  flowmeter  with  an  SVT-44  float.  After  expansion  to  low 
pressure,  the  nitrogen  was  heated  to  the  desired  initial  temperature  by 
flowing  the  gas  through  copper  coils  in  a  large  constant  temperature  water 
bath. 


The  pyrex  flow  tube,  1. 93  in.  diameter  x  48  in.  long,  has  six 
fine-wire  thermocouples  located  along  the  length  of  the  tube  with  the 
beads  positioned  at  the  tube's  centerline.  The  thermocouples  were  made 
from  Pt-10%  Rh  and  Pt  wires  0.  002  in.  in  diameter.  Uncoated  thermo¬ 
couple  beads  were  used,  based  on  experience  gained  in  the  previous 
program.  In  that  program,  temperature  measurements  were  made  in 
undiluted,  low  pressure  streams  of  N2O4/MMH  with  coated  and  uncoated 
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Pt-10%  Rh/Pt  thermocouples.  No  catalytic  effects  were  found  (1). 

The  cold  junctions  of  the  thermocouples  in  the  flow  tube  were  immersed 
in  a  single  ice  bath.  The  signals  of  any  four  of  the  thermocouples 
were  amplified  by  CEC  DC  Amplifiers  and  recorded  simultaneousl^- 
by  a  CEC  oscillograph.  The  system  provided  a  sensitivity  of  12C 
per  inch  deflection. 

The  flow  tube  was  wire  wound  (Fig.  25)  and  insulated  for 
runs  at  elevated  initial  temperatures.  The  flow  tube  exhausted  into 
the  low  pressure  housing  and  vacuum  system  discussed  in  Sec.  IV. 

A.  4  (Fig.  25). 

A  summary  of  the  experimental  parameters  and  their  ranges 
is  given  in  Table  XVI. 


TABLE  XVI 

SUMMARY  OF  EXPERIMENTAL  PARAMETERS  ANP  THEIR 
RANGES  IN  THE  KINETICS  STUDY 


Total  pressure,  mmHg 

20  - 

60 

Total  flowrate,  gr/sec 

1.  73 

-  3. 

26 

Reactants  to  total  flow  ratio,  wt. 

0.  05 

-  0. 

13 

Reactant  mixture  ratio.  O/F  by  wt. 

0.  58 

-  5. 

i 

Initial  temperature,  °C 

23  - 

75 

Flow  tube  diameter,  in. 

i.  93 

Flow  tube  length,  in. 

48 

Flov  velocity,  ft/sec 

35  - 

170 

ricw  Mach  number 

less 

than 

0.  15 

Reynolds  number  (diameter) 

2800 

-  4800 

Under  the  conditions  of  the  experiments,  the  minimum  equili¬ 
brium  percent  dissociation  of  ;he  oxiiizer  is  95*%  according  to  data  repo*  led 
in  (11).  Because  the  dissociation  and  assoc  ation  reactions  are  extremely 
fast,  the  equilibrium  mixture  is  a*tained  prior  to  oxidiaer  entry  into  the 
flow  tube  Further  dissociation  into  nitric  oxide  does  not  occur  at  the 
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temperatures  encountered  in  the  experiments  (80°C  maximum).  In 
view  of  the  above,  the  oxidizer  is  taken  as  being  NO^  only. 

The  general  procedure  for  making  the  kinetic  measurements 
at  room  temperature  was  as  follows.  Afte*  evacuating  the  apparatus, 
the  desired  nitrogen  diluent  flow  was  established.  The  pressure  in 
the  flow  tube  was  then  adjusted  to  the  desired  level  using  the  throttle 
valve  in  the  vacuum  line  downstream  from  the  low  pressure  housing 
(Sec.  IV.  A.  4).  Then,  the  signals  fr.m  the  thermocouples  were 
recorded  to  check  for  temperature  uniformity  along  the  flow  tube. 

Next,  the  desired  oxidizer  flow  was  established  and  the  pressure 
in  the  flow  tube  readjusted  if  necessary.  The  temperature  profile 
along  ihe  tube  was  again  recorded  to  check  for  uniformity.  Then, 
utilizing  the  bypass  line,  the  desired  fuel  flow  was  established. 

With  the  fuel  flow  set  and  the  recorder  on,  the  fuel  flow  was  shifted 
to  the  flow  tube  by  quickly  rotating  the  two-way  stopcock  in  the  fuel 
feed  system. 

From  the  continuously  recorded  thermocouple  signals, 
it  was  found  that  steady  state  conditions  were  attained  within  a  few 
seconds  after  shifting  the  fuel  flow  into  the  flow  tube.  The  short 
time  to  steady  state  conditions  is  due  to  the  low  temperature  rise 
(generally  5  to  10C°)  from  the  reactions  in  the  diluted  system. 

For  runs  at  elevated  initial  temperatures,  the  procedure  was 
the  same  although  pre-heating  of  the  flow  tube  and  the  three  gaseous 
streams  and  their  lines  was  of  course  necessary.  Considerable  care 
was  required  to  obtain  Variac  and  heater  settings  which  would  provide 
uniform  inlet  temperatures  for  the  three  streams  and  also  a  uniform 
temperature  profile  along  the  tube  axis  when  the  diluent  and  either 
reactant  were  flowed. 

2.  Ex;  erimental  Results 

A  total  of  68  runs  were  made  in  the  flow  tube  apparatus, 
including  the  initial  check-out  tests.  The  ranges  of  experimental 
parameters  covered  in  the  tests  are  given  in  Table  XVI. 

The  output  data  was  recorded  in  the  form  of  thermocouple 
signals  versus  real  time.  The  thermocouple  signals  were  converted 
to  temperature  by  manually  measuring  the  deflection  of  the  traces. 

As  a  rule,  temperature  values  were  obtained  from  the  records  when 
(1)  nitrogen  alone  was  flowing,  (2)  nitrogen  plus  oxidizer  were  flowing, 
and  (3)  nitrogen  plus  oxidizer  plus  fuel  were  flowing  at  a  point  in  time 
following  attainment  of  steady  state  condition  (usually  2  seconds). 
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For  the  cases  of  nitrogen  and  nitrogen  plus  oxidizer,  temperatures 
were  reduced  to  check  for  uniformity  of  the  axial  temperature 
profile  without  reaction. 

From  the  known  flowrates,  temperature,  pressure  and 
flow  tube  diameter,  gas  velocities  were  computed  for  each  run. 

Then,  using  the  computed  velocities  and  the  known  thermocouple 
positions  relative  to  the  final  mixing  point  of  the  reactants  (Fig. 

27),  the  time  from  mixing  to  the  thermocouple  stations  was  computed. 
The  data  is  presented  below  in  the  form  of  temperature  rise  due  to 
reaction  versus  time  from  mixing. 

Figures  28  and  29  present  results  of  similar  runs  at  room 
temperature  and  elevated  temperature  for  two  reactant  mixture 
ratios.  In  terms  of  complete  reaction  to  CO £,  and  Ng,  Fig. 

28  is  for  a  fuel-rich  mixture  and  Fig.  29  is  for  an  oxidizer-rich 
mixture. 


The  most  striking  feature  of  the  two  figures  is  the  small 
temperature  rise  of  the  high  temperature  runs  as  compared  to  the 
large  temperature  rise  of  the  room  temperature  runs.  The  signifi¬ 
cance  of  this  is  discussed  in  detail  in  the  following  sub- section, 
Discussion  of  Results.  In  short,  the  large  temperature  rise  of 
the  -oom  temperature  runs  includes  condensation  phenomena  which 
do  not  occur  in  the  high  temperature  runs. 

It  should  be  noted  in  the  two  figures,  particularly  the  high 
temperature  runs,  that  the  temperature  rises  immediately  from  the 
mixing  point  and  that  the  temperature  rise  increases  linearly  with 
time  during  most  of  the  stay  time  of  the  reactants  in  the  flow  tube. 

The  implications  of  this  behavior  are  three-fold.  First,  mixing 
is  accomplished  very  quickly.  Second,  reaction  starts  immediately 
without  any  induction  period.  This  point  is  important  with  regard  to 
the  mechanism  of  reaction  as  it  indicates  a  pure  thermal  mechanism, 
without  chain  branching.  Thirdly,  the  linear  temperature  rise 
indicates  that  only  small  amounts  of  reactants  are  consumed.  The 
local  reactant  concentrations,  then,  are  very  nearly  the  initial 
reactant  concentrations.  Advantage  of  this  point  is  taken  in  inter¬ 
preting  the  data  in  the  next  section. 

In  all  of  the  high  temperature  runs  and  in  some  room 
temperature  runs  (those  at  the  lower  pressures),  the  temperature 
rises  do  not  exceed  5C°.  Thus,  the  temperature  difference  between 
the  gaues  and  the  flow  tube  wall  is  everywhere  small  in  these  cases, 
making  the  system  very  nearly  adiabatic.  For  the  room  temperature 
runs  at  the  higher  pressures,  the  excessive  temperature  rises  lead 


reaction 


AFRPL-TR- 67-129 


Time  from  Mixing,  msec 


Figure  29.  Effect  of  Pressure  and  Initial  Temperature  on  Temperature 
Rise  due  to  Reaction  of  NO2/MMH  (O/F  =  3.  5) 
in  Flow  Tube  with  Diluent  Nitrogen. 
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to  non-adiabaticity  which  tends  to  cause  a  downward  bending  of  the 
A  T-time  curves. 

The  results  of  room  temperature  runs  at  four  reactant  mixture 
ratios  are  given  in  Figs.  30  through  33.  The  concentration  of  the 
reactant  not  in  excess  (in  terms  of  complete  combustion)  is  con¬ 
stant  for  all  runs  at  the  same  pressure.  The  interrelationship  of 
the  reactant  concentrations  in  these  tests  is  given  below,  where 
X  is  the  concentration  of  the  propellant  not  in  excess: 

Concentration 
MMH  NOz 


30 

0.  58 

1.7  X 

X 

31 

1.  0 

X 

X 

32 

2.  5 

X 

2.  5  X 

33 

3.  5 

X 

3.  5  X 

3.  Discussion  of  Results  ‘ 


no2  ' 

Figure  No.  MMHi 


A  comparison  of  the  runs  at  elevated  initial  temperatures  with 
similar  runs  at  room  temperature  (Figs.  28  and  29)  suggests  on  the 
surface  a  negative  temperature  coefficient.  This  is  not  the  case, 
however,  since  an  additional  process  was  operative  in  the  low 
temperature  runs  only.  Visual  observations  of  the  flow  tube  re¬ 
vealed  that  condensation  occurred  in  the  room  temperature  runs 
but  not  in  the  runs  at  elevated  initial  temperature.  Condensation, 
when  it  occurred,  was  readily  apparent  as  the  flow  tube  became 
quickly  clouded.  This  was  the  case  in  all  room  temperature  runs 
but  in  none  of  the  elevated  initial  temperature  runs.  The  temperature 
rise  of  the  lower  temperature  runs  includes,  then,  the  heat  of  con¬ 
densation  of  the  reaction  intermediate. 

Considering  first  the  runs  at  elevated  initial  temperatures  where 
condensation  effects  are  absent,  the  overall  order  of  reaction  for 
formation  of  vapor  phase  reaction  intermediate  can  be  determined 
from  the  effect  of  pressure  on  temperature  rise  due  to  reaction. 

For  a  volume  element  of  reactants  in  an  adiabatic  system,  the 
heat  equation  is: 


/5cp  ^  =  QW 


(25) 
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N02/MMH  =  0.  58 
Tc  =  25°C 
Reactants 


Total  Flow 


30  mmHg 


Time  from  Mixing,  msec 

Figure  30.  Effect  of  Pressure  on  Temperature  Rise  due  to  Reaction 
of  N02/MMH  (O/F  *  0.  58)  in  Flow  Tube  with  Diluent  Nitrogen. 


N02/MMH  =1.0 
T0  =  24°C 
Reactants 


Total  Flow 


30  mmHg 


Time  from  Mixing,  msec 

Figure  31.  Effect  of  Pressure  on  Temperature  Rise  due  to  Reaction 
of  N02/MMH  (O/F  =  1. 0)  in  Flow  Tube  with  Diluent  Nitrogen. 
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Time  from  Mixing,  msec 


F18U/kt^2'/ww»^  Pres8ure  on  Temperature  Rise  due  to  Reaction 
of  N02/MMH  (O/F  =  2.  5)  in  Flow  Tube  with  Diluent  Nitrogen. 


t 
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Figure  33.  Effect  of  Preeiure  on  Temperntur'e  Riee  due  to  Reaction 
of  NO2/MMH  (O/F  =  3.  5)  in  Flow  Tube  with  Diluent  Nitrogen. 
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For  a  perfect  gas,  Eq.  25  may  be  written 


PMcb  dT 
RT  3 T 


=  QW 


For  a  first  order  reaction,  Eq.  26  becomes 


PKlc],  dT 
RT  dt 


k,CiQ 


and  for  a  second  order  reaction,  Eq.  26  becomes 


PMcJj  dT 
RT  dt 


k2ctc2Q 


(26) 


(27) 


(28) 


Since  consumption  of  the  reactants  is  very  small  as  discussed 
in  the  preceding  section,  the  concentrations  of  reactants  can  be  considered 
constant.  Then,  expressing  concentrations  in  terms  of  initial  pressure, 
temperature  and  mole  fractions,  Eq.  27  for  a  first  order  reaction  becomes: 

il  =  Qki*'  dt 

T  Sfcyfo  d‘  (29) 

where  xj  is  the  mole  fraction  cf  specie  1  and  P  is  the  initial  pressure 
which  remains  constant  in  the  run.  Equation  29  can  be  readily  integrated 
if  kj,  the  first  order  ra£e  constant,  can  be  considered  constant.  By 
definition,  k  =  Aexp  |-  p'jj-  The  pre-exponential  factor  A  is  considered 
a  constant.  *  For  the  high  temperature  runs  now  under  consideration,  T 
varies  by  no  more  than  5  C°  during  the  runs  (Figs.  28  and  29).  Also,  a 
low  activation  energy  E  is  indicated  by  the  fact  that  the  reaction  proceeds, 
and  rather  rapidly,  at  room  temperature.  Under  these  conditions  of  a 
small  temperature  change  and  a  low  activation  energy,  say  4.  5  kcal/mole 
which  is  close  to  the  activation  energy  (5.  2  kcal/mole)  of  gas  phase 
ignition  reactions  of  NO2/MMH0),  the  rate  constant  k  vanes  by  less  than 
10%. 


*  Strictly  speaking  A  is  only  constant  in  the  case  of  a  reaction  of  the 
first  order.  When  the  order  of  reaction  is  greater  than  one,  A  is 
a  weak  function  of  temperature. 
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Taking  k  as  constant  and  noting  that  T  =  TQ  when  t  =  0, 
Eq.  29  can  be  integrated,  giving: 


(30) 


Expanding  ln(T/T0)and  neglecting  terms  other  than  the  first  since  T/TQ 
is  close  to  unity,  Eq.  30  becomes: 

_  _  Qk.x . 

T-T0  =  ATr..ction  =  *  (31) 

It  is  seen  therefore  that,  in  the  case  of  a  first  order  reaction,  the 
temperature  rise  as  a  function  of  time  is  independent  of  pressure.  Since 
the  high  temperature  runs  of  Figs.  28  and  29  indicate  a  significant  effect 
of  pressure,  it  is  concluded  that  the  overall  order  of  reaction  for  forma¬ 
tion  of  vapor -phase  reaction  intermediate  is  not  one. 

Again  expressing  concentration's  in  terms  of  initial  pressure, 
temperature  „-!d  mole  fractions,  Eq.  28  for  a  second  order  reaction 
becomes: 


dT 

T 


Uk|Xtx2 


Pdt 


(32) 


Equation  32  can  also  be  integrated  under  the  above  conditions  giving, 
after  expansion: 


T"T°  5  *Tre.ct.on  =  P'  133) 

in  this  case,  for  a  second  order  reaction,  the  temperature  rise  depends 
directly  on  pressure.  Also,  the  second  order  rate  constant  can  in 
principle  be  calculated  from  the  slope  of  the  A  T-t  curver  if  the  heat  of 
reaction  Q  is  known  for  the  vapor-phase  formation  of  the  reaction 
intermediate. 

Analogously,  it  can  shown  that  for  a  third  order  reaction, 
tne  temperature  rise  undei  the  above  conditions  depends  on  the  square 
of  pressure. 

As  a  consequence  of  the  above,  the  present  experimental  system 
provides  a  rather  sensitive  method  for  determining  the  overall  order  of 
reaction  if  condensation  effects  are  absent  The  latter  ondition  was  foui.d 
to  be  met  only  in  the  runs  at  elevated  initial  temperature.  From  these 
runs,  it  is  found  that,  since  the  temperature  rise  increases  directly  with 
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pressure,  the  reactions  for  formation  of  gas -phase  reaction  intermediate 
have  an  overall  order  of  approximately  two.  As  stated  above,  in  prin¬ 
ciple  the  second  order  rate  constant  can  be  calculated.  He’  in  the 

present  case,  the  required  value  of  Q  for  the  gas  phase  t  ./*  m,  with¬ 
out  condensation,  is  unknown. 

Activation  energy  is  conventionally  determined  fi  om  measure¬ 
ments  at  various  initial  temperatures.  The  same  reactions  muiv  of 
course  occur  at  the  various  temperature  levels.  In  the  present  case, 
an  additional  process,  namely  condensation,  occurred  1~  the  room 
temperature  runs  only.  Since  the  condensation  effects  cannot  be 
adequately  separated  in  the  present  limited  data,  activation  energy 
cannot  be  extracted. 

G.  EFFECTS  OF  FORMATION  OF  REACTION  INTERMEDIATE  ON 
PRE-IGNITION  CHAMBER  PRESSURIZATION 

Formation  of  the  reaction  intermediate  during  an  ignition 
delay  period  affects  ignition  delay  time  in  two  principal  ways.  Forma¬ 
tion  of  the  intermediate  liberates  heat  which  tends  to  increase  the  rate 
of  pre-ignition  pressure  rise  and  to  shorten  the  ignition  delay  time. 

Also,  formation  of  the  intermediate  consumes  reactants.  If  the 
reactants  are  vapors,  which  is  taken  to  be  t.ie  case  here,  then  pre- 
lgnition  chamber  nressure  tends  to  rise  more  slowly  and  the  ignition 
delay  time  tends  to  increase.  Thus  the  two  effects  are  in  opposition, 
the  heat  evolved  tendmg  to  increase  pre-ignition  pressure  rise  and 
the  vapor  mass  consumed  tending  to  slow  pre-ignition  pressure  rise 
Figures  1  l  and  12  indicate  that  the  dominant  effect  is  that  of  heat 
evolution.  However,  both  effects  are  analytically  treated. 

1 .  Mass  Sink 

Because  of  its  low  vapor  pressure  (less  than  I  rrmHg}  at  the 
temperatures  encountered  during  an  ignition  delay  period  (Sec  III  D). 
the  reaction  intermediate  is  treated  as  existing  wholly  in  the  condensed 
phase.  That  is.  the  reaction  intermediate  is  treated  as  forming  in  the 
vapor  phase  from  vapor  phase  reactants  but  then  immediately  condensing. 
The  controlling  process,  then,  is  formation  of  vapor  phase  reaction 
intermediate.  Surface  catalyzed  reaction*  on  the  droplets  formed  are 
assumed  to  be  negligible. 
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For  a  second  order  reaction  (Sec.  IV.  F.  3)  with  partial 
orders  with  respect  to  fuel  and  oxidizer  of  unity,  the  rate  equation 

it 


a  W  *  Acfcoxexp  (-i) 


Rt 


4) 


For  a  small  but  firite  time  interval,  A  t*r,  Eq.  34  becomes 


-  Acf  =  Ac fcox  A  tN  exp  ( 


(35) 


Expressing  concentrations  in  terms  of  mass  of  vapor  phase  reactants 
in  a  thrust  chamber,  rne,  for  compatibility  with  the  equaaons  of  Sec. 
Ill,  and  relating  the  amount  of  fuel  consumed  to  the  amount  of  reaction 
intermediate  formed  utilising  the  stoichiometry  given  in  Sec.  IV.  E.  3, 
Eq.  35  becomes: 


m 

nt 


1  0  ;f  '  M  '  ox 


AtN  exp  (  -  ) 


(36) 


Equation  36  is  used  to  calculate  the  mass  of  reaction  intermediate 
formed  in  each  time  interval  A  during  the  ignition  delay  time  calcu 
lation. 


Sir.ce  the  reaction  intermediate  is  formed  from  vapor  phase 
reactants,  the  mass  of  vapors  ol  each  reactant  are  correspondingly 
reduced  in  ea^h  time  interval.  6.  in  Sec.  111.  A  ,  is  there¬ 

fore  replaced  by 

VCM  ,  Pg.  - 

'mcfN  *  (f|)  *  9  imlilt)N 


(37) 


where  9-0.  571  and  0-  429  for  fuel  and  oxidizer,  respectively,  as  given 
by  the  stoichiometry  of  the  formation  reaction  derived  in  Sec.  IV.  E  3 


For  the  temperature  T  m  Eq.  36.  the  average  temperature  of 
the  reactant  vapors  in  each  t  ime  interval  is  used: 


(mcc 


Sav 


2&:s, 


If  *  (rrccpgTg)ox 


(mcCpg){  *  ?mccpg^ox 


{  38) 


T  for  each  reactant  is  obtained  from  Eq.  7  but  modified,  as  shown  below. 

g 

to  reflect  warming  of  the  vapors  by  the  heat  evolved  from  formation  of 
the  reaction  intermediate. 
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2.  Heat  Source 

Formation  of  the  reaction  intermediate  evolve?  heat  which 
raises  the  temperature  and  pressure  of  the  reactants  m  a  thrust 
chamber  and  thereby  tends  to  decrease  the  .ignition  delay  time.  Since 
the  reaction  intermediate  is  considered  to  exist  in  the  condensed  phase 
due  to  its  low  vapor  pressure,  the  heat  e.olved  is  both  tht  heat  of 
reaction  fc~  formation  of  the  gas  pha^e  react' on  intermediate  and  the 
heat  of  condensation  combined.  The  overall  reaction  is  taken  to  be 
the  one  given  in  Sec.  IV.  E.  namely: 

3NOz(g)  +  iCH3N2H3{g)  C4H24W906(1)  +  N^g) 

The  heat  of  this  reaction  was  calculated  to  be  A  Hr  =  -66.  2  kcal 
per  mole  of  MMH,  using  the  experimentally  determined  value  of  the 
hea*  of  formation  oi  the  liquid  reaction  intermediate  (Sec.  IV.  E.  2). 

Neglecting  sensible  heats,  the  heat  evolved  in  each  time 
interval  due  to  formation  of  liquid  reaction  intermediate,  Qr,  is 
therefore: 


Qr  =-(m;nt)N  AHr 


(39) 


where  AHr  is  the  heat  of  rr action  per  unit  mass  of  reaction  inter¬ 
mediate,  -0.900  keal/gr,  _md  (m.nt)N  is  given  by  Eq.  36. 


The  heat  evolved  in  each  time  interval,  Qr,  is  apportioned 
to  each  reactant  in  relation  to  their  vapor  phase  masses  in  the  thrust 
chamber  during  the  interval.  The  energy  which  raises  the  temperature 
of  the  oxidizer  vapors,  for  example,  is: 


Qr 

r  ox 


=  ( 


m. 


ox 


tnf  +  mnf 
Lox  ci 


)  Qr 


(40) 


Equation  7  for  the  temperature  of  the  vapors  of  one  reactant  at  the  end 
of  the  time  interval  is  modified,  giving  (for  the  oxidizer): 


-  r  ^(Tdi^N.j  +(Tg)N.lf(n,clN.!-(Gw)N-(mnoz)  1+ 

U1  cPg 

N  3 

.2,  Sj  Gij  +  (mc)N_i  '  (°W)N  '  (mno*>N 


The  effect  on  ignition  delay  lime  of  formation  of  reaction 
intermediate  is  given  in  Sec.  V.  which  follows. 


(41) 
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SECTION  V 

PHASE  III  -  REFINEMENT  AND  VERIFICATION 
OF  MATHEMATICAL  IGNITION  DELAY  MODEL 


In  this  section  of  the  report,  the  total  model  is  descriptively 
summarized  and  the  general  method  for  obtaining  solutions  is  outlined. 
Then,  computed  results  are  compared  to  experimental  results.  Finally, 
some  considerations  concerning  the  role  of  the  reaction  intermediate 
in  ignition  pressure  spikes  are  given. 

A.  INCORPORATED  REFINEMENTS 

In  the  present  program,  the  hypergolic  ignition  delay  model 
developed  under  the  preceding  program  has  been  refined  in  three  areas. 
Firstly,  the  initial  build-up  of  propellant  flowrate  in  a  thrust  chamber 
following  vulve  opening  has  been  incorporated  into  the  model  in  the 
manner  given  in  Section  III.  C.  Secondly,  heat  transfer  between 
thrust  chamber  walls  and  the  vaporizing  vapor/drop  system  has  been 
incorporated  as  given  in  Section  HI.  D.  These  first  two  refinements 
affect  pre-ignition  pressure  rise  due  to  vaporization  of  each  propellant 
independent  of  the  other.  The  third  area  concerns  reactions  between 
the  two  propellants  to  form  a  reaction  intermediate  which  also  affects 
the  pre-ignition  pressure  rise  in  a  thrust  chamber.  Formation  of  the 
reaction  intermediate  during  an  ignition  delay  period  is  accounted  for 
hi  the  manner  given  in  Section  IV.  G. 

B.  GENERAL  METHOD  OF  SOLUTION 

The  analytical  pressure-time  curves  given  in  Sec.  IH  were 
calculated  for  the  single  propellant  under  consideration  using  the 
propellant  vaporization  model  of  Sec.  III.  A.  ,  modified  as  described 
in  Sec.  III.  G  and  D,  For  calculating  ignition  delays,  the  vaporiza¬ 
tion  model  is  used  twice,  once  for  each  propellant,  and  the  calcu¬ 
lations  of  each  are  performed  simultaneously.  That  is,  the  calcu¬ 
lations  are  performed  for  one  time  interval  for  MMH  (Atjyj  =  50  p 
sec)  and  then  for  two  time  intervals  for  N2O4  (  Atj^  =  25  M-  sec). 

This  pattern  is  repeated  throughout  the  solution.  The  shorter  time 
intervals  in  the  case  of  N2O4  are  required  due  to  the  oxidizer's  h.gh 
volatility. 


At  the  er.J  of  each  50  microsecond  time  interval,  the  amount 
of  reaction  intermediate  formed  and  the  heat  evolved  during  that  interval 
are  computed.  As  discussed  in  Sec,  IV.  G.  .  the  evolved  heat  is  applied 
to  the  reactant  vapor  masses  in  the  chamber.  The  mass  of  reactants 
which  are  consumed  in  the  formation  of  the  reaction  intermediate  are 
also  determined  and  deducted  from  the  vapor  masses  remaining  in  the 
thrust  chamber. 
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Also  calculated  at  the  end  of  each  50  microsecond  interval 
is  the  fractional  degree  of  ignition  which  occurred  during  that  Interval. 
The  fractional  degree  of  ignition  is  simply  the  length  of  the  time 
interval  (50  microseconds)  divided  by  the  ignition  delay  time  (Tjg)N 
calculated  for  the  conditions  of  pressure,  temperature  and  composition 
of  the  reactant  vapors  prevailing  at  that  time.  The  fractional  degrees 
of  ignition  are  summed  over  succeeding  time  intervals  until  the 
fraction  becomes  unity,  indicating  ignition. 


The  analytical  expression  used  to  calculate  ( Tig)N  ^or 
each  time  interval  is: 

<^N  =  iri «*(i3sr) 

where  P  =  pf  +  pox,  3  =  pox/pf  and  Cpf  and  Cpox  are  molar  heat 
capacities.  T„  is  obtained  by  Eq.  38. 


Equation  42  is  derived  in  (1)  from  the  heat  equation  for  an 
adiabatic  volume  element  of  vapor  phase  reactants  undergoing  thermal 
bimolecular  ignition  reactions.  The  derivation  involves  the  Todes 
approximation  (12)  for  solution  of  the  differential  equation.  Although 
the  approximation  applies  to  values  of  Ejgn/RT  >17,  we  use  the  re¬ 
sulting  equation  42  in  the  present  case,  for  which  E-gn/RT  ~  9, 
in  the  absence  of  a  better  criterion  of  ignition. 


The  values  of  E-gn  and  (AQ)jgn  required  by  Eq.  42  were 
determined  experimentally  during  the  previous  program  (1).  The 
values  appropriate  to  the  present  propellant  combination  are  E^n= 

5.  2  kcal/mole  and  (AQ)^gP  =  3.  4  x  10*^  cal-ce/mole  -sec. 

The  fractional  degree  of  ignition  fer  each  time  interval  is 
simply  Atjq/(T-g)j^  where  (t  j g)isj  is  calculated  by  Eq.  42  at  the  end 
of  each  50  microsecond  interval.  The  ignition  fractions  are  summed 
over  successive  intervals.  The  ignition  delay  time  in  an  engine,  then, 
is  that  time  *r»  the  computer  solution  at  which  Tig)N|  becomes 

unity.  **  ■* 


C.  COMPUTED  RESULTS  AND  COMPARISONS  WITH  EXPERIMENTAL 
RESULTS 


1.  Ignition  Delay  Times 

Ignition  delay  times  were  calculated  for  N^O^/MMH  in 
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space-ambient  engines  in  the  manner  outlined  above.  The  actual 
computer  program  listing  is  reproduced  in  Appendix  A. 

The  value  of  activation  energy  for  formation  of  reaction 
intermediate  used  in  the  calculations  is  4.  5  kcal/mole.  Although 
the  activation  energy  could  not  be  determined  from  the  flow  tube 
temperature  rise  measurements  because  of  condensation  effects 
(Sec.  IV.  F.  3),  the  activation  energy  is  known  to  be  low  since  the 
reaction  intermediate  forms  at  room  temperature.  Furthermore, 
since  reaction  intermediate  forms  at  conditions  where  ignition  does 
not  occur  (in  the  flow  tube  experiments,  for  instance),  the  activa¬ 
tion  energy  of  the  former  must  be  lower  than  that  of  the  ignition 
reactions  (5.  2  kcal/mole),  A  value  of  4.  5  kcal/mole  was  there¬ 
fore  chosen  for  E^nt,  the  activation  energy  for  formation  of  the 
reaction  intermediate. 

With  Ejnt  =  4.  5  kcal/mole,  a  series  of  calculations  was 
run  in  which  Aint,  the  p re -exponential  factor  of  the  Arrhenius 
equation,  was  varied  in  order  to  obtain  for  Engine  Configuration 
No.  1  the  best  fit  between  (1)  calculated  and  experimental  pre- 
ignition  chamber  pressure-time  histories  and  (2)  ignition  delay 
times.  As  shown  in  Fig.  34,  a  good  fit  was  obtained  for  Ajnt  = 

1. 7  x  10^  *  cal-cc/mole  ^-sec.  This  value  was  then  used  in  a 
computer  run  for  Engine  Configuration  No.  2  and,  again,  as  shown 
in  Fig.  35,  good  agreement  resulted.  The  calculated  and  experi¬ 
mental  ignition  delay  times  (from  valve  signal  to  ignition)  for  the 
two  engines  are: 


Engine  Configuration  No,  1 
Engine  Configuration  No.  i 


Ignition  Delay  Time 

Calculated 

Measured 

msec 

msec 

5.  9 

5.  9 

6.  55 

6  6 

The  computer  input  data  used  in  the  calculations  for  Engine 
Configuration  No.  1  are  given  in  Appendix  B.  The  conditions  for  the 
corresponding  experimental  test  (Run  No.  37)  are  given  in  Table  IV. 

Also  plotted  in  Figs,  34  and  35  are  chamber  pressurization 
curves  calculated  without  permicting  formation  of  reaction  intermediate 
during  the  ignition  delay  periods,  i.  e.  with  A- nf  -  Eint  -  0.  The 
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Experimental 
(Run  No.  37) 
Calculated  with 
Aint  =  1.7  x  10n 
cc/mole-sec 
and  Eint  =  4.  5 
kcal/mole 
Calculated  with 

-^int  =  ^int  =  0 
(Reaction 
Intermediate 
not  permitted 
to  form) 


^Ignition 

(Analytical) 


“Ignition 
(Exp. ) 


Valve 

Signal 


4  6 

Time,  msec 


Figure  34.  Calculated  and  Experimental  Pre-ignit.on  Chamber  Pressui 
Histories  for  N^/MMH  in  Engine  Configuration  No.  1 
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60 

33  200 
6 
0 


Experimental 
(Run  No.  42) 
Calculated  with 
Aint  =  1.7  x  10I! 
cc/mole-  sec 
and  Eint  =  4.5 
kcal/ mole 
Calculated  with 

Aint  =  Eint  =  0 
(Reaction 

Intermediate 
not  permitted 
to  form) 


—Ignition 
^  (Exp. ) 
[♦Ignition  - 
(Analytical) 


, 

Ilf 


Valve 

Signal 


0  2  4  6  8 

Time,  msec 

Figure  35.  Calculated  and  Experimental  Pre  -Ignition  Chamber  Pressure 
Histories  for  in  Engine  Configuration  No.  2 
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difference  between  these  curves  and  the  other  calculated  curves  of 

the  figures  indicates  the  significance  of  formation  of  the  reaction 

intermediate  on  ignition  delay  time.  Although  formation  of  the 

reaction  intermediate  removes  mass  from  the  vapor  phase,  the 

heat  which  is  evolved  more  than  compensates  for  the  vapor  mass 

removal  and  finally  results  in  markedly  shorter  ignition  delay’ 

times.  In  both  cases  in  which  A-  f  -  E.  =0,  the  calculated 

ini  lnt 

ignition  delay  times  were  greater  than  10  msec. 

The  heat  evolved  when  the  reaction  intermediate  forms 
raises  the  vapor-phase  reactant  temperatures  substantially,  as 
shown  in  Appendix  C,  Printout  of  Computed  Results  for  Engine 
Configuration  No.  1*  In  the  case  given  in  Appendix  C,  the  average 
gas  temperature  TGAV  reaches  800°F  as  computed  by  the  present 
model.  Even  at  temperatures  well  below  this  level  (i.  e.  above 
approximately  150°F),  condensation  of  the  reaction  intermediate 
does  not  occur.  Consequently,  it  is  improper  to  use  the  heat  of 
reaction  for  formation  of  liquid  phase  reaction  intermediate 
(  AHr  =  -66.  2  kcal/mole  MMH,  Sec.  IV.  G.  2.  )  throughout  the 
entire  ignition  delay  period.  A  further  consequence  of  the  indicated 
high  temperatures  during  the  latter  part  of  the  ignition  delay  time 
i3  that  the  reaction  intermediate  which  formed  in  the  earlier  part  of 
the  ignition  delay  period  will  decompose,  as  indicated  by  the  heating 
curves  of  Fig.  22,  or  possibly  undergo  reaction  with  the  warmed 
oxidizer  present  in  the  chamber.  In  summary,  the  good  agreement 
shown  in  Figs.  34  and  35  between  calculated  and  experimental  ignition 
delays  and  also  pre-igrition  chamber  pressure  histories  is  misleading. 
Our  understanding  of  events  occurring  throughout  the  entire  ignition 
delay  period  is  not  complete  and  further  study  of  the  latter  part  of  the 
ignition  delay  time  is  needed. 

2.  Fuel  Vapor  Condensation  onto  Oxidizer  Drops 

Ir  accordance  with  Eq.  4  (in  Sec.  III.  A.  ),  propellant  vapor 
will  condense  onto  a  surface  if  the  partial  pressure  of  the  propellant 
vapor  exceeds  its  vapor  pressure  corresponding  to  the  temperature  of 
the  surface.  The  computer  program  predicts  that  many  N2O4  drops 
evaporatively  coo)  to  a.  level  such  that  MMH  vapor  will  tend  to  condense 
on  these  cold  N2O4  drops.  Expressed  in  another  way,  the  partial 
pressure  of  MMH  vapor  exceeds  the  vapor  pressure  of  MMH  corres¬ 
ponding  to  the  temperature  of  many  of  the  N2O4  drops.  Thus,  MMH 
vapor  will  tend  to  cordense  onto  the  colder  oxidizer  drops.  The  reverse 
situation  was  not  found  to  apply.  1.  e.  the  fuel  drops  do  not  cool  enough 
to  cause  condensation  thereon  of  oxidizer  vapors. 
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The  computer  program  as  it  presently  exists  does  not 
allow  for  fuel  vapor  condensation  onto  oxidizer  drops  but  the  com¬ 
puted  results  do  predict  that  this  phenomenon  will  tend  to  occur. 

A  question  is  therefore  raised  as  to  what  chemical  reaction  will 
occur  if  fuel  vapor  does  penetrate  the  evolving  oxidizer  vapor 
surrounding  the  cold  drops.  The  reaction  can  lead  to  forma¬ 

tion  of  additional  reaction  intermediate  or  to  complete  oxidation  of 
the  condensed  fuel.  Studies  of  the  reactions  of  a  fuel  droplet  (N?H4) 
in  an  oxidizer  atmosphere  {NO2)  are  presently  underway  (13);  how¬ 
ever,  the  reverse  situation  of  a  (cold)  oxidizer  drop  in  a  fuel  vapor 
atmosphere  has  not  been  investigated  to  the  authors'  knowledge. 

D.  IGNITION  PRESSURE  SPIKING  CONSIDERATIONS 

Some  qualitative  insight  into  the  possible  role  of  the  reaction 
intermediate  in  ignition  pressure  spiking  can  be  gleaned  from  the 
sensitivity  tests  discussed  in  Sec.  IV.  B.  To  obtain  further  insight, 
the  energy  release  upon  sudden  decomposition  of  the  reaction  inter¬ 
mediate  has  been  calculated  as  well  as  the  energy  release  upon  re¬ 
action  with  additional  oxidizer. 

Using  the  empirical  formula  deduced  in  Sec.  IV.  F.  3 
for  the  reaction  intermediate,  an  equation  for  self-decomposition 
to  H2O,  CO2,  N2  and  excess  MMH  may  be  written: 

5  C4H24N906  (1) - >  6  C02(g)  +  18  H20(g)  +81/2  N2(g)  +  14  CH3N2H3(g) 

Using  a  heat  of  formation  of  the  reaction  intermediate  of  AHf  =  -0.  509 
keal/gr  (Sec.  IV.  F.  2),  the  heat  of  reaction  for  the  above  decomposition 
is  computed  to  b:4Hr  =  -0.  365  keal/gr  cf  reaction  intermediete, 
indicating  a  relatively  unenergetic  reaction. 

For  reaction  with  additional  oxidizer,  the  following  equation 
can  be  written: 

C4H2<N9Ofe  (1)  +  7N02(g)  - »  4C02(g)  +  12H20(p1  +  8  N  Jg) 


The  heat  of  this  reaction  is  computed  to  be  AHr  =  -3.  32  keal/gr  of 
reaction  intermediate  which  is  similar  to  the  heat  of  combustion  of  TNT 
with  oxygen  (10).  Thus,  the  reaction  of  the  reaction  intermediate  vith 
additional  oxidizer  possesses  considerable  energy  but  the  rate  at  which 
the  energy  is  rel"ased  is  not  known  at  present.  This  parameter  is  of 
course  all  important  for  ignition  pressure  spiking  and  for  TOT's 
(oxidizer  manifold  explosions)  too. 
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The  total  quantity  of  reaction  intermediate  formed  during 
the  ignition  delay  period  of  Engine  Configuration  No.  1,  according 
to  the  present  computer  model,  is  given  in  Appendix  C  under  the 


>-5 


lb- 


does 


column  labelled  TM1NT,  LB.  If  this  quantity,  2.  33  x  10'  wm, 
react  with  additional  oxidizer  in  accordance  with  the  preceding 
equation,  2.  55  x  10°  lbm  of  oxidizer  are  required.  As  given  in 
Appendix  C,  the  amount  of  oxidizer  in  the  vapor  state  in  the  chamber 
at  ignition  is  2.  84  x  10”^lbm,  which  exceeds  slightly  the  required 
amount  for  complete  reaction  of  all  of  the  reaction  intermediate 
formed.  Because  the  reaction  possesses  sufficient  energy  and 
because  sufficient  oxidizer  is  present  to  permit  the  reaction  to 
proceed,  it  is  strongly  suspected  that  the  reaction  intermediate 
plays  a  critical  role  in  severe  ignition  pressure  spikes.  Further 
work  is  required  in  this  area. 
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SECTION  VI 
CONCLUSIONS 

A  workable,  fundamental  model  of  hypergolic  ignition  in 
■pace  ambient  engines  has  been  developed  and  verified.  The  com¬ 
puter  model  satisfactorily  predicted  the  pre-ignition  chamber  pressure 
histories  and  ignition  delay  times  of  two  different  engine  configurators. 

In  developing  the  model,  considerable  insight  has  been 
gained  into  the  physical  and  chemical  processes  operative  during  the 
ignition  delay  time.  Both  theoretical  and  experimental  evidence  has 
established  that  a  reaction  intermediate  of  N2O4/MMH  forms  during 
the  ignition  delay  period  and  markedly  affects  (shortens)  the  delay 
time. 


The  composition  of  the  N2O4/MMH  reaction  intermediate 
is  very  complex.  It  appears  to  be  a  mixture  consisting  of  at  least 
four  compounds,  the  principle  ones  of  which  are  MMH'  H2O  and 
MMH-  HNO3.  The  material  is  more  sensitive,  both  thermally  and 
to  impact,  than  MMH-  HNO3  prepared  in  the  laboratory. 

Based  on  its  fuel  contei  ,  the  liquid  reaction  intermediate 
still  contains  approximately  79%  of  the  energy  of  complete  combus¬ 
tion.  The  reaction  of  the  material  with  additional  oxidizer  which  i3 
present  in  a  thrust  chamber  m  sufficient  amounts  at  ignition  is  very 
energetic,  being  similar  to  that  of  TNT  burning  with  oxygen.  By 
comparison,  self-decomposition  of  the  material  is  relatively  un- 
energetic. 


Due  to  gradual  decomDosition  and  its  hygroscopic  tendency, 
careful  sampling  techniques  and  prompt  chemical  analyses  are  required 
to  obtain  meaningful  composition  'ata. 

The  computer  model  predicts  that  during  the  ignition  delay- 
period  MMH  vapor  tends  to  condense  onto  the  colder  N2O4  drops. 

The  kind  of  reactions  which  occur  (formation  of  reaction  intermediate 
or  complete  oxidation)  if  the  fuel  vapor  does  penetrate  the  evolving 
vapor  from  the  cold  oxidizer  drops  has  not  been  determined.  The 
model  predicts  that  the  reverse  situation,  condensation  of  oxidizer 
vapor  onto  cold  fuel  drops,  cioes  not  occur. 
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SECTION  VII 
RECOMMENDATIONS 


Several  factors  pertinent  to  hypergolic  ignition  of  N2O4/ 
hydrazine -type  fuels  in  space-ambient  engines  that  were  either 
outside  the  scope  of  the  present  program  or  require  further  study 
beyond  that  which  was  permitted  during  this  six-months  effort 
are  enumerated  below: 

1.  Specie  analyses  of  the  N2O4/MMH  reaction  intermediate 
identified  approximately  85%  of  the  material.  Further  work  is 
required  not  only  to  identify  the  remaining  15%  but  also  to  eliminate 
anomalies  that  exist  between  the  present  deduced  composition  of  the 
material  and  its  measured  thermal  stability,  impact  sensitivity  and 
vapor  pressure. 

2.  Because  the  reaction  intermediate  is  a  mixture  of  compounds, 
the  effects  of  initial  reactant  mixture  ratio,  pressure  and  temperature 
on  its  composition,  properties  and  sensitivity  should  be  investigated. 

3.  Further  work  is  required  to  determine  formation  and  con¬ 
densation  rates  of  the  reaction  intermediate  over  the  temperature 
range  encountered  during  ignition  delay  periods  in  space -ambient 
engines.  It  should  be  noted  that  the  composition  of  the  reaction 
intermediate  undoubtedly  varies  from  the  start  of  an  ignition  delay 
period  to  ignition. 

4.  To  further  define  the  role  of  the  reaction  intermediate  in 
ignition  pressure  spiking,  reaction  rates  of  the  reaction  intermediate 
with  additional  oxidizer  need  to  be  determined. 

5.  To  obtain  a  basis  from  which  to  develop  preventive/corrective 
measure:  to  eliminate  ZOT's  (oxidizer  manifold  explosions),  it  is 
necessary  to  determine  the  mechanism  by  which  fuel  enters  the 
oxidizer  manifold  and  to  determine  the  state  of  the  fuel  (unreacted  or 
partially  reacted  in  the  form  of  reaction  intermediate)  just  prior  to 
entry  of  the  oxidizer. 

6.  Considerable  progress  has  been  made  in  identifying  the 
N2O4/MMH  reaction  intermediate  and  in  determining  its  physical 
and  chemical  properties.  Similar  characterizations  of  the  reaction 
intermediate  formed  by  N2O4  with  UDMH,  N2H4  and  50-50  need  also 
be  determined. 
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7.  The  effects  of  NO-inhibited  on  the  formation  and 

properties  of  the  reaction  intermediates  remain  to  be  determined. 

It  has  already  been  determined  that  as  little  as  1%  NO  in  at 

least  doubles  the  ignition  delay  time  of  ^O^/UDMH  in  an  uncon/ined 
impinging  stream  apparatus  (Z). 

8.  A  mathematical  model  of  ignition  pressure  spiking  during 
start  transients  should  be  developed  to  permit  evaluation  of  chamber 
geometry  and  ignition  hydrauiics  from  a  pressure  spiking  viewpoint. 

The  present  hypergolic  ignition  delay  model  becomes  a  major  part 

of  such  a  pressure  spiking  model  since  the  ignition  delay  model  gives 
the  conditions  in  a  thrust  chamber  at  ignition,  that  is,  the  conditions 
from  which  ignition  pressure  spikes  result. 

9-  The  present  hypergolic  ignition  delay  model  should  be 
exploited  to  develop  design  information  by  first  isolating  the  main 
"arm  .iters  and  then  determining  the  rate  of  change  of  ignition  delay 
:e,  mass  of  intermediate  formed,  etc.  Parameters  such  as 
propellant  lead  and  lag,  transient  flow  increase  from  time  zero 
;feed  system /injector  effect),  initial  propellant  temperature,  initial 
hardware  temperature,  etc.  can  be  investigated  analytically  using 
the  present  model  to  economically  assist  engine  development  programs. 
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APPENDIX  A.  COMPUTER  PROGRAM  LISTING 


C  HYPEHGOUC  IGNITION  SIMULATION 

01  ME  NS  I  ON  R<  3,2111  ,TO(3,2<.l)  .XI  3. 20  1 1 
DOUBLE  PRECISION  OELTN, 

1  TOC ,WD0T,VC , AST  AS , AC, V,LC,PVW,M, ALPHA, RHOL.LMBDAE  , LMBOAS, LMBDAF , 

2  CPL ,CPG,TDFP,K,ME WG.KG.Ol , 02 , 03 , 04 , 05 , 06 , D 7 ,DB , 09 , 01C . 01 1 . 

3  fH2,T!ME,PG,TG,TEMP,GN,HN,UN,C< 3 ) • P V , GAM , Q , 7  *  OMEGA,  SIGMA. 

4  CHI  ,F*  NC.GW,NNOZ,TGPRM,PGPRH,SQPTRE,MINT 
DIMENSION  OKI  3 , 4f’  1  I  *0T0(  3, AO  1  (.OKI  3, 40  II 
DOUBLE  PRECISION  OOELTN, 

1  OTG',.OWOOT*OVC*OASTAR,nAC*GV«nLC»OPV'W*C’N*OALPHA,  ORHOL  .OLNBOE* 

2  OLMBOS.OLMODf , UCPL , OCPG , OTOF P ,0K ,OME  wG ,  QKO ,  DO  1 , 002  ,  t)D3 , 0D4 ,  DOS , 

3  OD 6. 007* 008* 009*001 0*0011  , 00 1 2 * OPG , OTG , OT E NP, OGN, OHN, OON ,OPV. 

4  OGAM,00,OZ,OOMEGA,OSIGNA,OCH1 ,0E , UMC , OGW , OMNOZ , OTGPRM , OPGPRM, 

5  OSCRT 

100  FORMAT!  nES.4) 

IC1  FORMAT! IDEA. 4/2F16. 41 

K 2  F  UK  M  A  f ( I  5, OPE  If .3,f 12.3, F 14.3, 2X, 1P4E15.5, JPF8.2.F 10.2.F9.4/) 

If  3  FORMAT! I  15, 1P3E IS. 5,31 OPF 13. 3,2X1, 1PE 15.51 
1C4  FORMAT! //IX, 'TIG**  ,F  7,4,  l*.»MWHt ) 

1C5  FORMAT! 1H1 ,6X,' INPUT  DATA'//I 

106  F0PMAU6X,*  INITIAL  CONDITIONS  -  FUEL'//) 

1106  FORMAT! /6X, • INITIAL  CONDITIONS  -  OX*//! 

1C  7  FORMAT! 6X, 'OELTA-T • ,9X, « TOO • • 9X  * • R ! I  *  1 >',M,'RI2,l)*,8X,»Rl3,ll', 

1  PX,  'TO!  l,H'#X,,TL'!2,l»,.7X,*TO<3,n«//l 
IC8  FORMAT! 1P9E14.5I 

109  FORMAT! //6X, 'CHAMBER  OATA'I 

1 10  FORMAT!  7/7 X,'  MOOT'  , 1 1 X , • VC ' , 1 2 X, • A* • , 1 2 X , • AC  * , 1 3 X, • V* • l 3X  ,  •  LC  ' , 
l  12X,  '  P  VW*  ,BX,  'SSFLT'  »9X»  • FRCTH'  * 

111  F0RMATI//6X, 'PROPELLANT  OATA'I 

112  FORMAT! //9X, 'H' ,1 IX, 'ALPHA* ,9X, 'RHOL ’ ,8X, • L AMBOA-E • , 6X , • L AMBDA- S* , 
1  6X,  'LAMBDA-F • ,9X , 'CPL * , 1 1 X , *C PG ' , 1 1 X, • TOP P*  I 

113  FOP  MAT!  /  /RX,  'X'  ,12X,  'll-f.  •  •  1 1  X  ,  'K-G*  ,  1 1  X  ,  '  0 1  •  ,  1 2X  # 
l  •02S12X,'03S12X,,04',12X,  '05'  I 

114  FORMAT!  ///9X, 'N*  ,4  X  ,  •  T  I  M£  ,  MS  •  ,  S<  .  •  PG  ,  m«h  *,  •  .  5  X  ,  •  T  G,  OEG  .  R  •  ,  8  X  , 

1  •MN0Z,LB',8X,'NC,LB,,ir<,'GN,Lft'//l 

115  FORMAT!  7//UX,' J*  ,5X,  'P  (  J,  1  ),FT',6X,'R!J,2I,FT',6X,»R!J,3I,FT',7X, 
l  'TO!J,ll',8X,'TO(J,2>',ex,'TO(J,3)',e»,'X(j,li'//l 

116  FORMAT! 1H1.3* ,'N' ,4<, *TI ME , MS • , 5 X , • PG , MMHG' , 5X  ,  • TG.OEG.R' ,8X, 

1  'MNOZ.LB' ,ax,»MC,LB' ,1CX, 'GN.LB'.BX, • TM|NT,LB',3X,'TIGN,NS* , 

2  5X,'TGAV,R',3X,'$Mir,FR*//l 

117  FORMAT! //7X, *06* ,12X, '07' ,l?Xf 'HA', l 2 X , • 09 • , 1 2  * , '010', 1 IX, 

1  'Oil', 11*. *012'. 11X. '013'  I 

118  FORMAT!  /  /7X,  '01W  ,9X.  'FRCU'  ,9X, • TOHT  *  I 

119  FORMAT! lE7.2ttX.lET. 21 

l  Z*  FORMAT! 1PE9. 2,1*, lP£9.?/< 

REA0(5,119»  AOIGN.AEINT 
WR1  T E  ( 6, 12') I  AQIGN.AEINT 

READ! 5, 100)  RCAL,THETA,OTHETA,OFLHR,EAIGN,EA|NT 
WRITE (6,  1081  RCAL, THETA, OfHETA.OFL HR ,EAIGN, FAINT 
C  READ  CHAMBER  OATA 

1  REA0I5, 10(  »  OELTN, TG>, MOOT, VC, AST -k, AC, V.LC.PVW, SSFLT, FRCTH, 
l  OIW, FRCTL. TOHT 
C  READ  PROPELLANT  OATA 
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READ! 5, I'M  I  M, ALPHA , RHOL »LMBOA£ , LMftD AS , LMBOAE.CPL »CPG» TDFP , K, 
l  MEWG,KG 

C  SET  INITIAL  CONDITIONS 

READ!  5,  IDO  »  mj,l»,J»l,3>,<TD<J,ll,J*l,l) 

01*1546. 0 
02*12.7 

CLOSURE (M*l)2/l6.28  31853*01 1 1 
04-ALPHA*AC*DELTN*03 

D6.OSORT(K*(2.L/(K»I.0n**(lK*l.T  I/IK-1.1111 
07*ASTA«*0FLTN*D6*0SyRTIM*0?/ni» 

OH  VC4M/D1 

09*2.0*V«M/(MEWG*Dll 
OK*ICPG*M£wG/KCI**<  1.G/3.CI 
Oll*LC/( V*OElTNI 
Ol2*2.o *RHOL*CPG/ I  3. }*KG*DELTNI 
D13-0IW/V/0ELTN 
f.  READ  CHAMBER  DATA  N?04 

IOC  l  READ!  5, 10(  »  OOELTN.IJTG'I.  JWIXIT ,OVC , OAS T A» , OAC . OV, OLC » OPVM, 

1  OSSELT.OFRCTH.ODIM.OPRC  tl.otoht 

C  READ  PROPELLANT  DATA  N?CK 

RE  AO  I  5, 101 J  OM,nALPHA.nRMf)L.nLMCH>c.OLJ<BOS.OL*BDE,Oef»l.!K'?G,OTOfP. 

2  OK i CMf WG.OKG 

C  SET  INITIAL  CONDITIONS  N204 

READ!  5,  lOfl  <U«»J.l»,J*l,3l,(0T0U.t>,J.l,}) 

003*  OSQR T(0M*02/I6.253I  853*0111 
004»  DAL  PHA*IJAC*tlOE  L  TN*Dt) } 

U06*DSQRTI')K*<  2.1  / I  OK ♦ 1 . 0 H ♦  * II DK ♦ 1 . 0  I / IDK -l ,C» II 
0O7-OA5TAE  AODf  LTN*ni)6«OS0PT(  ()M*0  7/0I  I 
Q08*0VC  *OM/oi 
009»2.»*r'V*0M/(UME  «0*0l  I 
ooit»cocpcaom£m&/okgi**i l .o/n.ri 
onii-nLC/invpfiDtiTM 

0012*2.  '•OBH0L*«lCP(./ 1  i.:  ■''*>. •f’l'f  I  f  n  J 

oni'.df*-  4  . 

N«l 

NN-1 

TIME*'.  . " 

Pr,rf  ,n 
TG-TG1' 

OpO*r  .r 

UTG«OT&o 
MINT  *'j  ,r 
TM|M|.f  ,o 
SN|(.f  R*".r 
,r 

nys  •  r.o 

DO  5  >  11*1,3 
DO  4'  I2-1.2T1 
OKI  I  1  » I  2l«f'.C 
0«l  I  l.  I2»?''0|.r  .  ) 

50  <1 tl ,I?I*C." 

C1SI-'  .3' 

Cl  2  • 

C 1  3  >  -<* .  3  ' 
rf»  I  T  E  1 6 . 1 ' 5 » 
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WRI TE  ( 6, 1C6 1 

WR ! Te ( ft , 107 » 

WRITE (6, lOeiOELTN, TGO  ,R  ( 1 ,  1  > ,  R  <  2 , 1 J  ,R  (  3,  l ) .  7D< 1,11,70(2,11,70(3,1) 
WRI TF( 6, 1091 

wri re  c  ft, lioi 

WRI IE  16,1081  WD0T,VC,AST4R,AC, v,  ic , pvw, ssfl t,frc th 
WR  I  ?  E  ( 6  ,  1 1 8 1 

WR I TE ( 6  , 108  I  DIW*FRCTl, TOHT 

WRIT£(6,111I 

WR  I  T£  (  6, 1 121 

WR'THo,lC'8»  H» ALPHA  ,  RHQL  *LN80Afc  , LM9DAS  *LN80AF*CPL*CRG»  TOFP 
WRITE (6,1131 

WRITE (6, 108)  K,M£WG,KG,D1,02,D3,04,D5 
WRITEI6, 1171 

WRITE (6, 1C  81  06, 07,08, 09, QIC, 011 ,012 ,013 
WRITEI6, 11061 
WRITE! 6, 1071 

WRITE (6,1081  OOELTN.OTGC  « OR (1,11* OR (2*11,  OR  (3, II, 0TO(1, II, 

1  OTO ( 2 , 1 1 ,0T0( 3,11 
WRITEI6, 1091 
WRI TE (6 , 1 10  I 

WRITE(6,l(,8»  OWDOT ,OVC ,OASTAR «OAC,OV , OLC ,0RVW,05SFLT,0FRCTH 
WRI TE ( 6, 1 18  I 

WRITE(6,1E>8I  OOIW.OFRCTL.OTOHT 
WRI TE ( 6: 1 1 1 1 
WRI TE ( 6 , 1 171 

WRI TE(6, 108.  OM,OALPHA,ORHOL,OIMBOF,OLMBOS,OLWBDF,OCPL,OCPG,OTOFR 
WR I TE ( 6, 1 1 3 1 

WRI TE( 6,K 8  I  OK, OWE WG » OK G, 01, 02, 003, 004,003 
WRITE (6,1171 

WRITE (6, 108 1  006  j  007 , 008 ,006, 00 1^,0011 ,0012,0013 
WRITE (6,1161 

2  GN«<  ,o 
HN«0 .0 
ON«C.O 
J«N*l 

310  IE  (J-Oll-211  3,3,55 
55  J»D1 1*21 

3  J*J-l 
I 

05«W00T*0ELTN 

IE  (W-(  SSFIT/0ELTN*01  IM  300,4,4 
3 CO  IE  <  J-N*S4FU  /  OELTNI  4,4,305 
315  05,n5*IN-J*0.5l*0€lTN/SSFlT 

4  I  - 1  ♦  1 
G4N«  h ( ( , J) 

TENP«TO(  I  ,  J I 

IF  ( TENP-TOEPI  60,65,65 
60  TEWP«TOF» 

65  PV-PFUNCTI TtrPI 

IF  (J-OUI  70,70,5 
70  |F  (  TENP-TDFP  )  75,75,8'’ 

73  IF  I K(  1 ,31-1.01  80,5,3 

80  Q.3.**(PV-PGI*ALPMA*OfLTN/(KMnL*GAH»*O3*r>S0RT( l.O/TfNp| 


AFRPL-TR- 67-129 


85  IF  (PG»  5*5,9C 
9 C  $QRTRE*0SQRr(09*GAM*PG/7GI 

Z«012*GAM**2*G/<2.Cf J,6*01C*SaRTREI 
IF  (Z-30.M  20 5* 205 »  2r,i 
2<"  Z «  30  .  * 

2X*>  0MEGA*CPG'MTG-r£MP|/(0FMP(  Z I -1.0 1 

sigma-»om£oa*o5*c(  n*o 

(JN«ON«- SIGMA 

5  R(  l,JU»“GAM*(l.0-Ql**(l.9/3.*l 
TOfI  ,  J>U«rEMP-lMB0AE/CPL*8 
IF  (XII, Jit  210,21C,297 
2*7  TOC  l,JMI«TD FP 

IF  (0.6O-XC 1 • J I >  225,225  »4C0 
210  IF  (Q>  215,225,215 
215  IF  (PGI  225,225,220 
220  fO( I , J*1 l-TEMP»0/CPL*(0HEGA-LM50AE» 

4CC  IF  (J-D13I  225,401,401 
401  If  (X(I,JI»  4*5,405,236 
4C5  IF  <*-(Oi3*TOHT/OELT#ll  410,415,415 
410  FRC  T»(  FRCTH-FRCTL  )•(!-(  N-0 13 1  *OFLTN/ Tf)HT  I  ♦•FRCTL 
GO  10  420 
415  FRCT.FPCTt 

420  T0(  1,J*1  >«T0(  !,JMl4FRCT*f  TGO-TOI  t,j4lH 
225  CMI-O.f 

!0  ( T0( I ,J*l l-TOFP I  230,230,235 
230  CH|-0/LM80AF*(LMB0AS-0MEGAI 

235  X(I,J^ll*Xll,JKCHI 

IF  (XU.Jtlll  236,237,237 

236  X(l,jFll*o.O 

237  IF  (*(I,J»ll-l.JI  239,239,238 

238  X( I « J+l  I*  1.9 

239  F«C (  1 1 *05*0 
GN-GN+F 
MN*HN*F*)f MP 

: '-31  4,24*,24C 
"  i r  (j-tl  245,245,3 
245  MC*l)8*PG>'T0- THETA*M|NT 
GW**.0 

IF  (PG-PVWI  252,252,250 
25*  GW*f)4*(PG-PVWt/0S0RT(  fGI 
25 2  MNOZ «07*PG/0SQRT (TGI 

rGPRM*  ( HN»TG*(  MC-GW-MFJOZ  I  ♦(  FOR -OH  I  /C PC)  71  GN*MC-GM-*NOZ  • 
PGPRM«TGPRM*( (GN-Gtt-MNOZ l/08*PG/TGI 
If *2  OGN-c.f 
OMN«*.* 

OQN«*.* 

JJ„NN*1 

1310  IF  I J J-D01 1-2 1 1  1003, IOC  3, 1055 
1055  JJ-001U21 
19*3  JJ.JJ-l 
1 1«0 

OO5-OW0GT*OO£ILT« 

(F  ( NN-IOSSFLT/OOF tTN*001 1 1 1  13(0,1004,1004 
13**  IF  ( JJ-NN4OSSFir/00EtTNI  1004,1**4,13*5 
13*5  O05«0O5*(**-JJ»).5l*OOEI.TN/nS$Fl  T 
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nGA*«OP(  1 1,  jj) 

orf«p.oro( it ,jj> 

cnlZZ?0”1 

><>5  0PV'V4P(f)TE*tPJ 
Ott-t.O 

,p  U  73,1^70,1^5 

^T^p-0T0fPi  75,1 0*0 

>75  JF  <  0  >(  (  i  I  *  J  J  f  - 1 ,  I  f  IOC5 1 10*5 

>«5  jf  <oPG{0m5^i5o^i}^**OOELTN/<ORMOL*OGAM,*i,,:,5*oso#r<i.f/orE*»p> 

7'J'  OSOP  T*DS0RT  I  DDO*QGAM*OPG/0TG ) 

07*00/ ?*f)r.AM**2*00/ ( 2. c+ ).6*00i' 

12’5'12«-'»0 

?f5  00Wfr.A.{rOG*,£,rG_0T£Hp)/| 

0S(CMA-00M6GA*005*a l f\ino  1*C) 

00N*0QN»usigma 

ll  S  1  i  •  JJ*ll*0GAH*(  1 ,0-00>**(  1.0/1 
070( (I, >J*ll*OrEMP-OLMRD£/OCPL*Oo’ 
r7  l7Kfl21f ,  1207 

f.7  orom,jj*u«c’TOFp 

IF  n.frO-OXU  I, JJH  1225, 1225, 140" 
ft  IF  ((J0»  1215,1226,1215 
15  IF  (fjPGI  1225, 1225, 122° 

-1  IF  (OXin.JJtJ  14"6,14 >6, 1236 

,F  (ft">-<OOn*(mJHT/f)OELTNM  1410  I4is  ifif 
n  OFPC  T*(OFftCTH-OFRCtLl*<  J -«  i  i  !  *"nci  rii/nrnur, 

GU  ru  1427  . . . OtirN/orOHT|*OF*CTL 

15  OFRC I *GFRC  TL 

‘5  OK  I  1  1  ,  JjFl  )  «nx  |  I  1  ,  J  J  )  ,f)fH  i 

16  Imu.'AlliilT 

9  0F*t  (  I  j  |  *005*rjg 

or,N-.rr,Fi»(jF 

UHN-  CHN»OP*Ut£MP 

IF  MI-11  nr  4,1240,1241 

1  IF  (JJ-ll  1245,1245,1  "7, 3 

5  OMC-OOfl*fjPG/OTG-r)rHFr«*wiwr 
no*c'',o 

IF  inPG-OPVWl  1252, 1 2‘.'>2,  12  y 

OGMOO4*(0P0-UPyttl/0SgR7(f)fC ) 

?  OMN' 7«007*0Pf,/0SttH  noroi 
0Pr,p»*»«orcp«K*j  (0GN-0GW-0HN02  »/npg*opr,/rm;i 
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r|Hfc*TlMe«-00£CTN 
r iMf*s»iflco.r *tjm£ 

PGMMHG»C.35913l*0PGPRM 

WRITE  (6#  102 1  NN,  TIMERS.  PGWMHG , OTGPRM, OMNOZ , OMC , OGN 

NN*NW*l 

nor,«nPGPRW 

OTG-OTGPRM 

IF  <N«-1-2*N)  1302,3'}*VMOr2 

33CC  PG«»»HG*0.359i3l*PGPRW 

TG«V«<OMC*OCPG*OTG^HC*CPG*TGPRW»/{OIC*OCPC*MC*CPG» 

M!N7*C •  J3«W6*0MC*WC*A£  fNT*OELTN*F  KP(-),8*EM NT/PCAL/TGAVI / 1 OM* 

.  thetapvc* 

TMINT»THINT*RINT 

OR«-WINT*OELHR 

OOR«QMC/IOHC*MCI*OR 

F0R»,«C/(O**C»‘'CI*0« 

TIGN-l9.24*RCAt>Ol*TGAV**3/<  I OPGfPGPR* » *A0 I GN*FA 1GN > *<  l.OPGPR"/ 

1  OPGI*(W*CPG*OR/2*OCPG/l'i*nPG/PGPPM|*EKPI 1 . 8*E A IGN/RCAl/TGAV J • 

2  lorc.r 

FRC 1  I G *  l *OFL  TN/T I GN 

$H1GFR«$MIGFRfFACTIG 

WRITE! 6,  U  21  N.TIMCMS.PGWPHr,,  TGPP«,«NnT,Mr  ,GN,TR|NT,1  IGM.TGAV, 

1  SWIGFP 

IF  ISMIGFR-l.ri  264,2*, 28t 
269  N«N*1 

265  IF  IN-I'.lf  27^,28'  ,28J 
2  7'.'  PG«=PGPFW 
T  G*  TGPRW 

GO  TO  2 

28G  H°ITE(6,lf4l  TIMERS 
GO  TO  1 
END 
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APPENDIX  B  -  cont'd 


Propellant  Data  -  Fuel 


4—  - 

M  = 

46  lbm/lb-mole 

cpl 

=  0.694  Btu/lbm"°R 

a  = 

0.5 

cpg 

=  0.  34  Btu/lbm'°R 

A  - 

55.  65  lbm/ cu.  ft. 

TdFP 

r  397.  7  °R 

\  = 
e 

389  Btu/lbm 

k 

=  1.13 

it 

«■< 

486  Btu/lbm 

=  4.91  x  10"6  lbm/ft-sec 

\F  = 

97.  5  Btu/lbm 

Kg 

=  1.9xl0’6 

Btu/sec-ft2-(°R/ft) 

Initial  Conditions  -  Ox 

AtN 

To 

r  (1.  11 

= 

2.  5  x  10'5  sec 

537°R 

5.83  x  10"5  ft 

Calculated  from 

r  (2.  M 

_ 

1  A  .* 

1  *  .  A  a. 

mas  s -median 
i  drop  sire  of  50 

r  {3.  1) 

_ 

3.84  x  10'4  ft 

microns  as 

discussed  n 

Td  (1.1) 

=  Td  (2,  1)  =  Td  (3.  1) 

= 

J 

537°R 

Sec.  HI.  A. 

Chamber  Data  *  Ox 

*  =  0.  120  lbm/ sec 

^yw 

=  2558  ?SF 

Vc  = 

8. 16  x  10’*  cu.  ft. 

t*s 

*  1.6x10* 

*  sec 

A*  = 

9.  61  x  10’*  sq.  ft. 

*max 

=  0.15 

Ac  = 

4.  60  x  10’^  sq.  ft. 

*min 

*  0.02 

Yd  = 

69.0  ftfaec 

=  0. 0852  ft 

Lc  = 

0. 167  ft 

ldh 

r  2.5x10- 

*  sec 

44* 
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APPENDIX  B  -  cont'd 


Propellant  Data  -  Ox 


M 

= 

92  lbm/Ib-moie 

cpl 

0.  36  Btu/lbm-°R 

a 

= 

0.  1 

Cpg  = 

2.  0  Btu/lbm-°R 

r 

92.5  1b  /cu.ft 
m 

TdFP 

471.  8  °R 

^  e 

= 

173  Btu/lb 

m 

k 

1.  08 

Xs 

= 

237  Btu/lb 

m 

Ug 

7.  9  x  10*6lbm/ft-sec 

XF 

- 

68.  5  Btu/lbm 

Kg  = 

2.  4  x  ID"5  Btu/sec-ft 

Reaction 

Intermediate  Data 

9f 

= 

0.  571 

eox 

= 

0.  429 

AHr 

- 

-1619  Btu/lbm 

^int 

- 

1.7  x  10^  cc/molc 

-sec 

^int 

= 

4500  cal/ g-mole 

Ignition 

Data 

(Aj)jgn  ’ 

3.4  x  10!4  cal-cc /mol**'- sec 

^ign 

' 

5200  cal  / g-mole 

Rra  1 

- 

1.987  cal/mole-°K 

(°R/£t) 
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APPENDIX  C.  PRINTOUT  OF  COM!  UTED  RESULTS  FOR 
ENGINE  CONFIGURATION  NO.  1 


The  first  column  of  the  printout  which  follows  gives  the 
number  N  of  time  intervals  through  which  the  solution  has  proceeded. 

At.  any  given  time  into  the  calculation,  second  column,  the  N  for  the 
oxidizer  is  twice  the  N  for  the  fuel  since  the  time  interval  step  used 
for  the  oxidizer  is  one-half  that  of  the  fuel.  The  times  given  in  the 
second  column  do  rot  include  the  time  from  valve  signal  to  initial 
emergence  of  propellant  from  the  injector  into  the  thrust  chamber. 

This  time  was  determined  from  schlieren  movies  and  oscillogram 
traces  to  be  3.  5  msec  for  the  feed  system/propellant  valve/l-on-1 
doublet  injector  assembly  used  in  the  program.  Thus,  3.  5  msec 
must  be  added  to  the  times  given  in  the  second  column  of  the  print¬ 
out  to  give  the  total  time  from  valve  signal. 

The  third  column  of  the  printout  gives  the  partial  pressure 
of  the  appropriate  reactant  and  the  fourth  column  gives  analogously 
the  temperature  of  the  reactant  vapor  in  the  chamber.  The  fifth,  sixth 
and  seventh  columns  give  respectively:  the  mass  efflux  of  vapor  phase 
reactant  through  the  nozzle  during  the  time  interval  in  question,  the 
mass  of  reactant  vapor  present  in  the  chamber  at  the  end  of  the  time 
interval,  and  the  net  mass  vaporized  during  the  time  interval. 

The  calculated  results  are  given  in  the  printout  in  groups 
consisting  of  three  rows.  The  first  two  rows  in  each  three-row  group 
pertain  to  the  oxidizer  for  which  smaller  time  intervals  are  required 
as  a  consequence  of  its  high  volatility.  The  third  row  of  each  three- 
row  group  pertains  to  the  fuel,  MMH.  This  row  also  gives:  (a)  the 
total  mass  of  reaction  intermediate  TMINT  formed  from  the  start  of 
the  calculation  to  the  time  interval  N  in  question,  (b)  the  ignition  delay 
time  TIGN  calculated  for  the  pressure,  temperature  and  vapor  composi¬ 
tion  prevailing  at  the  end  of  t:.ve  time  interval  in  question,  (c)  the 
average  gas  temperature  TGAV  ;or  the  fuel  and  oxidizer  vapors  in  the 
chamber  at  the  end  of  the  time  interval  in  question,  and  (d)  the  jum  of 
ignition  fractions  SMIGFR  from  the  first  time  interval  of  the  calculation 
to  the  time  interval  in  question. 

Additional  calculated  quantities  which  have  been  printed  out 
in  oth;r  computer  runs  (and  which  made  the  printout  voluminous) 
inclua>_.  the  radius,  temperature  and  frozen  fraction  of  each  drop  class 
(3)  in  each  ensemble  (i,  2 . N)  in  each,  or  every  tenth,  time  interval. 
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The  computer  printout  for  the  ignition  delay  calculations 
for  Engine  Configuration  No.  1  is  reproduced  on  the  following  pages. 
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